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Water hammer constant in the direction of flow - m (ft) 

Water hammer constant moving against the flow - m (ft) 

Gas sound speed - m/sec (ft/sec) 

Liquid sound speed - m/sec (ft/sec) 

Constant determined by pipe and screen restraining fixtures 

Pipe diameter - m (in) 

Pipe equivalent wall thickness - m (in) 

2 2 

Elastic modulus - N/m (Ibf/in ) 

Equivalent screen fraction 

2 2 

Gravitational acceleration - m/sec (ft/sec ) 

9. 8067 m/sec^ (32. 174 ft-lbm/lbf-sec^) 

Static pressure in equivalent feet of liquid - m (ft) 

2 2 

Liquid bulk modulus of elasticity - N/m (Ibf/in ) 

2 2 
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Water hammer constant a/gA in direction of flow - sec/m 
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Water hammer constant a/gA against flow - sec/m (sec/ft ) 
Time - seconds 

Temperature, °R or wave travel time - seconds 
Velocity - m/sec (ft/sec) 

Axial coordinate, gas volume fraction 
Ratio of specific heat 
Stress - N/m^ (Ibf/in^) 

Poisson's ratio 

Density - kg/m^ (Ibm/ft^) 
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Traveling against the flow 
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INTRODUCTION 


As the nation progresses toward more routine space transportation opera- 
tions, presently embodied by the development of the Space Shuttle, continuing 
technology improvement is required to increase payload and reduce in-orbit 
costs. As the Shuttle approaches the operational phase, renewed emphasis 
will be placed on improving the existing system by improving subsystems, 
such as using subcritical cryogen storage for life support and fuel cell 
reactant supplies to take advantage of significant weight savings and payload 
improvements. In addition, the Orbital Transfer Vehicle (OTV), as currently 
envisioned, requires orbital fueling, and plans for future manned orbiting 
stations include replenishment of fluid supplies. All of these systerns have 
a common requirement: the ability to acquire and transfer subcritical fluids 
in low gravity. 

Brute-force systems such as engine- accelerated fluid settling and transfer, 
and for some applications, supercritical fluid storage, can be used with 
attendant weight penalties. However, the potential of fine-mesh screen 
devices for low-gravity propellant acquisition and transfer is starting to be 
realized in some space vehicle systems, including the Shuttle orbital maneu- 
vering system being developed by McDonnell Douglas Astronautics Company 
(MDAC). 

These passive screen acquisition devices allow liquid flow while preventing 
gas flow through the screens by using the surface tension forces in the small 
pores of the screen. The bubble point (the pressure head of vapor which can 
be resisted by the screen) for even the finest-mesh screens is of the order 
of a few hundredths of an atmosphere. Consequently, the screen device must 
be designed for the worst expected combination of gravity head, dynamic head, 
and head losses due to friction. In addition, the device configuration and oper- 
ating characteristics must be selected so the worst head combination imposed 
on the screen is below the bubble point by some defined safety factor. Further, 
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there are a variety of dynamic and transient effects which must be considered 
such as vibration, heat transfer, transient pressure surges at flow startup or 
shutdown caused by valve or pump operation, and surges from liquid boiling 
in warm feedlines. 

In order to design screen devices which will operate reliably for both steady 
and dynamic operating conditions, realistic analyses, based on experimental 
data, are needed. Such analyses are available for steady operation of screen 
devices, accounting for real bubble point, flow losses, non-uniform flow, 
and heat-transfer effects (References 1-4). Analyses are in development for 
vibration effects (Reference 5) and through the program reported herein, for 
transient flow effects on screen retention devices. 

The problem is the ability to analytically predict the effects on screen devices 
of pressure surges caused by the internal environment, configuration, and 
operating characteristics to ensure the successful operation of the unit. This 
requires an evaluation of current and potential screen system designs, and 
an understanding of the critical configurational and operational characteris- 
tics. A complete, versatile analysis of transient feed system dynamics is 
required which can predict pressure and flow transients anywhere within a 
specified feed system and can account for the interaction of components, 
screen device, and fluid. Further, and of considerable importance since the 
transient pressure surges may be several orders of magnitude above the 
screen device bubble point, the analysis must be able to define the pressure 
surge attenuation occurring in lines, fluid, junctions, attenuating components, 
and the screen device. A comprehensive experimental program is required 
which will verify the ability of the analysis to predict screen and feed system 
behavior, or define experimental factors required for successful system 
analysis. 

Limited previous work attempted to define and understand this problem. 

Gluck, et al. (Reference 6) evaluated many of the fluid dynamic effects on 
screen devices of transient pressure surges from valve operation. However, 
this work did not consider the screen device application or potential environ- 
ment, and did not evaluate the potential for complete device destabilization 
(failure). Further, the effects of pressure surges from liquid vaporization 
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in warm lines were not investigated. Although a complex and cumbersom 
compressible flow transient analysis was developed, the program was 
terminated before the experimental data could be used to verify the analysis 
and, therefore, produce a useful design tool. It was determined that 
significant gas ingestion at flow startup and liquid spillover at flow shutdown 
occurs. On the other hand, a limited program in which transient flow effects 
were studied (Reference 7) found no gas ingestion. However, to save cost, 
an existing screen device configuration was used in the latter program without 
evaluation or regard for how well the device simulated potential space 
acquisition systems. Also, pressure surges due to liquid vaporization were 
not evaluated. This work did not result in a comprehensive analytical tool, 
experimentally verified, which could be used for screen system design. 

The objective of the program described herein is to develop a realistic 
design analysis, based on experimental data, to predict the effects of tran- 
sient flow and pressure surges (cau,sed either by valve or pump operation, 
or by boiling of liquids in warm lines) on screen device liquid retention per- 
formance. This will be accomplished. by performing a survey of screen 
liquid acquisition system appiications. and determining the appropriate tran- 
sient pressure and flow environment, analyzing the relations between tran- 
sient flow effects and sc reeb. 'liquid retention characteristics, and experi- 
mentally verifying the analysis. 
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ANALYTICAL PROGRAM 


To ensure that the analysis to be developed would have the widest possible 
application and would reflect the design conditions for realistic systems, a 
comprehensive survey of current and potential space vehicle screen acquisi- 
tion systems was performed. Following this, two existing MDAC computer 
codes, H672- Transient Effects Analysis, and P4557- Cryo-line Pressure 
Surge, were adapted to analyze the screen device effects. These codes 
were then used to simulate a representative surveyed system, and additional 
test cases run to determine the important parameters affecting screen device 
transient performance and aid in developing the experimental program plan. 

SURVEY OF TRANSIENT EFFECTS 

All current and potential space vehicle screen acquisition systems were 
surveyed and documented to determine, where available, data on systems 
involved, fluids, tankagej acquisition system and transfer line configuration, 
valve and pump operating characteristics, acceleration, vibration and thermal 
environment, and fluid quantity and mass flow requirements. A total of 
twelve systems were identified, and their salient characteristics are sum- 
marized in Table 1. Details on the configuration and characteristics of each 
system, where available, are contained in Appendix A. 

The systems include currently developed devices, those in the preliminary 
stages of hardware development, and others which have been studied for 
potential applications, but which don't really exist in the form of hardware. 
The systems tabulated can be conveniently divided into two broad categories: 
A, Small, localized screen devices used for engine restart. These 
devices (e. g. , start baskets, start tanks, traps, etc. ) are char- 
acterized by (1) having relatively large flow rates for a short period, 
and (2) experiencing relatively large g-levels (during which the 
device may fail, then refill). 
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Table 1 

CURRENT AND POTENTIAL USES OF FINE-MESH WOVEN SCREEN IN LIQUID ACQUISITION SYSTEMS 



No. 

Application 

Fluids 

Acquisition System 
C onfigu ration 

Transfer Line 
Configuration 

Valve and Pump 
Operating 
Characteristics 

Ref. 

• 

CRYOGENIC TUG 







1 

Cryogenic APS 

LH^, LO 2 

Multiple channels 

Valves, turbopump, heat 
exchanger, accumulator 

TBD 

4 


2 

Advanced Space Propulsion 
Module (ASPM) Concept 

LH 2 

Numerous screen cylinders 
within rese rvoir with 
containment screen 

Valves, turbopump, 
accumulator 


4 


3 

Centaur (Advanced) 

LH^. LO 2 

Reservoir with containment 
screen 

Subcooler, boost pump, 
engine pump 

RL-10 Pump 

3 


4 

LH 2 Tank Concept 

LH 2 

Annular cylindrical screen 
device in start tank 

Valves, turbopump, 
accumulator 


8 

• 

CRYOGENIC TANKER 







5 

LO 2 Tanker for S-IIB 

LO 2 

Screen liner entire tank 
(5. 54m sphere) 



9 


6 

LH^/LO^ Tanker 

LH 2 , LO 2 

Complete and partial wall 
screen 1 ine r 

Valves, Q/D's, lines 


2 


7 

S-IVC 

LHz 

LO 2 

Reservoir plus channels 
Reservoir with containment 
sc reen 

Valves, turbopump 
Valves, turbopump 


9, 10 
10 

• 

SPACE SHUTTLE 







8 

Space Shuttle OMS 

N^O^. A-50 

4 channels 

Valves, lines 

No Pump 

11 


9 

Space Shuttle RCS 

N 2 O 4 , A-50 

4 channels 

Valves, lines 

No Pump 

12 


10 

Space Shuttle Fuel 
Cell Reactant Supply 

LH 2 , EO 2 

Complete pleated wall 
line r 

Line, regulator 

No Pump 

2 

• 

OTHER 







11 

Improved Agena 

Primary Propulsion Sys 
Secondary Propulsion Sys 

N^O A-50 

N,0^, A-50 
2 4 

Reservoir with several 
containment sc reens 
(common to both systems) 

Valves, lines, turbopump 
Valves, lines, pneumatic 
bellows pump 

8096 Mod Pump 

13 


12 

USAF Unmanned Satellite 

N2H4 

4 channels, 1.57m 
spherical tank 

Valves, lines 

No Pump 

13, 14 



B. Large distributed devices used for extended periods of liquid feed 
for auxiliary propulsion thrusters, life support system and fuel 
cell supply, or fluid transfer. These devices (e. g. , distributed 
channels, screen liners, etc. ) are characterized by ( 1 ) having 
relatively low flow rates for long durations under low-g, and ( 2 ) not 
usually allowing device failure or refill during use. 

The systems tabulated are defined by number as shown in Appendix A. The 
systems in each category are shown in Table 2. In the first category, sys- 
tems 2, 3, and 4 are Cryogenic Tug engine restart systems. System 3, the 
Centaur acquisition system, although conjectural at this time, is fairly well 
defined. Because the L 1 H 2 / LO^ propellants in this system have the potential 
for thermal pressure surges (as well as hydraulic surges), this system was 
initially recommended for further analytical study using the H672 code, as a 
system which is particularly representative of the Cryogenic Tug. System 7, 
for S-IVC engine restart, is also well-defined, but is a more conjectural 
system. However, because it represents a large-scale increase with LH-,/ 
LO^ over the Centaur system, it was also initially recommended for detailed 
analysis using H672. The Improved Agena (No. 11) is also a conjectural 
system which is not well-defined since development was halted in 
mid -program. 


Table 2 

SCREEN ACQUISITION SYSTEM CATEGORIES 


Local! 

zed Engine Restart 

Distributed Long-Term Transfer 

2 

ASPM 

1 

CSS/APS 

*3 

Centaur D- IS 

5 

S-IIB LO^ Tanker 

4 

IDU 

6 

Tug Transfe r Module 

*7 

S-IVC Stage 

*8 

Shuttle OMS 

11 

Improved Agena 

*9 

Shuttle RCS 

^Initially recommended for further 

10 

12 

analysis 

Advanced Fuel Cell Supply 
Satellite Orbit Adjust 


7 





In the second category. System 1 is not as well-defined as other APS systems 
in this category and was not recommended for further analysis. Systems 5 
and 6 are poorly defined with respect to valve and pump characteristics and 
downstream transfer line characteristics and could not be analyzed effectively. 
On the other hand. Systems 8 and 9, the Shuttle OMS and RCS systems, are 
very well-defined and differ by an order of magnitude in scale, and thus 
were both initially recommended for further detailed analysis with H672. 

System 10, the Shuttle Fuel Cell Reactant Supply, was of current interest, but 
detailed investigation into the performance requirements for this system indi- 
cated the extreme likelihood that no transient problems could occur. The rea- 
sons for this are that the control valves are actuated when the Shuttle is on the 
launch pad and the tank and acquisition system are full. Further flow demand 
is accommodated through rather slow regulators on demand from rather large 
volume (high compliance) fuel cells. Thus the probability for large transient 
flow surges with this system appears minimal, and this system was not 
recommended for further analysis. 

The final system. No. 12, the satellite orbit adjust system with the 1. 57m 
(62-in) ^ 211 ^ tank, was the only currently operational screen system found. 
However, the system is used in a sensitive application and detailed system 
data is limited. Because this system was of the same general order for flow 
rate as the Shuttle RCS, it was not recommended for further analysis. The 
four systems initially recommended for further analysis, using the H672 
code, included both cryogenic and storable, local and distributed, and large 
and small systems, and were thought representative of all applications for 
screen acquisition systems. Howevei; the complexity of modeling of these 
systems, plus the limited utility of the data derived from their analysis, 
resulted in only one system, the Centaur D-IS, actually being modeled, as 
described further below. 

TRANSIENT EFFECTS ANALYSIS 

The analytical tool to be developed to analyze screen device response to 
transient pressure surges must also be able to handle the dynamic simulation 
of the entire feed system This is because the pressure surges generated by 
downstream valves, pumps, or line boiling, will be attenuated as they travel 


8 



back toward the screen by the line material compliance, gas bubbles, if 
any, in the line, bends, constrictions, etc. , and by the screen device itself. 
Our approach, therefore, was to use an existing feed system dynamic analysis 
computer code, modify it to include screen devices, and separately analyze 
the line boiling pressure surge, using it as an input to the dynamic analysis 
code. 

MDAC had developed under its IRAD program a comprehensive computerized 
analysis of transient flow in propellant feed systems. This computer code, 
H672, is described in Appendix B, and is unique in that it uses a simple build- 
ing block technique to model complex engine feed systems by retaining non- 
linear descriptions of valves, tanks, and other engine components at discrete 
junctions in the total system, but solves the flow equations between compo- 
nents by an exact technique, using the method of characteristics. In this 
analysis, a propellant feed system is portrayed as variable sets of compo- 
nents located at discrete junctions within the system. Pressure and flow 
information is transmitted from component to component by characteristic 
waves, generated and modified by perturbations and boundary conditions at 
the components. 

The simulations of feed system components which have been developed include 
feed tanks, valves, surge tanks, spring-loaded accumulators, bellows, 
pumps, and liquid injectors. The component simulations were designed to 
use the characteristic wave equations, together with the specified wave 
modification associated with each component, to determine the pressures 
and flows in the system. The simplicity with which the characteristic equa- 
tions are developed for lines, junctions, aid components is described in 
Appendix B. 

This analysis has been used to analyze transient effects in actual launch 
vehicles, and has been correlated with flight data. For example. Figure 1 
shows the H672 computer simulation of the MDAC Delta fuel pump inlet pres- 
sure transient, which compares well with the flight data. 

In order to analyze screen device dynamic response using H672, it was 
necessary to develop new subroutines for inclusion in the H672 code. The 
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nomenclature employed is that used in H672, as described in Appendix B. The 
H672 code was written in British engineering units, as are the new analyses 
for the screen device performance, in order to be compatible with the existing 
code. The International System of Units (SI) are only used for the final results 
following complete system simulation. Mechanization of the screen device 
model described below consists of the creation of two new subroutines, plus 
the addition of suitable "common blocks" and instructions in the main program. 
The first new subroutine, SCRPRP, is a preparation routine for the actual 
screen system input. The SCRPRP routine takes the screen system input data, 
processes it for correct dimensions, manipulates it as required, and stores it 
for use in SR26. The SR26 routine performs the iterative calculations, and 
stores the solutions and characteristics of each model section of the screen 
device. Development of the SR26 screen model equations is described in de- 
tail in Appendix C. The important relations and features of the SR26 subrou- 
tine are described in the following section. 


CR21 



Figure 1. H672 Simulation of Delta Start Transient 
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SR2 6 SCREEN DEVICE MODEL DEVELOPMENT 

A fine-mesh screen used on a propellant acquisition device exhibits character- 
istic behavior when subjected to a positive or negative pressure pulse in the 
liquid retained within the screen device. With a positive pressure pulse, 
liquid outflow occurs through the screen, which resists flow and also exhibits 
compliance in the direction of the flow through the screen. If the liquid outflow 
through the sci:een possesses sufficient kinetic energy to overcome the 
surface tension forces tending to hold the liquid to the screen, the liquid will 
leave the vicinity of the screen. This may be important to screen behavior 
when a negative pressure pulse arrives later. With a negative pressure 
pulse, and with liquid on both sides of the screen (including liquid which may 
have been pushed through the screen by a previous positive pressure pulse), 
liquid inflow will occur, again resisted by the screen which again exhibits 
compliance. With vapor on the outside of a screen experiencing a negative 
pressure pulse, vapor inflow through the screen will occur when the local 
pressure pulse magnitude exceeds the screen bubble point. Gas will con- 
tinue to flow through the screen, resisted by the compliant screen, until the 
negative pressure pulse magnitude drops below the bubble point (assuming a 
wet screen) or to zero (assuming a dried screen). The pressure pulses are 
propagated along the screen device at a velocity determined by the bulk 
modulus of the liquid (modified by bubble ingestion) and the compliance of the 
screen . 

The first step in the development of the screen-channel model was the deter- 
mination of the equation for acoustic velocity in the screen device which 
included properties of the liquid, gas bubble, channel (supporting structure), 
and screen. It was assumed that initially the screen "pipe" was full of liquid. 
Fine-mesh twilled- weave screens are basically orthotropic in nature, i. e. , 
they have different properties in the right angle directions of the weaving axes. 
As shown in Figure 2, the kind of screen-channel configuration which was 
amenable to analysis was a uniformly distributed screen-channel structure 
which could be characterized as a pipe having properties based on the proper 
combination of screen and channel properties. 
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PIPE PROPERTIES ELASTIC MODULUS, Ep 

POISSON'S RATIO, 

WALL THICKNESS, 6p 

SCREEN PROPERTIES ELASTIC MODULUS, AXIAL DIRECTION, Es,, 
ELASTIC MODULUS, CIRCUMFERENTIAL, Esj 
POISSON'S RATIO, AXIAL, 

SCREEN THICKNESS, Bg 

Figure 2. Orthotropic Structure Nomenclature 


Following the usual derivation of the continuity equation for one-dimension, 
unsteady flow (shown in detail in Appendix C) and accounting for the elasticity 
of the fluid and pipe, and the effect of gas in the pipe, the final equation for 
the acoustic velocity in a liquid-gas mixture in a composite orthotropic pipe 
was; 


a' 



11 

C 



2 


+ X (1 


X) 


( 




K g, 


+ 


(1 - x) 




DCj , 
eE, + e 


1/2 


( 1 ) 


where x is the volume fraction of gas; C and Ct are the acoustic velocities 

g 

in the gas and liquid, respectively; p and p, are the corresponding densities; 

g tj 

P . is the absolute pressure in the pipe; V is the ratio of specific heats for the 


12 



gas; g is the gravitational constant; K is the compressibility of the liquid; D 
and e are the pipe inside diameter and equivalent thickness; and E2 are 
the orthotropic elastic moduli; and and are constants depending on the 
form of pipe constraint. 


The complete screen device was divided into a series of sections, one of 
which is shown in Figure 3. The nomenclature used in Figure 3 is that used 
in the H672 analysis (see Appendix B and "Nomenclature") and the upstream 
conditions of head and flow, HI and Ql, and downstream conditions H3 and 
Q5 are assumed known from the previous timestep so that the waterhammer 
constants can be defined; 


CX = HI + SX . Ql 
CY = H3 + SY ■ Q5 


( 2 ) 



Figure 3. Screen Device Model 
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where 


SX 


a' 
g A 


P 


SY = - 


g A 
c p 


With the waterhammer constants defined, the characteristic identities can be 
evaluated: 


Q2 = (CX - H2)/SX 

Q3 = (CY - H2)/SY (3) 

Q4 = Q2 - Q3 


Some of the Q4 flow, flows through the screen (Q6) because of the head 
difference: 


H2 - H4 


A'Q6 ^ B'(Q6)‘ 


(4) 


where A' and B' are the experimentally determined liquid flow-loss 
coefficients for the screen, is the screen area, and H4 is the tank pressure. 

The rest of the flow, Q4 - Q6, expands the pipe slightly such that the change 
of volume is: 

AT (Q4 - Q6) = L 27tR AR (5) 

The above equations were combined and the outflow Q6 was solved explicitly 
(see Appendix C). The outflow Q6 could actually be liquid outflow, liquid 
inflow, or gas inflow, depending on the relative magnitudes of H2, H4, and 
the proper boundary conditions and flow-loss coefficients. For liquid outflow, 
Q6 is: 
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a = A'/A^ 
s 

b = B’/A ^ 
s 

c = 2e E 2 /a' ( 1 - P^) f 


From Q6, H2 is computed from Equation 4 and the flows from Equation 3. 

Q3 and H2 then become the upstream values used to compute the new con- 
ditions in the next section. 

In order not to permanently deform the screen, the equivalent stress in the 
screen during the pressure pulse must be below some critical stress, e. g. , 
the proportional limit of the screen. Another criterion which was checked 
when gas surrounds the screen device was whether liquid expelled through 
the screen by a pressure pulse would leave the screen. 

Assuming negligible potential energy, the criterion for liquid breakaway was 
for the globule kinetic energy to exceed the surface tension energy. The 
critical velocity, which resulted in liquid globule breakaway was found 

to be: 



If the liquid velocity was larger than Vlq the liquid globule would leave the 
screen. In addition, if the quantity of liquid outflow was sufficient to wet the 
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entire screen surface, it was assunned that screen pore surface tension 
forces no longer existed, so that the liquid also left the screen. The quantity 
of liquid (QLSUM) held by the screen is kept current by correction (if any) 
in each time step. If QLSUM is positive, then inflow during a negative pres- 
sure pulse is liquid, until QLSUM goes to zero. This may retard gas ingestion 
for several time steps. 

The model for gas ingestion and quality assumes that the quality in the device 
is uniform (i. e., the gas bubbles are not concentrated near the screen, but 
are distributed uniformly throughout the liquid in the screen). 

The quality in the device is equal to the gas volume divided by the liquid 
volume (screen device volume minus gas volume. The gas volume (QSUM) 
is corrected in each time step by the amount of entrained gas that outflows 
(or inflows) from the screen device. The acoustic velocity a', from 
Equation 1 is updated in each time step by the quality from the preceding time 
step. 

The most important boundary condition for the screen segment is whether 
gas or liquid surrounds it — this condition is input for the segment. The 
basic SR26 routine has three main branches as shown in the logic flow 
diagram in Figure 4. The code first checks if gas or liquid surrounds the 
screen — if liquid, it goes to Branch 1 which determines liquid inflow or out- 
flow, If gas surrounds the screen, the code checks the previous time- step 
screen pressure-ullage pressure difference. If positive (screen pressure 
higher than ullage pressure) it goes to Branch 1 for liquid outflow. If nega- 
tive, it checks whether the pressure difference exceeds the bubble point; if 
so, it goes to Branch 2, which determines gas ingestion. If the pressure 
difference is less than the bubble point, it goes to Branch 3, which sets the 
inflow (Q6) to zero. 

In Branches 1 and 2, the screen stress based on the device section pressure 
difference is computed and compared to the maximum allowable stress 
(input), and if excessive, a warning message is printed. In addition, for 
liquid inflow /outflow the QLSUM is computed in Branch 1 (and Branches 2 
and 3, if applicable). In all branches, the quality within the device section 
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Figure 4. SR26 Logic Flow Diagram 
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is computed in each time step; if the quality exceeds 1. 0, it is set to 1. 0, and 
a warning message is printed. 


The final configuration of the SCRPRP subroutine is shown in Appendix D., 
The data contained in SCRPRP is shown in Table 3. Parameters HI through 
N must be input for each screen device section; SURX (17) through (ZO) are 
computed within the SCRPRP subroutine. 

Table 3 

DATA CONTAINED IN SURX BLOCK 


HI 

Screen Device Location 

SURX(l) 

GAMMA 

ullage gas y 

SURX(2) 

RTGAS 

ullage gas R- T (Ft/*R-‘R) 

SURX(3) 

PGAS 

initial ullage pressure (Psia) 

SURX(4) 

BPS 

screen bubble point (Ft) 

SURX(5) 

EE 

screen/pipe (eE^/D^) (Lb/In"^) 

SURX(6) 

STEM 

maximum stress in screen (Psi)* equiv. screen thick- 
ness (In) 

SURX(7) 

ASl 

screen flow-loss coefficient laminar- liquid (Sec) 

SURX(8) 

BSl 

2 

screen flow-loss coefficient turbulent-both (Sec /Ft) 

SURX(9) 

AS2 

screen flow-loss coefficient laminar-gas (Sec) 

SURX(IO) 

DPS 

screen pore size (Ft) 

SURX(ll) 

DS 

screen channel diameter (In) 
/DCj DC^X 2 

SURX(IZ) 

DC 

screen/pipe 1 (In /Lb) 

SURX(13) 

P0 

screen fraction of open area 

SURX( 14) 

F 

screen/pipe screen fraction 

SURX(15) 

LS 

screen/pipe length (In) 

SURX(16) 

PW 

boundary code for screen/pipe 

0. 0 if dry 

1. 0 if wet 

N 


Propellant type 1 or 2 

SURX(17) 

RH0G 

P/RT 

SURX(18) 

CG 

(ygRT)^/^ 

SURX(19) 


P^/y'Pg 

SURX(ZO) 


Pc/Kg 


tB 





In order to integrate the screen device section(s) into, the total H672 feed 
system simulations, coded control numbers for head, flow, etc. must be 
input for each screen device section. These are called "J- BLOCK'' and are 
shown in Table 4, where the nomenclature refers to that in Figure 3. 

Table 4 

J- BLOCK FOR SCREEN DEVICE 

J- Block 
Location 

1 SPTY Screen Device Identification Code (1 through 4) 

2 HI Pressure Head Control Number for station 1 

3 H2 Control Number for pressure head at screen device junc- 

tion, station 2 

4 H3 Pressure head Control Number for station 3 

5 Q1 Flow Control Number for station 1 

6 Q2 Control Number for flow at screen device in branch 1-2 

7 Q3 Control Number for flow at screen device in branch 2-3 

8 Q5 Flow Control Number for station 3 

9 TT Total Time Control Number for screen device 

10 T1 Control Number for wave travel time between stations 1 

and 2 

11 T2 Control Number for wave travel time between stations 2 

and 3 

12 SX Control Number for pipe constant traveling with the flow 

between stations 1 and 2 (1 through 30 for system 1; 31 
through 60 for system 2) 

13 SY Control Number for pipe constant traveling against the flow 

between stations 2 and 3(1 through 30 for system 1; 31 
through 60 for system 2) 

14 A1 Pipe Type Control Number (1 through 30); the area cor- - 

responding to this pipe type is used to compute the static 
pressure at station 2 

15 PRTY Control Number for propellant type in screen device (1 or 

2); this coincides with the system identification code 

16 DAT Datum Control Number for screen device 

17 INFL Initial Flow Control (1 or 2): 1 signifies that flow initializa- 

tion at the screen device proceeds from station 1; 2 desig- 
nates procession from station 3. 

18 BYP Bypass Flow Control (1 or 2): 1 signifies that the flow 

direction of Q5 follows the assumed convention; 2 designates 
that the direction of Q5 is opposite the assumed convention. 
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Five constants are carried by the H672 program, as shown in Table 5. The 
final configuration of the screen device subroutine SR26 is shown in 
Appendix E. 


Table 5 

H672 PROGRAM CONSTANTS 


CNST 

(1) = 144. 

CNST 

(2) = 576. 

CNST 

(3) = 3.141593 

CNST 

(4) = 32. 174 

CNST 

(5) = 12. ^i^SQRT (62. 4) 


The BLOCK output is presently arranged to output screen device screen- 
ullage pressure differential, quality, local sonic velocity, screen inflow/ 
outflow and QLSUM for up to four screen device sections, together with up 
to 20 heads (pressures) and 10 flows throughout the feed system, and other 
data such as time, valve position, outlet (sink) pressure, etc. 

CRYOGEN BOILING PRESSURE SURGE ANALYSIS 

Pressure surges caused by boiling liquids in warm lines is a significant 
problem which is not specifically modeled in the H672 code but which can be 
easily accommodated by the engine and injector subroutine. This subroutine 
allows a pressure pulse or ignition spike occurring in an engine to be trans- 
mitted to the rest of the feed system through the injector. Our approach was 
to use a zero- resistance injector and model the engine pulse to conform to 
the shape of the pulse generated by boiling fluid, which was predicted from 
the MDAC computer analysis, P4557, described below. 

A number of analyses exist to determine the pressure and flow history which 
occurs with boiling of cold fluid in a warm line during flow, and to concur- 
rently determine the line cooldown time. The most comprehensive was 
developed by Steward, et al. (Reference 15). The MDAC version of that 
analysis, P0734, is large, complex, and expensive to run, and also is much 
too comprehensive for the results needed for this study, which only requires 
the initial (worst) pressure surge and flow reversal. Steward earlier 
developed a much simpler analysis to determine initial pressure surge 
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(Reference 16), and under its IRAD program MDAC developed a computer 
code, P4557, based on this analysis which is simple, inexpensive to run, 
and has been correlated with Steward's experimental LN 2 data. The 
cryopressure surge analysis is described in detail in Appendix F, Important 
features of the analysis are: 

• It accounts for any initial temperature distribution in the line. 

• It accounts for the degree of subcooling of the entering liquid. 

• It determines the pressure and flow history throughout the pipe. 

• Heat-transfer mechanisms are forced convection until boiling 
occurs, then an empirical boiling heat flux. 

• Fanno flow is used to determine the vapor flow characteristics. 

The P4557 program was originally set up for LN^ to correlate with the 

experimental data of Steward (Reference 16), The experimental configuration 

and operating parameters are shown in Figure 5. The major unknown in the 

correlation is the selection of an empirical boiling heat flux. It was found 

for LN, that the Breen- Westwater film boiling correlation shown in Refer- 
^ 2 

ence 17 (3.2 Watt/cm at AT = 220K) gave excellent correlation with the 
experimental data as shown in Figure 6. (The two solid lines represent 
extremes of several experimental tests.) Since use of the film boiling heat 
flux is reasonable from a physical sense, this heat flux, based on the Breen- 
Westwater correlation, was also used for the other fluids of interest to this 
program (LH2~9 Watt/cm^, L02~3 Watt/cm^, Freon 114 -1.5 Watt/cm^), 
together with the appropriate physical properties. The pressure surges for 
these fluids with the same physical system are also shown for comparison in 
Figure 6. The because of its higher density and inertia, peaks at a 

pressure slightly higher, and somewhat later than LN^- Conversely, LH 2 
with its very low viscosity, density, and inertia, peaks very rapidly and at- a 
low pressure value. The Freon 114 was assumed to boil at the minimum 
film boiling heat flux, due to the small difference between the pipe and 
liquid temperatures. The pressure rises slowly toward the saturation pres- 
sure at the pipe wall temperature (26.2 N/cm [38 psig]). 

The P4557 code was used to predict the boiling pressure surge for both 
experiment planning and for experiment data correlation, as described in 
the following paragraphs. 
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FLUID, LNj 

FLUID TEMPEFIATURE, 75.7°K 

TANK (DRIVING! PRESSURE, 51.6 N/CM^ (74.8 PSIAl 

AMBIENT PRESSURE, 8.5 N/CM^ (12.3 PSIA) 

VALVE OPENING TIME, 0.5 SEC 
EQUIVALENT COLD PIPE LENGTH, 127 CM 
DISTANCE TO VALVE, 15 CM 


Figure 5. NBS Experimental Configuration and Operating Parameters 
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Figure 6. Correlation of Transient Pressure Surge 

EXPERIMENTAL PROGRAM ANALYSIS AND PLANNING 

The experimental program analysis and planning was accomplished in several 
steps. First, a simple test model was simulated using H672 (and P4557) to 
check out the SR26 subroutine and also determine the important influences on 
screen device transient behavior. Next, the Centaur D-IS system was 
analyzed in detail, and finally the results of the above simulations were used 
to develop a matrix of test conditions which would evaluate those effects with 
the strongest influence on screen device transient behavior, and provide data 
usable for realistic analytical correlation. 

Two test cases were assumed, with conditions and configuration as shown in 
Figure 7. The fluid for both cases was LH^, saturated at 10.34 N/cm 
(15 psia). The first test case assumed a totally cold system (giving only 
hydrodynamic effects, no cryogenic boiling pressure surge). The system 
shown in Figure 7 consisted of a tank, a screen device, 0. 61 m (2 ft) of 
2 . 5 - cm (1. 0-in) diameter pipe, a valve, 3.05 m (10 ft) of 2 . 5-cm ( 1 . 0-in) 
diameter pipe, and an injector. The tank pressure was assumed constant at 
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© 

SR05 TANK 

FLUID: 

LHj SATURATED AT 15 PSIA 


SR26 SCREEN 


P » 4.4 LB/FT^ 


SR01 JUNCTION 


K - 9750 PSI 

© 

SR11 VALVE 


0 ■» 102X10"* LB/FT 


SR12 INJECTOR 

(SR18 - PRESSURE SURGE) 

ULLAGE GAS: 

GHeAT40®R 
7 = 1.67 


RT - 1546/4 -40 » 15,460 FT 


TANK: 48-IN DIA SPHERE X 0.03 WALL VOLUME •= 33.61 FT3 

ULLAGE HEIGHT, 1 FT. ULLAGE PRESSURE, 25.0 PSIA 

SCREEN: 2.0-IN OIA X 1 .0 FT LONG 326 X 2,300 BPS - 2.42 FT 

AS1 -=1.44 BS1 = 2.27 AS2 = 0.297 DPS = 0.00001 FT 
OS - 2J)IN P<> - 0.060 DC » 0.0936 IN^/LB LS - IX) FT 
EE = 2.46LB/IN^ STEM ■= 15.75 PSI IN F - 0.999 PW - 0.0 

PLUMBING: PIPES: 1 .0 IN DIA X 0.02 WALL 

SCREEN: 2.0 IN DIA X 0.00000071 EQUIV. WALL FOR CORRECT a 
E = 30 X 106PSI (STEEL) 

VALVE: 1-IN DIA BALL VALVE UO = 1 AP - 0.003 PSI AT 08 LB/SEC LH? 

Cv ■= 34.4 VALVE OPEN/CLOSE TIME = 0.05 SEC 

INJECTOR: (EQUIVALENT SAFETY FACTOR = 20.) Cy “ 108 


Figure?. H672 Model for Test Cates 
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2 

17. 24 N/cm (25 psia). The pressure downstream of the valve was 
2 

10. 41 N/cm (15. 1 psia). The valve was assumed to be a low-loss 
(Cy = 34.4) ball valve with an opening time of 0. 05 sec. The injector was 
used to simulate a regulator or higher-loss plumbing downstream. The 
injector Cy was chosen so that the resulting flowrate would give a screen 
safety factor of 2. 0, in order to evaluate the performance of a screen device 
designed in accordance with a commonly accepted criterion. The screen 
device was assumed to be an unsupported tube of 325 x 2,300 dutch twill 
screen, 5. 08 cm (2 in) in diameter by 0. 305 m (1 ft) long, surrounded by 
gaseous helium at 22. 2 K (40®R). 

The second test case simulated the pressure surge caused by cryogenic 

boiling in a warm line. The P4557 code was used to determine the pressure 

surge-time history. With only 3. 66 m (12 ft) of pipe, the LH_ flowed freely 

^ 2 

out the 222 K (400*R) pipe with a pressure surge to only 11. 7 N/cm 

(16. 9 psia). In order to simulate the resistance caused by a regulator or 

additional plumbing downstream, an additional 24. 8 m (81. 5 ft) of 2. 5 cm 

(1 in) diameter line was added (giving the same resistance as the assumed 

2 

test case injector), which resulted in a surge to 29.3 N/cm (42. 5 psia) in 
0. 26 sec. A "combustor" with the cryogenic pressure surge trace was added 
to the injector for the second test case (identical to the first case in all other 
respects). 

Hydrodynamics Effects Test Cases 

A series of 10 test cases was run to evaluate the influence of various system 
configurational and operational conditions on the response of the screen device 
and the occurrence of gas ingestion following valve opening (startup) and valve 
closing (shutdown). The test cases, operating conditions, and results are ' 
summarized in Table 6. Details and discussion of each case are found in the 
following paragraphs. 

The first test case was run with the following valve operational mode: the 
valve was started open at 0.001 sec, was wide open at 0.051 sec, was 
started closed at 0. 1 sec, and was completely closed at 0. 15 sec. The case 
run time was limited to 0. 2 sec to reduce computing time while adequately 
evaluating screen performance. The H672 code was run with the output 
plotting option, and the figures that follow were produced by the SC4020 
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Gas Ingestion 


Figure Test Case Startup Shutdown 


8 

Baseline 

Moderate 

Severe 

9 

Screen Properties 

(P0) 

More Severe 

Very Severe 

10 

Screen Properties 
(E) 

Very Severe 

Extreme 

11 

Valve Open/Close 
Time 

Moderate 

Less Severe 

12 

Valve Open/Close 
Time 

More Moderate 

, Much Less 
Severe 

13 

Screen-Valve 

Distance 

Zero 

Ve ry Seve re 

14 

Screen Area Increase 
(Pleating) 

Zero 

Mode rate 

15 

Screen Area Increase 
(Pleating) 

Zero 

Mild 

■f 

16 

Screen Diameter 

Ze ro 

Mode rate 

17 

Screen Flowrate 

Very Mild 

Moderate 


^Pressure pulse reduced-gas quantity ingested increased 


Table 6 





plotter. The screen-ullage pressure differential was chosen as being 

representative of screen performance for plotting purposes, although the 

complete H672 output was available. The baseline test case is shown in 

2 

Figure 8. The tank pressure is 17.24 N/cm (25 psia) and the pressure 

2 

downstream of the valve is 10.41 N/cm (15. 1 psia), or an equivalent head 

differential of 99.7 6 m (327. 3 ft) of (compared with the screen bubble 

point of 0.74 m [2.42 ft] of LH,). During opening, the screen-ullage 

^ 2 

pressure differential is seen to vary between 0 and 0. 07 N/cm (0. 10 psi). 

The sharp vertical lines indicate points of ullage gas ingestion. It can be 

seen that relatively severe gas ingestion occurs during valve opening, but 

much more severe pressure surges and gas ingestion occur following valve 

closure. Figures 9 and 10 show the effect of screen properties on device 

performance. Originally, the percent open area (P0) was input as equal to 

the void fraction (0. 245); this may not have been accurate. The percent open 

area is not analytically definable for dutch twill weave screens. For square- 

weave screens, the open area is analytically defined, and generally is 

approximately equal to the void fraction squared. This may not be true for 

dutch twill weave screens, but use of the void fraction squared is probably 

closer to reality than just the void fraction. Therefore, for parametric 

purposes, the P0 was input as 0. 06. Figure 9 shows the baseline case but 

with P0 = 0. 245, which indicates a rather strong effect on gas ingestion — 

in fact, the quality in the screen device after 0. 2 sec was 25.43% for 

P0 = 0. 245 compared to 9.49% for P0 = 0. 06. Similarly, Figure 10 shows 

the performance of the baseline case with the screen modulus of elasticity 

2 

equal to twice the previously assumed value of approximately 3, 860 N/cm 
(5, 600 psi). Gas ingestion is so severe that the screen device is completely 
full of gas (quality equals 1.0) by 0. 188 sec. Clearly, screen properties can 
have a strong effect on screen device performance and must be accurately . 
accounted for. 

Considering the relatively slow valve operation time (0. 050 sec), the sever- 
ity of pressure pulse following valve closure was unexpected. Examination 
of the flow characteristics of the system revealed a very interesting effect. 
The valve does not control the flow rate in the system except to turn it "on" 
or "off". The flow rate is controlled by the higher resistance plumbing 
(injector) downstream of the valve. Hence, when the valve is only 4% open, 
it is at essentially full flow. This means that the effective valve opening and 
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closing time is not 0. 050 sec but about 0. 002 sec, which is very fast indeed. 
This is a real system effect and will occur whenever the primary flow rate 
control is not the valve, but rather the downstream system, be it a long 
transfer line and plumbing in the case of a propellant transfer system, or an 
engine injector in the case of an engine restart system. Thus, it is the 
effective valve operation time, relative to the total system, rather than the 
gross actuation time, with which the designer must concern himself. 

In Figure 11 the effective valve opening /closing time was doubled. There is 
essentially no effect on screen performance during valve opening, but the 
severity of the shutdown pulses is lessened considerably. Similarly, in 
Figure 12, the effective valve opening/closing time was increased by a 
factor of 5. It can be seen that the screen performance during valve opening 
was improved, but some gas ingestion still occurred. At shutdown, the 
pressure spike was considerably attenuated, but gas ingestion still occurred. 

It was thought that increasing the length of pipe between the screen and valve 
should have the same beneficial effect on screen performance as an increase 
in the effective valve opening/closing time. Figure 13 shows the effect of 
increasing this pipe length from 0.61m (2 ft) to 3.05m (10 ft) for the baseline 
case: much reduced valve opening pressure surges. Hence, the effective 
valve opening /closing time relative to the screen device -plumbing configuration 
and line length are important parameters. 

It was felt that increasing the screen area would allow higher inflow/ outflow 
rates, and increased screen deflection, with consequent decrease in the 
screen pressure pulse. Figure 14 shows the performance of a screen device 
with 2 to 1 pleating (which increases screen area without increasing screen- 
device cross-sectional area or safety factor). Clearly, the screen-device 
performance is enhanced, and no gas ingestion occurs during valve opening, 
and in Figure 15, with 5 ,to 1 pleating, the pressure pulses are attenuated 
even more. However, it should be noted that gas ingestion still occurs 
during shutdown with 5 to 1 pleating, and because there is five times the 
screen area, even more gas volume is ingested than for the 2 to 1 pleated 
screen (or the baseline screen) despite the reduction in pressure pulse 
intensity. This assumes that gas ingestion occurs over the entire screen 
surface, when in fact experience indicates that there are preferential bubble 
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LH2 SCaCEN DEVICE 


Figure 11. Screen Device Performance, 2 X Effective Valve Open/Close Time 
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Figure 15. Screen Device Performance, 5/1 Pleated Screen Area 
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breakthrough locations (missing wires, pinholes, etc. ) on the screen. How- 
ever, the total gas ingestion required to satisfy the system flow and head 
balances should not be affected by whether ingestion occurs at a few points 
or over the entire screen surface. Because of this area effect, pleating 
is most effective when used in conjunction with other attenuation methods 
so that the pressure pulse is reduced to a value less than the bubble point and 
zero gas ingestion occurs. 

Another method of providing increased screen area is to increase the 
diameter of the screen device, but this decreases the velocity and dynamic 
head in the device, thus increasing the screen device safety factor. Figure 16 
shows the performance of the screen device with the diameter increased by 
the fourth root of 2, so that the safety factor is doubled to 4. Gas ingestion 
is eliminated during valve opening, and much reduced during shutdown. 
Another method of increasing safety factor, and one over which the designer 
may have little control, is to reduce the flow rate. Figure 17 shows the 
performance of the screen device with the flow reduced by the square root of 
2, so that the safety factor is again doubled to 4. Note that, although gas 
ingestion is reduced, flow-rate reduction is not as effective as diameter 
increase because no increase in screen area is obtained. 

Cryogenic Boiling Test Cases 

The results for the second test case for cryogenic boiling pressure surge 
are shown in Figure 18. The magnitude of the pressure surge (to a peak 
of 29. 3 N/cm^ (42. 5 psia) in 0. 26 sec) is shown on the top of the figure. 

The screen-ullage pressure differential is shown on the bottom of the figure. 
The initial results for the first 0. 06 sec (during valve opening) are essentially 
identical to those shown in Figure 8. It can be seen that the screen can 
easily accommodate the reverse flow caused by the downstream pressure 
surge with no problems of excessive stress or bubble ingestion. 

It was thought that LO^, with its higher surge pressure, might be a more 
severe case for screen device performance. However, although the pressure 
surged to 75. 66 N/cm^ (109.73 psia) in 4. 11 sec, no gas ingestion occurred 
during startup, shutdown, or during the cryogenic pressure surge, because 
of the higher surface tension of l-iO^ (compared to BH^). Therefore, it was 
judged that BH2 propellant represented a more severe test of screen device 
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Figure 16. Screen Device Performance, Increased Screen Device Diameter 
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Figure 17. Screen Device Performance, Decreased Screen Device Flowrate 
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Figure 18 . Cryogenic Pressure Surge - Screen Pressure Transient 

performance. The screen device responded to the cryo-surge in the 

same manner as for the — it out flowed easily during the pressure surge 

backflow and the maximum screen-ullage pressure differential reached was 
only 0.0169 N/cm^ (0. 0245 psi). 


These preliminary results implied that cryogenic pressure surges may not 
be as serious a problem as other dynamic effects. This would be investigated 
further in the experimental program. 


Centaur D-IS Test Cases 

The Centaur D-IS Acquisition System (Reference 3) was analyzed using the 
P4557 Cryosurge Code and the H672 Transient Analysis. A significant 
increase in the complexity of the model was required to simulate a screen 
device surrounded by both liquid and gas, the subcooler, and the branched 
lines to the two RL-IO engines. The model and the required parameters are 
shown in Appendix G for both the LH 2 and L<02 feed systems. The P4557 
Code was used to determine the cryogenic pressure surge for the two pro- 
pellants. It was assumed that the worst case would be the engine chilldown 
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and startup with the sump and boostpump chilled and filled. The sump/pump 

chilldown was at very low flowrates (0. 159 kg/sec (0. 35 Ib/sec) for and 

0.0658 kg/sec (0. 145 Ib/sec) for L'02 while the engine chilldown and startup 

rates were much higher: 5. 08 kg/sec (11. 2 Ib/sec) of and 25. 58 kg/sec 

(56.4 Ib/sec) of The temperatures of the feedlines were specified in 

Reference 3 as 144. 4 K (260°R) for and 150 K (270°R) for LO^. The 

engine chilldown valve opening times were assumed at 0.05 sec (Reference 3). 

The engine feedline diameters were assumed to be 6. 35 cm (2. 5 in) for both 

the and systems (Appendix A) and in order to obtain the necessary 

restriction caused by the LO^ valve /injection and the chilldown valve, 

these large-diameter lines were made long to obtain the proper flowrate. 

2 

For the LH, side, the pressure peaked at 14.62 N/cm (21.2 psia) in 0.08 sec. 

^ 2 
compared with the tank driving pressure of 14. 13 N/cm (20. 5 psia). This 

pressure surge is so small and so close to the driving pressure that it was 

judged to have a negligible effect on the screen dynamic performance (recall 

2 

that the H672 test case surge of 29. 3 N/cm (42. 5 psia) relative to the 
2 

17.24 N/cm (25 psia) driving pressure had no adverse dynamic effect, as 
discussed previously). The LO 2 pressure surge case was complicated by 
limitations of the model: the higher-density LO 2 surged much farther into 

the artificially long pipe than the actual distance. The pressure surged to 

2 2 
18. 27 N/cm (26. 5 psia) — relative to the 21. 72 N/cm (31. 5 psia) driving 

pressure — at 0. 13 sec, subsided, again reached 18. 27 N/cm^ at 3. 27 sec, 

and continued to climb due to the increased heat transfer from the artificially 

long pipe surrounding the L 1 O 2 . Accordingly, the heat transfer was reduced 

to approximately the correct value by correcting the boiling heat fl\ix and 

2 

pipe wall temperature. This resulted in a maximum pressure of 18.75 N/cm 

(27. 2 psia) in 0.74 sec, which then subsided and did not reach the driving 

2 

pressure of 21,72 N/cm (31.5 psia) until 41 sec. This was longer than the 
engine start transient time (Reference 3); thus it was judged that cryogenic 
pressure surge was not a significant problem for the Centaur startup, and 
hence it was ignored in the Centaur dynamic simulation. 

The simulation setup shown in Appendix G was run for engine valve startup 
and shutdown, assuming that the half of the screen device surrounded by 
vapor was next to the tank outlet (probably true for valve opening). The LH 2 
screen device was modeled as a 16. 1-cm (6. 34 in) diameter tube, 2. 44-m 
(8 ft) long with the screen fraction, f, set at 1. 5 to give the proper screen 
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area. The four channels were combined into one 15. 2 cm (6. 0 in) 

diameter tube, 20. 8 cm (8. 2 in) long, with £ equal to 6. 2 to give the proper 

screen area. For the LH- case, the maximum screen-ullage pressure 

^ 2 

difference during opening was 0. 0234 N/cm (0. 034 psi) at 0. 93 sec in the 

2 

screen surrounded by vapor (well below the bubble point of 0. 051 N/cm 

(0. 074 psi), and in the part of the screen surrounded by liquid, only 0. 0022 
2 

N/cm (0. 0032 psi). After valve closing, three pressure pulses were strong 
enough to ingest vapor, with a maximum quality of 1.465%. However, this 
was based on the assumption that the screen was surrounded by vapor at the 
tank outlet; actually, at engine shutdown, the screen will almost certainly be 
covered with liquid at the tank outlet because of thrust settling. Therefore, 
the case was rerun with the screen half covered by liquid at the tank outlet. 
The maximum pressure differential in the screen surrounded by vapor was 
0.000324 N/cm^ (0.00047 psi). 

For the case of LO,, with the screen surrounded by vapor at the tank outlet, 

^ 2 
the maximum pressure difference during opening was 0. 00331 N/cm 

2 

(0. 0048 psi) (compared with a bubble point of 0. 31 N/cm (0,45 psi). At 

2 

shutdown, with the same model, the difference was 0. 029 N/cm (0. 042 psi). 
It therefore appears that the design of the Centaur acquisition system is 
more than adequate for engine startup/ shutdown dynamics. Because the 
complex Centaur system simulation yielded little additional design informa- 
tion, the other system simulations were not performed. 

Test Matrix Development 

The Centaur D-IS simulation and the results of the LH^ test cases were 
reviewed and the relevant parameters and their effects are shown in Table 7. 
In addition to preliminary coupon tests to determine screen properties of 
elastic modulus, Poissons ratio, proportional limit, and a bubble point 
check, the original test matrix proposed in Reference 18 included five basic 
horizontal specimens and two vertical specimens (see Table 2-4 of 
Reference 18), together with two methods of attenuation: an accumulator and 
a diameter- ratio attenuator. Because of the importance of screen/structure 
properties and pleating in their effect on device response, and in concurrence 
•with the NASA project manager, it was decided to eliminate attenuation 
devices and substitute seven additional test specimens to evaluate these 
influences. The test specimens are shown in Table 8, and are representative 
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Table 7 

PARAMETER SENSITIVITY ANALYSIS RESULTS 



Parameter 

Effect on Screen-Ullage 
Pressure Difference 

1. 

Scr een/Structure Properties, Effective E 

Strong 

2. 

Pleating to Increase Screen Area 

Strong 

3. 

Effective Valve Opening Time 

Moderate 

4. 

Effective Valve Closing Time 

Strong 

5. 

Screen-Valve Distance 

Moderate 

6. 

Device Diameter, Increased Safety Factor 

Moderate 

7. 

Device Flowrate, Increased Safety Factor 

Moderate 

8. 

Saturated or Subcooled Fluid 



A. Freon 114 

None 


B. LH^ 

Weak 


Table 8 

TEST SPECIMEN CONFIGURATIONS AND NOMENCLATURE 


Specimen 

i 

Vertical Horizontal 


Screen 

Plain 

Plain 

Pleated 

Coarse Screen 
Backup 

Perforated 
Sheet Backup 
No. 1 

Perforated 
Sheet Backup 
• No. 2 

325 X 2300 

AV 

AH 





200 X 1400 

BV 

BH 

BP 

BS 

BPl 

BP2 

720 X 140 


CH 




- 

165 X 800 


DH 

DP 

DS 

BPl 


500 X 500 


EH 






of the construction methods generally employed in screen devices: plain 
screens, pleated screens, screens backed up with coarse screen (e. g. , 

14 X 14 mesh), and screens backed up with perforated sheet. It was antici- 
pated that the degree of screen preload against the backup perforated sheet 
could be significant, so that two degrees of preload were evaluated, using 
200 X 1400 mesh (see Table 8). 
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The experiment matrix covers five screens of different types. The usual 
fine-mesh dutch twill screens were selected, represented by 325 x 2, 300 
mesh and 200 x 1,400 mesh, plus the plain dutch 165 x 800 mesh. These 
weaves are often used for screen devices and considerable data on their 
properties are available. A reverse dutch weave, 720 x 140 mesh, was 
selected because it was anticipated that it would have a larger elastic modulus 
than Dutch twill screens. The fifth screen selected was 500 x 500 mesh 
twilled square -weave screen which probably would never be used in a screen 
device because it is too flimsy, but which has the advantages of equal 
orthotropic elastic modxili and an analytically predictable open area (which 
will be an advantage for later analytical correlation). All of the screens were 
of stainless steel because fine-mesh screens of aluminum are costly and 
difficult to obtain. The bubble point and flow loss coefficients of all these 
screens have been determined previously (Reference 1) using LH^ at 50 psi. 
The screens employed for the screen/ structure integration evaluation were 
the dutch twill 200 x 1,400 mesh and the plain Dutch 165 x 800 mesh. 

From Table 7, the screen-valve distance had a moderate effect on valve 
opening and closing pressure surges, while the effective valves opening/ 
closing time had a moderate to strong effect. Therefore, these were selected 
as additional parameters in the test matrix. Since the device diameter is 
fixed, the.most convenient way to vary the safety factor is to vary the 
flowrate, again a matrix parameter. The initial value of flowrate will be set 
to give a safety factor of 2, the next value will be double or halve the safety 
factor, depending on the results of the first test. 

Selection of the proper test fluids was an important consideration for the 
experimental program. Isopropyl alcohol, an inert fluid with known charac- 
teristics and which allows good visibility, was selected to evaluate the fluid 
dynamic transient effects. MDAC has had considerable experience in testing 
screens with isopropyl alcohol (References 1, 4, and 8) and has found it an 
excellent bubble point simulant for To determine the combined fluid 

dynamic and boiling pressure surge transient effects, a simulant fluid with a 
lower boiling point was required. Freon 114, with a normal boiling point 
of 276.7K (38®F), was selected for this fluid. The Freons are characterized 
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by inertness, nonflammability, low toxicity, good transparency, and a wide 
range of fluid properties. The other potentially usable Freon was Freon 21, 
which boils at 282. 2K (48®F). However, because this is a powerful solvent 
which dissolves plexiglass (our proposed test apparatus material, see below), 
it was rejected. Freon 114, while compatible with plexiglass, exhibits many 
other desirable properties, principally that it is an excellent simulant for 
LiO^. It has virtually the identical bubble point, a somewhat higher laminar 
flow loss coefficient (due to the low viscosity of LO^) and only a 24% lower 
speed of sound than selected as the third test fluid to evaluate 

the combined fluid dynamic and thermal pressure transient effects on screens 
using an actual propellant with a high probability of use in space vehicle 
screen acquisition systems. 

It was predicted that with Freon 114, there would be essentially no effect of 
subcooling (or pipe temperature) on the boiling pressure surge because the 
Freon 114 pressure slowly converged on saturation with no pressure spike. 

Two data points were planned to check this nondependence on subcooling and 
pipe temperature using the less attenuated vertical configuration. The balance 
of the tests would be run at saturated conditions in a room-temperature pipe. 

For on the other hand, subcooling was predicted to be more important 

(but still a relatively weak effect due to the rapidity of the pressure spike), 
but the wall temperature was not predicted to be an effective variable. The 
LH 2 tests, therefore, would be run at both saturated and subcooled conditions 
at one wall temperature, except for one point as a check. 

The final recommended test matrix is shown in Table 9 and consists of 178 
tests. It was anticipated that the choice of pipe length and valve open and close 
time for the Freon 114 tests would be influenced by the results of the alcohol 
tests, considering the lower surface tension of Freon 114. For the tests, 

only the short line was planned because of apparatus limitations and the 
requirement to keep up to the valve. Also, because of the rapidity of the 
cryogenic pressure surge, only the relatively fast valve open and close 
time value was planned for the matrix. 
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Table 9 (Page 1 of 2) 
PLANNED TEST MATRDC 

Fluid 


Fluid Specimen 

Line 

Valve T 

Condition 

Flow Rate 

Alcohol AH 

S 

F 


Ql, Q2 

92 tests AH 

S 

SL 


Ql. Q2 

AH 

L 

F 


Ql, Q2 

AH 

L 

SL 


Ql, Q2 

AV 

. L 

F 


Ql, Q2 

AV 

L 

SL 


Ql, Q2 

BH 

S 

F 


Q3, Q4 

BH 

S 

SL 


Q3, Q4 

BH 

L 

F 


Q3, Q4 

BH 

L 

SL 


Q3, Q4 

BV 

L 

F 


Q3, Q4 

BV 

L 

SL 


Q3, Q4 

BP 

S 

F 


Q3, Q4 

BP 

S 

SL 


Q3, Q4 

BP 

L 

F 


Q3, Q4 

BS 

S 

F 


Q3, Q4 

BS 

S 

SL 


Q3, Q4 

BS 

L 

F 


Q3. Q4 

BS 

L 

SL 


Q3, Q4 

BPl 

S 

F 


Q3, Q4 

BPl 

S 

SL 


Q3, Q4 

BPl 

L 

F 


Q3, Q4 

BPl 

L 

SL 


Q3. Q4 

BP2 

S 

f/sl1 


Q3, Q4 

BP 2 

L 

F/SLl 


Q3, Q4 

CH 

S 

F 


Q5, Q6 

CH 

S 

SL 


Q5, Q6 

CH 

L 

F 


Q5, Q6 

CH 

L 

SL 


Q5, Q6 

DH 

S 

F 


Q7, Q8 

DH 

S 

SL 


Q7, Q8 

DH 

L 

F 


Q7, Q8 

DH 

L 

SL 


Q7, Q8 

DP 

S 

F 


Q7, Q8 

DP 

S 

SL 


Q7. Q8 

DP 

L 

F 


Q7, Q8 

DS 


F 


Q7, Q8 

DS 

S/L^ 

SL/f2 


Q7, Q8 

DS 

L 

F/SL2 


Q7, Q8 

DPI 

S 

F 


Q7, Q8 

DPI 

S/l3 

SL/F3 


Q7, Q8 

DPI 

L 

F/SL3 


Q7, Q8 

EH 

S 

F 


Q9. QIO 

EH 

S 

SL 


Q9. QIO 

EH 

L 

F 


Q9. QIO 

EH 

^Extreme case of BPl 
, Extreme Case of BS 
Extreme Case of DS/BPI 
Specimen See Table 8 

L 

SL 


Q9, QIO 


Line Valve T 

S = Short F = Fast 

L = Long SL = Slow 
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Table 9 (Page 2 of 2) 
PLANNED TEST MATRIX 


Fluid 

Specimen 

Line 

Valve T 

Fluid 

Condition 

Flow Rate 

Freon 114 

AH 

S 

F 

SAT, T1 

Qll, Q12 

60 tests 

AH 

S/L4 

SL/F4 

SAT, T1 

Qll, Q12 


AH 

L 

F/SL4 

SAT, T1 

Qll, Q12 


AV 

L 

F 

SAT, T1 

Qll, Q12 


AV 

L 

F 

SUB, T1 

Qll, Q12 


AV 

L 

SL 

SAT, T1 

Qll, Q12 


AV 

L 

SL 

SUB, T2 

Qll, Q12 


BH 

S 

F 

SAT, T1 

Q13, Q14 


BH 

S/l4 

SL/f4 

SAT, T1 

Q13, Q14 


BH 

L 

F/SL4 

SAT, T1 

Q13, Q14 


BV 

L 

F 

SAT, T1 

Q13, Q14 


BV 

L 

SL 

SAT. T1 

Q13, Q14 


BP 

S 

F 

SAT, T1 

Q13, Q14 


BP 

S 

SL 

SAT, T1 

Q13, Q14 


BP 

L 

F 

SAT, T1 

Q13, Q14 


BS 

S 

F 

SAT, T1 

Q13, Q14 


BS 

S/l4 

SL/F4 

SAT, T1 

Q13, Q14 


BS 

L 

F/SL4 

SAT, T1 

Q13, Q14 


BPl 

S 

F 

SAT, T1 

Q13, Q14 


BPl 

S/l4 

SL/F4 

SAT, T1 

Q13, Q14 


BP1/BP25 

L 

F/SL4 

SAT, T1 

Q13, Q14 


CH 

S 

F 

SAT, T1 

Q15, Q16 


CH 

S/L4 

SL/F4 

SAT, TI 

Q15, Q16 


CH 

L 

F/SL4 

SAT, Tl 

Q15, Q16 


DH 

S 

F 

SAT, Tl 

Q17, Q18 


DH 

S/l4 

SL/F4 

SAT, Tl 

Q17, Q18 


DH 

L 

F/SL4 

SAT, Tl 

Q17, Q18 


EH 

S 

F 

SAT, Tl 

Q19, Q20 


EH 

S/l4 

SL/f4 

SAT, Tl 

Q19, Q20 


EH 

L 

F/SL4 

SAT, Tl • 

Q19, Q20 

LH2 

AH 

S 

F . 

SAT, Tl 

Q21, Q22 

26 tests 

AH 

s 

F 

SUB, Tl 

Q21, Q22 


AH 

s 

F 

SUB, T2 

Q21, Q22 


AV 

s 

F 

SAT, Tl 

Q21, Q22 


AV 

s 

F 

SUB, Tl 

Q21, Q22 


BH 

s 

F 

SAT, Tl 

Q23, Q24 ' 


BH 

s 

F 

SUB, Tl 

Q23, Q24 


BP 

s 

F 

SAT, Tl 

Q23, Q24 


BP 

s 

F 

SUB, Tl 

Q23, Q24 


BS 

s 

F 

SAT, Tl 

Q23, Q24 


BS 

s 

F 

SUB, Tl 

Q23, Q24 


BV 

s 

F 

SAT, Tl 

Q23, Q24 


BV 

s 

F 

SUB, Tl 

Q23, Q24 

4 

^Extreme based on alcohol data 




If BP2 effect from alcohol data 




Specimen Se 

e Table 8 





Line 


Valve T 

Fluid Condition 

S = Short 


F = 

Fast 

SAT 

= Saturated 

L = Long 


SL = 

Slow 

SUB 

= Subcooled 
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EXPERIMENTAL PROGRAM 


The experiment planning and analysis described in the previous section 
indicated that the screen structural properties have a profound effect on 
screen device response to transient pressure surges. The only data on 
screen structural properties (see Reference 18) are of obscure origin, and 
therefore the first experimental task was to determine the structural prop- 
erties of the screens to be used. Following this, the test apparatus, includ- 
ing screen specimen configuration, instrumentation, components, and 
apparatus arrangement, was designed and analyzed. Exploratory tests were 
performed to deterrnine valve operating characteristics and instrumentation 
response and define required line length parameters, and the test apparatus 
was fabricated and installed. Test operational philosophy and procedures 
were developed, and the test matrix for all three test fluids was performed. 

The test results are discussed in this section, but the data correlation using 
the H672 analysis is presented in the section entitled "Data Correlation. " 

SCREEN STRUCTURAL PROPERTIES 

The specimen screen mechanical properties of proportional limit stress and 
effective elastic moduli in both the warp and shute directions were deter- 
mined using an Instron tensile testing machine. The stress was determined 
using the load and the actual screen wire cross-sectional area, and the 
elastic modulus by dividing by the actual strain as determined by a 5. 08-cm 
(2 in) Instron strain gage. The screen specimens were generally about 2. 3 
by 25 cm. For all five specimen screens, the load-strain curve was deter- 
mined in both the warp and shute directions. For the 14 x 14 backup mesh, 
the curve was determined only in the shute direction (the material was too 
narrow to provide a specimen of sufficient length in the warp direction). 

For dutch weave screens, the shute direction is parallel to the shute wires, 
which are the fine wires bent over and under the warp wires. The warp 
direction is parallel to the warp wires, which are the heavier wires laid 
straight along the length of the screen fabric. Exceptions are the 
reverse dutch, 720 x 140, which has heavier shute wires that are 

PRECEDING PAGE BLANX NOT flLME? 
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essentially straight, and finer warp wires bent over the shute wires, and 
the square weave, where both warp and shute wires are bent over one another. 
For the dutch weave screens, the load- strain curves were of different shapes 
in the shute and warp directions. In the shute wire direction, shown typi- 
cally in Figure 19» the initial region of the curve was straight, leading to a 
definite inflection point. It is believed that the springiness of the bent shute 
wires was responsible for this initially straight section of the curve. Com- 
plete strain-free recovery was observed along this initial section of the 
curve. At the inflection point, the wire springiness apparently disappeared, 
as the wires reached full (springy) extension, and strain deformation of the 
wires appeared. From this point on, strain offset appeared as shown by the 
reversed load, and the strain increased enormously, due probably to wire 
displacement. In the shute direction, therefore, the proportional limit 
stress was assumed to be at the inflection point shown, and the effective 
elastic modulus as the stress over strain at the inflection point. For the 
warp wires, the typical load-strain curve is shown in Figure 20. Because 
the warp wires are essentially straight, there was no springiness, but only 
simple tension. The proportional limit is the point shown where permanent 
set begins to occur. The screen properties are shown in Table 10. For the 
springy shute wires, the effective elastic modulus was one- sixth to one-third 
of the stainless steel modulus, while for the "straight" warp wires, the 
effective elastic modulus was two-thirds to three-thirds of the stainless 
steel modulus. Why the warp wire modulus was less than the stainless mod- 
ulus is not known; one explanation may be that the warp wires are not really 
straight, but are slightly bent during the weaving process. 

For the twilled square weave 500 x 500 screen, the effective elastic modulus 
was nearly identical in both the warp and shute directions, as expected, but 
the proportional limits were somewhat different. The reason for this is 
unknown, but it may have been an idiosyncrasy of the weaving process; i. e. , 
the warp wires may have been kept tighter during weaving, thus having less 
available spring deflection (and load) before strain deformation occurs. One 
specimen, the 720 x 140 /shute, was stressed to failure, which occurred at 
a stress of 57, 000 N/cm^ (82,700 psi), which is in agreement with the text- 
book value of tensile strength of 304 stainless steel (Reference 19). 
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It is clear from Table 10 that the unsubstantiated data on screen elastic moduli 
shown in Table 2-2 of Reference 18 were completely erroneous. The impli- 
cation of the much higher, values of elastic moduli was that screen devices 
would not attenuate pressure pulses to the degree anticipated during experi- 
ment planning and may be more susceptible to breakdown. 

The maximum attenuation achievable with a cylindrical screen device occurs 
when the minimum value of effective modulus is oriented in the circum- 
ferential direction. Therefore, the screen specimens for testing with dutch 
weave screens were oriented with the shute wires in the circumferential 
direction, the reverse dutch screen with the warp wires in the circumferen- 
tial direction, and the square weave screen with the shute wires in the 
circumferential direction (taking advantage of the higher proportional limit 
stress ). 

TEST APPARATUS DESIGN 

The basic approach to the design of the test apparatus was to design a con- 
figuration which could be used interchangeably in all three test fluids and 
to use the same dewar, components, lines, instrumentation, and data equip- 
ment (where possible) for all three fluids. The advantage of this approach 
is that random configurational variations could be eliminated and an accurate 
determination made of the interaction of screen properties (such as bubble 
point, flow loss characteristics, and moduli of elasticity) and fluid proper- 
ties. The basic test tank selected was a 0. 13 m^ (35 gal) LH 2 dewar of 
34. 5 N/cm^ (50 psig) working pressure which has been used by MDAC in 
many previous screen test programs (References 1, 2, 4, and 8). The 
dewar has five quartz windows which allow lighting and observation of the 
test specimen. Excellent viewing, even in LH 2 * had been experienced in 
the past. 

The instrumentation to be used for data acquisition during the test program 
may impact both the design of the screen specimen and the arrangement of 
the components within the test apparatus; hence, the instrumentation was 
evaluated and selected first. One of the most important requirements of 
the test program was to assure that adequate data would be obtained in order 
to achieve analytical correlation. The basic transient data are pressure. 
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gas-liquid differential pressure, and flowrate. For the isopropanol and 
Freon 114 tests, obtaining adequate transient pressure data was not antici- 
pated to be a problem: available Endevco piezoelectric transducers could be used 
for dynamic pressure measurement, and Statham strain-gage transducers of 
0-0.7 N/cm (1.0 psid) could be used for screen gas-liquid pressure differential 
measurements. These transducers would be close-coupled to the test specimen 
to obtain maximum dynamic response. 

However, these transducers are temperature-limited, and cannot be close- 
coupled (or submerged) in the LH 2 tests. In fact, there are very few pres- 
sure transducers usable at LH2 temperature, and with these the electronics 
are invariably located at ambient temperature. The usual problems with 
these transducers are inaccuracy, zero shift, and requirement for frequent 
recalibration. To circumvent these difficulties, redundant pressure and AP 
transducers could be installed on long sensing lines and situated outside the 
dewar. During the isopropanol and Freon 114 tests, the data from these 
redundant transducers could be compared to the close-coupled transducer 
data to determine the effective lag, and this could be correlated with fluid 
properties. During the LH2 tests, only the remote transducers could be 
used, and the data could be corrected for the lag. This technique has prob- 
lems, however; variations in the predicted lag may occur because of 
variations in the condition of the fluid in the sensing line, and further, remote 
sensing lines in LH 2 are prone to purging and freezing problems. 

Because of the potential problems with LH 2 pressure and AP measurements, 
dynamic flow-measurement techniques were also investigated thoroughly. 

Flow measurements would be required both for steady flow (to regulate the 
flow velocity to values consistent with screen device design requirements) 
and for transient flow surges and reversal in lines and screen device. Since 
the LH 2 tests would be the most severe, it was deemed desirable that the 
same flowmeter should be used for all of the tests. An exhaustive search 
was made to determine the basic kinds of flowmeters suitable for use with 
LH 2 , which also promise to give reasonable dynamic response. There 
were only three basic kinds of flowmeters found: turbine flowmeter with 
magnetic pickup, turbine flowmeter with nonmagnetic pickup, and vortex- 
shedding flowmeter. The turbine flowmeter with magnetic pickup was found 
to be generally available and has been completely developed for use with LH 2 . 
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The turbine is made of magnetic material (ferrite) and has two severe 
shortcomings; the magnetic drag results in a high-velocity sensing thres- 
hold for repeatable data (~35 Hz at 3 m/sec for a 2. 54-cm (1. 0 in) diameter 
meter), and the turbine is rather brittle at LHz temperature and is suscep- 
tible to damage from reverse flow, gas ingestion, and strong flow pulses. 

The vortex- shedding flowmeter has been used with LN 2 , but not with LH 2 . 

The manufacturer sees no reason why it would not be satisfactory with LH 2 . 

To be used with all test fluids, two sets of thermistor pickups must be used: 
one for near-ambient temperatures, and one for LH 2 -temperature. The 
startup vortices do occur with laminar flow, giving a frequency threshold of 
'^22 Hz at 0. 3 m/sec for 2. 54-cm diameter. Reverse flow capability is still 
undeveloped, although the meter is extremely rugged, and would not be 
damaged by reverse flow, gas ingestion, or flow surges. 

The turbine flowmeter with nonmagnetic pickup was found to eliminate the 
problems associated with the magnetic turbine and has been completely 
developed for use in LH2'. The lack of magnetic drag results in high response 
(-'-500 Hz at 0. 12 m/sec for 2. 54-cm diameter) and the meter pickup could 
be configured to explicitly record reverse flow. Because of the high 
response, this capability may not be necessary for our system as reverse 
flow could be implicitly determinable. In addition, the turbine is made of 
high-strength material compatible with LH 2 temperature, and could handle 
reverse flow, gas ingestion, and severe flow pulses without damage. 

Unfortunately, the nonmagnetic-pickup turbine flowmeter was found to be 
costly, and its procurement was not within the scope of the contract. An 
exhaustive search of NASA/DoD facilities to find an existing and available 
unit of this type was unsuccessful. For these reasons, it was decided to 
evaluate the performance of a magnetic-pickup turbine flowmeter in con- 
junction with a high-response Ramapo drag-body flowmeter in simulant 
(water) flow tests. At the same time, the response of remotely located 
(compared to close-coupled) piezoelectric pressure transducers, and the 
operating characteristics of the flow-control valve could be determined. 

The cryogenic valve selected to start and stop the flow during the test series 
was a 2. 54-cm (1. 0 in) diameter Flowmatics ball valve with a solenoid- 
operated pneumatic actuator. This was a clean valve (equivalent L/D of 3) 
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and hence did not offer significant flow resistance. The steady- state flowrate 
was controlled by dewar pressure and by adjusting a CCI cryogenic hand- 
operated globe valve of 2. 54-cm diameter. This valve, wide open, had an 
equivalent L/D of about 340, and thus imparted much of the flow circuit 
resistance. 

The pneumatically actuated cryogenic Flowmatics ball valve was checked 
out with water-flow tests to determine the effective valve open and close 
time as a function of actuation pressure. Large solenoid valves were 
installed on the pneumatic actuator and a 0. 014 m^ (0. 5 ft^) accumulator was 
used to reduce actuation time. The CCI globe valve was installed downstream 
of the ball valve, and the steady- state flowrate was regulated to about 
6. 3 X 10 m /sec (1 gpm). The effective ball valve open/close time was 
determined by measuring the pressure rise and fall time immediately down- 
stream of the ball valve (upstream of the control valve). The results of the 
tests are shown in Table 11, together with the calculated minimum flow 
line length (upstream of the valve) necessary to develop the maximum water- 
hammer pressure surge following valve closure. Isopropanol, with its high 
sonic velocity, would require the longest line at 9. 2 m (30. 2 ft), which 
defines the longest line used with the isopropanol (and Freon 114) test setups. 

During the same tests, the flowmeter and pressure transducer response was 
evaluated. The instrumentation was set up with the Flowmatics cryogenic 
ball valve and CCI flow control valve as shown schematically in Figure 21. 

The tests were run with water at about 31 N/cm^ (45 psi) and with the flow 
control valve set to regulate the steady- state flowrate at about 6. 3 x 10 " 5 

3 

m /sec (1 gpm). The dimensions shown in Figure 21 are in centimeters.. 

The coding for the Endevco piezoelectric pressure transducers is: upstream- 
close (UC), upstream-far (UF), downstream-close (DC), and downstream- 
far (DF). An example of the data for flow startup is shown in Figure 22. 

The bottom curve is the turbine frequency and the top curve is flowrate 
(increasing downward) for the Foxboro magnetic turbine flowmeter. The 
second curve is the flowrate (increasing downward) for the Ramapo drag 
body flowmeter. Note that the Ramapo flow follows exactly (but inversely) 
the pressure traces (increasing upward) for the UC (= Pjjq) and DC (= Pqux^ 
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Table 11 

VALVE OPEN/ CLOSE CHARACTERISTICS 


Maximum Pressure Surge 
Line Length, m (ft) 

Actuation Pressure Effective* Times 


N/cm^ 

(psig) 

(sec) 

Isopropanol 

Freon 114 

LH 2 

51. 7 

(75) 

Open: 0. 01 9 
Close: 0. 019 

11. 6 (38. .2) 

6. 3 (20. 6) 

9.4 (30. 9) 

103, 4 

(150) 

Open: 0. 01 5 






Close: 0. 01 5 

9. 2 (30. 2) 

5. 0 (16. 3) 

7.4 (24.4) 


^ - 5 3 

*Flow Controlled to ~6. 3 x 10 m /sec (1 gpm) 


transducers only (the UF and DF transducers were not used here). This was 
in accordance with the flow characteristics equation; 

AH+^AQ = 0 (8) 

Clearly, the Foxboro flowmeter did not respond to the transient flow/pressure 
pulses. The Endevco piezoelectric pressure transducers indicated the 
change in pressure, and clearly showed the pressure attenuation across the 
Ramapo flowmeter (both pressure transducers had about the same calibration 
factor). The pressure peak-to-peak time constant agreed exactly with the 
wave travel time from pressure transducer to pipe outlet and return. 

Testing with both the close-coupled and remote (UF and DF) pressure trans- 
ducers (see Figure 21) revealed the following not- unexpected results: 

1. Severe pressure surges could be recorded if the long (122 cm) 
sensing lines were not properly bled free from gas. 

2. Small sensing lines (0. 32 cm, 0. 125 in, diameter) attenuated the 
pressure signal and masked details revealed by the close-coupled 
transducers. 

Number 1 above was anticipated to be a potentially severe problem with the 
LH 2 tests, since gas will of necessity be present in the sensing lines. 
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The results of these exploratory tests, as they apply to the instrumentation 
design for the transient test apparatus were as follows: 

1. The turbine flowmeter did not have adequate response for transient 
flow testing; therefore, for the isopropanol and Freon 114 tests, 
the Ramapo flowmeter would be used. For the LH 2 tests, a flow- 
meter would not be used, but pressure and differential pressure 
data would be acquired. 

2. The close-coupled piezoelectric pressure transducers would be 
used for the isopropanol and Freon 114 tests, but remotely located 
Statham pressure and differential pressure transducers would be 
used for the LH 2 tests. The remote sensing lines in the LH 2 tests 
would be 0. 64- cm (0. 25 in) diameter to reduce line attenuation. 

The original screen test specimen configuration proposed in Reference 18 
was envisioned as a cylindrical plexiglass channel with screen bonded across 
an opening along the top. This configuration was reviewed and modified to 
provide directly correlatable data. Use of a partially cylindrical plexiglass 
channel with flat screen bonded on top would probably lead to correlation 
problems. Such a configuration may be analytically difficult to characterize, 
since the channel properties may dominate the screen/structure prop- 
erties (see analysis section), with the result that the differences in response 
for different screen/structure configurations may be masked. The 
plexiglass channel would have the advantage of viewability, and gas ingestion 
may be able to be directly observed (although the quantity would still have 
to be determined from the bubble trap). It appeared likely that the partial- 
screen channel would have lower damping than a full-screen channel, so 
that a strong pressure pulse could crack it, especially in L.H 2 . Further, 
the channel would have a different coefficient of expansion than the screen,' 
so that the flat- screen tension would change from the simulant fluid tests to 
the LH 2 tests, leading to additional correlation difficulties. 

On the other hand, a cylindrical screen specimen would be a more uniform 
structure which could be analyzed directly. This type of specimen would 
provide direct correlation of the interaction of screen and structural backup 
and their effect on specimen response. Because the specimen would all be 
of the same material, contraction effects would not change screen tension. 

It was probable that the occurrence of gas ingestion with this specimen. 
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while not directly observable, could be inferred fronn the time of pulse travel 
through the specimen. Further, at shutdown, the amount of liquid spill, 
which would be observable, could probably be related directly to gas 
ingestion. 

It was, however, believed to be more important to obtain response data in 
support of screen/structure configuration analysis, than to simply observe 
gas ingestion. Therefore, the recommended design for all specimens, and 
especially those with structural backup, was that the screens be configured 
as cylinders, bonded to an open tubular support frame, as shown in 
Figure 23. The fluid inlet would be through the screen on the right, and 
the gas pressure would be imposed on the screen on the left, surrounded by 
a plexiglass enclosure. The close-coupled bubble trap would also be made 
of plexiglass for direct observation of gas ingestion quantity during startup. 


The plexiglass walls would be thick enough to withstand the predicted pres- 
sure surges in isopropyl alcohol and Freon 114, but there was some question 
as to the use of plexiglass with LH 2 . MDAC has run tests in LH 2 using 



Figure 23. Screen/Structure Specimen Test Configuration 
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plexiglass apparatus (References 1 and 2) where it was discovered that if 
the plexiglass was unstressed by point loads, such as steel bolts, fittings, 
etc, and if the dewar was slowly filled with LH 2 to avoid sudden excessive 
chilling, the plexiglass survived without cracking. Both the plexiglass 
bubble trap and the plexiglass gas enclosure would be bonded to the specimen 
in a strain-free fashion, and the large gas enclosure would not be directly 
subjected to a pressure pulse, thus giving hope of survival in the LH 2 tests. 
For each specimen, the screen structure, transducer mounting base, and 
metal end plates would be bonded to the central tube with high-strength LH 2 - 
compatible polyurethane adhesive. Note that the horizontal specimen shown 
could also be directly used for the vertical configuration tests by reversing 
the gas inlet and liquid drain connections. 

The construction of the screen/structure specimens was straightforward: 
the typical plain screen specimen was a simple cylinder of screen, seam- 
welded longitudinally, and sized to slide over the support tube and be bonded 
to it. The coarse mesh backup specimen was made in the same fashion, 
with the coarse mesh welded first, then covered with the fine mesh screen, 
and TIG-welded longitudinally. The coarse mesh backup was sized to slide 
over the support tube, and both screens were bonded to the tube. The 
pleated screen specimen had approximately 3/1 area ratio pleats running 
in the direction of flow. The ends were coined flat and bonded to the tube. 

The two perforated sheet backups were formed into open cylinders, as shown 
in Figure 24, with two values of angular opening, The screen was seam- 

welded into a cylinder, and the perforated sheet cylinder was compressed, 
slipped inside the screen cylinder, and allowed to reopen, exerting tension 
on the screen, as shown in Figure 24. 

The degree of tension exerted will be a function of the initial and final angu- 
lar openings, as shown in Figure 25. The final angular opening was deter- 
mined by examination of the specimen after installation of the perforated 
sheet. Both the screen and perforated sheet were bonded to the support 
tube. The configuration of a typical plain screen specimen, as fabricated, 
is shown in Figure 26. 
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a = TENSILE STRESS AT GAP 
Figure 25. Tensile Stress on Deflected Perforated Sheet 
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Figure 28. Plain Screen Specimen Construction 


TEST APPARATUS INSTALLATION 

The arrangement of the dewar, test specimen, instrumentation, lines, valves, 
and other components is shown schematically in Figure 27. The solid-line 
representation of the on/off valve and throttle valve in Figure 27 shows the 
"short-line" configuration; the dashed-line representation of the valves shows 
the valve position for the "long-line" configuration. The helium pressuriza- 
tion gas used to pressurize both the dewar (for outflow) and the gas enclosure 
around the screen specimen is cooled to fluid temperature with in-tank heat 
exchanger coils as shown in Figure 27. 

The horizontal test apparatus (as installed within the dewar during testing) 
is shown in Figure 28. The dewar is in the background, and the apparatus 
components are identified. The apparatus configuration with the vertical 
specimen attitude is shown in Figure 29. The 325 x 2300 screen was also 
tested in the vertical position without the flowmeter installed in order to 
impose the minimum pressure pulse damping. Figure 30 shows the test 
apparatus arranged with the "short" line length of 1. 2 m (4 ft), and Figure 31 



shows details of the valve/dewar arrangement with the 1.2 m line. The 
outflow line for all tests was 2. 54-cm (1 in) diameter by 0. 089-cm (0. 035 in) 
wall stainless steel tubing. Figure 32 shows the test apparatus arrangement 
with the "long" line length of 10.4 m (34 ft). The arrangement was achieved 
by moving the valve complex from a position adjacent to the tank (see 
Figure 31) to the other end of the 9. 2 m (30 ft) line, next to the supply/ 
receiver tank. The apparatus is shown connected to a drum of isopropyl 
alcohol; the large horizontal cylindrical tank is Freon 114. 

The apparatus configurations shown were used for the isopropyl alcohol and 
Freon 114 tests. For the LH2 tests, the test apparatus was modified by 
removing or remoting instrumentation sensitive to low temperatures, 
arranging the outflow valve to provide for the "short" line and "fast" valve 
operation (see test matrix of Table 8), and providing for apparatus venting, 
purging, vacuum, and line insulation. 
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RAMAPO FLOWMETER 


WIRING J-STRIP 


IN-TANK 

HEAT EXCHANGERS 


OUTFLOW LINE 

(2.54 CM (1 IN) DIAMETER) 


DEWAR WINDOW 
INNER SCREEN 


OUTER SCREEN 


PLEXIGLASS GAS ENCLOSURE 


AP TRANSDUCER 
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Figure 28. Horizontal Test Spedfnen as Installed 
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Figure 30. Short Line Apparatus Configuration 
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Figure 32. Long Line Apparetus Configuration 


TESTING PHILOSOPHY AND PROCEDURES 

For screen acquisition devices used for low-g fluid transfer, the velocity- 
head, flow-through and friction head losses imposed by fluid flow usually 
constitute the bulk of the heads which must be resisted by the screen (since 
the gravitational hydrostatic head is usually very small). Under these con- 
ditions, the fluid flow surges will have the maximum effect on screen reten- 
tion performance. In order to simulate these conditions in 1-g testing, care 
must be taken to minimize the imposed hydrostatic head. The cylindrical 
screen specimen must be surrounded by gas (to ensure that gas ingestion 
or liquid spillover can occur if conditions so warrant) and the surrounding 
gas is imposed by draining the fluid from around the specimen. Complete 
liquid draining was detected by the evolution of small gas bubbles, as shown 
in Figure 33. The maximum gas pressure imposed occurs at the top of the 
cylindrical specimen, and the liquid drain opening was situated (see Fig- 
ure 26) so that the imposed hydrostatic head (Hqr) was minimized to 
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Figure 33. Screen Specimen Head Loss Nomenclature 




5. 1 cm (2. 0 in). The head and flow losses are shown in Figure 33. The 
sum of imposed heads is: 


BP 

S. F. 


^ ^FL ^ ”dY 



where BP is the screen bubble point head capability and A' and B' are the 
screen liquid flow-loss coefficients defined by Equation 4. The safety 
factor (S. F. ) is selected to give an adequate screen device design perform- 
ance margin; some studies (References 1, 2, 4, and 8) have selected a safety 
factor of 2. 0. Kg is the loss coefficient for flow turning and frictional losses. 
The screen approach velocity, Vg, is related to the pipe (specimen) velocity, 
Vp, by the ratio of screen area to cross-sectional area: 



After substitution of Equation 10, Equation 9 was solved for the appropriate 
pipe velocity, Vp, (or baseline outflow rate) which results in a safety factor 
for the screen specimen of 2. 0. The baseline outflow rates were imposed 
through a combination of dewar pressure and control valve setting. As 
defined previously in Table 9, the following test parameters were varied 
during the tests: 

• Fluid; isopropyl alcohol. Freon 114, and LH2. 

• Screen device horizontal (all screens) or vertical (325 x 2300 and 
200 X 1400 plain only). 

Screen weave, ranging from 325 x 2300 dutch twill to 500 x 500 
square 

• j Screen device mounting /support method, ranging from plain and 

pleated screen to coarse screen and perforated sheet backup. 

• Outflow line length from screen to valve of "short" (1.2 m, 4 ft) 
or "long" (10.4 m, 34 ft). 
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• Outflow rate at two baseline values; to give retention safety factors 
of 2. 0 and 1. 0, with minimum possible static gas head of about 

5. 1 cm (2. 0 in) of alcohol imposed. 

• Outflow valve open/close rate of either open/close = 0.019 sec or 
open/close = 0.015 sec. 

The following general procedure was followed during alcohol testing for a 
given specimen and line-length configuration; 

• With the dewar vented, the alcohol storage tank was pressurized 
and the dewar was filled to about 0.2 m (8 in) above the screen 
specimen. 

• The screen specimen static bubble point was determined using a 
water manometer. 

• The storage tank was vented and the dewar was pressurized to an 

2 

appropriate pressure from about 1 to 1 1 N/cm (1.5 to 16 psig) as 
required to obtain the proper steady- state outflow rate (determined 
by stop watch and internal scale). 

• Outflow was terminated while maintaining the proper dewar pressure 
the screen gas head was imposed and the bubble trap bled and filled 
with liquid. 

• With instrumentation and oscillograph operating, outflow was 
initiated (valve opened) and terminated (valve closed) while the 
occurrence of gas ingestion or liquid spillover was observed 
visually through dewar windows. 

• The valve actuation pressure (hence, closing speed) was changed 
and the test repeated. 

• The dewar pressure was changed to achieve the other required 
steady-state outflow rate, and the test repeated, as above. 

Essentially the same procedures were followed for the Freon 114 tests, 
except that after filling the dewar, the Freon was allowed time to saturate at 
atmospheric pressure. (Some early Freon 114 tests were run at saturation 
pressures of 2. 1 N/cm^ (3. 0 psig) and 0. 7 N/cm^ (1. 0 psig); see test results 
which follow. 


72 



Freon 114 tests performed with the long-line configuration revealed that 
continuous backflow occurred due to heat transfer and line boiling, and that 
liquid could not effectively be kept in the long line up to the valve. There- 
fore, most of the Freon 114 tests were run with the short-line configuration 
(see test results). 

The same problem was anticipated with the LH 2 tests, which were all per- 
formed with the short-line configuration, and 0. 015-sec valve actuation time, 

2 

as planned. For the LH 2 tests, the dewar was cyclically purged to 7 N/cm 
(10 psig) about 10 times with helium after specimen changes, before LHg 
filling took place. The L.H 2 was allowed to saturate at either 1.4 N/cm^ or 
0. 0 N/cm^ (2. 0 psig or 0. 0 psig) except for two tests at 0. 7 N/cm^ (1. 0 psig). 
Since there was severe line boiling and surging, it was not practical to 
establish a steady outflow rate; rather, all tests were run with dewar pres- 
sures of about 3. 5 N/cm^ (5 psig) or 6.9 N/cm^ (10 psig) which were 
anticipated to bracket the appropriate flow rate range (see test results). 

TRANSIENT TEST RESULTS 

Isopropyl Alcohol Tests 

As each test specimen was installed for the alcohol tests, the first test made 
was to determine the static alcohol bubble point of the specimen. So that gas 
breakthrough could be observed, the specimen, as installed, was submerged 
in alcohol within the dewar, and alcohol was introduced into the plexiglass 
gas enclosure, surrounding the inner screen. The specimen was then pres- 
surized with gaseous helium on the inside of the screen, and the pressure 
increased until gas bubbles emerged from the inner or outer screen. The 
bubble point of the inner screen was recorded because it was the critical 
parameter for gas ingestion during subsequent tests. The bubble point pres- 
sure was determined with a self-calibrating water manometer, and was 
corrected for the depth of alcohol above the specimen. The data for each 
screen mesh is shown in Table 12 and for comparison previous alcohol 
bubble point data from Reference 1 is also shown. The bubble point pressure 
for each test fluid, as derived from the water manometer data, is shown. 
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STATIC ALCOHOL 


Sc reen 

Corrected 
Alcohol Bubble 
Point 
(in. H^O) 

Reference 1 
Alcohol Bubble 
Point 
(in. H 2 O) 

325 X 2300 

24. 8 

23. 9 

200 X 1400 

16. 5 

16 . 6 

720 X 140 

9. 2 

• 

00 

165 X 800 

7.9 

6. 7 

500 X 500 

7. 7 

7. 7 
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POINT-SCREEN SPECIMENS 

Bubble Point Pressure For Test Fluids 


Alcohol Freon 114 L'H 2 


N/cm^ 

PSI 

N/cm^ 

PSI 

N/cm^ 

PSI 

0.618 

0. 896 

0.411 

0. 596 

0. 052 

0. 076 

0.411 

0. 596 

0. 273 

0. 396 

0. 035 

0. 051 

0. 229 

0. 332 

0. 152 

0.221 

0. 019 

0. 028 

0. 197 

0.285 

0. 130 

0. 189 

0. 017 

0. 024 

0. 192 

0. 278 

0. 128 

0. 185 

0. 017 

0. 024 



For the 200 x 1400 and 165 x 800 mesh screens, there were several 
specimens fabricated. For some of these there was minor leakage from the 
outer screen which made it difficult to obtain accurate bubble point data for 
the inner screen. This leakage was not repaired, since it was thought that 
it would not materially affect the inflow through the outer screen, or the gas- 
resistance capacity of the inner screen. All of the 200 x 1400 (and 165 x 
800) mesh specimens were made from the same piece of screen, and the 
bubble point data which were obtained for different specimens were quite 
close in value; hence the bubble point shown in Table 12 for the 200 x 1400 
and 165 x 800 mesh was assumed as the same for all specimens of that mesh. 

The experiment matrix was set up to explore the sensitivity of the screen 
device response to various parameters as described in Table 7. The screen 
device sensitivity to some of these parameters, such as the effects of screen 
specimen structure (percent open area, screen fraction, effective elastic 
moduli, pleating, perforated sheet preload, device diameter, etc. ) would be 
resolved by correlations with screen model simulations generated by the 
H672 analysis, by comparing gas ingestion, screen-ullage pressure differ- 
ence, etc. in response to given conditions of initial imposed gas pressure 
and flowrate. Screen sensitivity to other parameters, such as valve open 
and close time effects, screen-valve distance, and flowrate, are responses 
to hydrodynamic effects which were evaluated with the alcohol tests (and to 
some degree with the Freon 114 tests). The screen sensitivity to parameters 
such as saturated or subcooled fluid conditions was anticipated to be resolv- 
able by examination of the F reon 114 and LH2 tests and analysis of the cryo- 
surge pressure pulses obtained during testing. 

For each test in the matrix, both startup and shutdown occurred. Screen - 
retention failure or gas ingestion during startup may be of prime concern 
for screen device design since gas could cause problems with downstream 
equipment (e. g. , rocket engines or fuel cells) or, in low-gravity, could lead 
to liquid flow blockage in the screen device and complete retention failure. 

On the other hand, gas ingestion following shutdown may not be as severe a 
design problem, since the mission may have been completed (except when 
multiple flow cycles are required) unless complete screen destabilization 


75 



occurred (the screen device emptied completely, prohibiting further liquid 
flow). For the 177 tests performed, liquid spill (construed as equivalent to 
gas ingestion) after shutdown occurred only twice--once with Freon 114 
(Run 148, Table 13) and once with LH 2 (Run 155, Table 13). In both cases, 
the spill was of the order of I 6 cc ( 1 . 0 in^) and gross gas overpressure had 
inadvertently been imposed (due to circumstances described in detail below). 
Complete screen destabilization never occurred, even under the most severe 
operating conditions imposed during testing, and therefore is believed to be 
an extremely unlikely occurrence. 

The startup performance of the screen specimen was of particular interest 
because ingested gas could be trapped, observed, and measured, and 
because of the implications of startup gas ingestion on orbital screen device 
design. The startup transient and initial flow pulse data for every test are 
shown in Table 13. The table shows the conditions for each test, the initial 
gas -liquid pressure difference applied to the inner screen, the approximate 
steady-state flowrate as recorded by the Ramapo flowmeter, and the time 
and magnitude of the initial flow pulse at startup. Also shown in Table 13 is 
whether gas ingestion was observed, either visually into the bubble trap or 
by flowmeter spiking as recorded on the oscillograph. Figure 34 shows a 
typical oscillograph trace for flow startup with gas ingestion. The spiking 
of the flowmeter trace is apparently caused by springback of the dragbody 
when it encounters a gas bubble. (Trace ID is shown in Figure Z7. ) 

It had been anticipated (see Table 7) that the effective valve open time would 
have a moderate effect on screen startup performance and potential gas 
ingestion. In fact, there was essentially no difference in startup performance 
for the 0. 015-sec effective valve open time or 0. 019-sec effective valve open 
time. This is seen by comparing the oscillograph record (Figure 35) for 
startup of Run 21 (0. 01 9- sec valve, short line) with Figure 36 for startup 
of Run 22 (0. 015-sec valve, short line). The shutdown traces for these 
tests were also identical. This is believed to be because the dynamics of the 
system are dominated by the large flow resistances which control the startup 
transients, rather than the effective open time differences for the low- 
resistance valve. The startup performance for the long line tests (e. g. , 

Run 25; see Figure 37) was also the same for both valve actuation 
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Table 13 


INITIAL FLOW PULSE TEST DATA 


Test 

No. 

Screen 

Configuration 

Line 

Length 

(m) 

Valve 

Opening 

(1) 

' — 

Imposed 
Screen* Ullage 
Pressure Difference 

N/cm^ psi 

Steady Flowrate 
(Approx) 

m^/Sec X 10^ KP *’' 

Time 

After 

Start 

(Sec) 

Flow Pulse 
Screen-Ullage 
Pressure Difference 

N/cm^ psi 

Gas Ingestion 
Observed 

Bubble Oscillograph 
Trap (Flowmeter) 

Isopropyl Alcohol — Horisontal Specimen 











1 

325 X 

2300 

Plain 

1 

2 

F 

0. 26 

0. 36 - 

11. 1 

1.8 

0. 57 

0.41 

0.60 

No 

Yes? 

2 







S 

0.23 

0. 33 

11. 1 

1.8 

0. 74 

0. 40 

0. 58 

No 

Yes? 

3 







F 

0.21 

0. 30 

.10. 7 

1.7 

0.67 

0. 41 

0. 60 

No 

Yes? 

4 





1 

2 

S 

0.21 

0.31 

8.9 

1.4 

1. 1 

0.41 

0. 59 

No 

No 

5 





9 

2 

s 

0.21 

0. 30 

11. 1 

1.8 

0. 33 

0. 37 

0. 54 

No 

Yes? 

6 







F 

0. 17 

0. 24 

12.2 

1.9 

1.75 

0. 51 

0. 74 

No 

No 

7 







F 

0. 08 

0. 11 

10.2 

1.6 

2. 30 

0.42 

0.61 

No 

No 

8 

325 X 

2300 

Plain 

9 

2 

S 

0. 10 

0. 15 

9.3 

1.5 

2. SO 

0.41 

0. 59 

No 

No 

9 

200 X 

1400 

Plain 

9 

2 

F 

0. 17 

0.25 

n . 1 

l.B 

2. )0 

0. 33 

0.46 

No 

No 

10 







S 

0. 17 

0.25 

11. 1 

l.B 

2. 10 

0. 34 

0. 49 

No 

No 

11 







S 

0. 18 

0. 26 

10.2 

1.6 

1.92 

0. 34 

0. 49 

Yes? (2) 

No 

12 





9 

2 

F 

0. 19 

0.27 

10. 2 

1.6 

1.92 

0. 34 

0. 49 

Yes? 

No 

13 





1 

2 

F 

0. 17 

0.25 

17.7 

2. 8 

0. 36 

0. 34 

0. 50 

No? 

Yes? 

14 







S 

0.21 

0. 31 

15. 3 

2.4 

0.90 

0. 39 

0. 56 

No? 

No 

15 







s 

0. 17 

0.25 

12.5 

2. 0 

0.90 

0. 34 

0.49 

Yes 

No 

16 







F 

0. 17 

0.25 

13.2 

2. 1 

0.90 

0. 34 

0. 50 

Yes 

No 

17 







F 

0. 17 

0.25 

15.6 

2. 5 

1. 10 

0. 38 

0. 55 

Yes 

No 

18 

20U X 

1400 

Plain 

1 

2 

S 

0. 17 

0.25 

16. 5 

2.6 

0.90 

0. 39 

0. 56 

Yes 

No 

19 

200 X 

1400 

Pleated 

1 

2 

F 

0. 14 

0. 20 

24.7 

3.9 

0. 32 

0. 58 

0. 84 (3) 

Yes? 

Yes? 

20 







S 

0. 15 

0. 22 

23.6 

3.8 

0. 59 

0. 28 

0. 40 

No 

No 

21 







S 

0. 14 

0.20 

31.7 

5.0 

0. 36 

0.28 

0.40 

No 

No 

22 







F 

0. 14 

0. 20 

31.7 

5.0 

0.40 

0.28 

0.40 

No 

No 

23 





1 

2 

F 

0. 14 

0.20 

31.8 

5. I 

0. 47 

0.26 

0. 38 

No 

No 

24 





9 

2 

F 

0. 14 

0.20 

27.2 

4. 3 

0. 46 

0.21 

0.31 

Yes? 

Yes? 

25 







F 

0. 12 

0. 18 

31.7 

5. 0 

0. 33 

0. 23 

0. 34 

Yes? 

No 

26 

200 X 

1400 

Pl€ 

a ted 

9 

2 

S 

0. 12 

0. 18 

31.8 

5. 1 

0. 32 

0.24 

0. 35 

No 

No 

27 

200 X 

1400 Screen Backup 

9 

2 

F 

0. 21 

0. 30 

12. 5 

2.0 

0.92 

0. 34 

0. 49 

Yes? 

Yes? 

28 







S 

0. 17 

0. 24 

12. 5 

2.0 

0.95 

0. 30 

0. 44 

Yes? 

Yes? 

29 







S 

0. 15 

0.22 

6. 3 

KO 

1.55 

0.27 

0. 39 

Yes? 

No 

30 





9 

2 

F 

0. 17 

0.24 

7.2 

1. 1 

0. 75 

0.25 

0. 36 

Yes? 

No 

31 





1 

2 

F 

0. 17 

0.25 

14.5 

2.3 

0.60 

0. 34 

0. 49 

Yes? 

Yes? 

32 







s 

0. 15 

0.22 

14.2 

2.2 

0. 69 

0. 33 

0.4 B 

Yes? 

Yes? 

33 







s 

0. 17 

0. 24 

11. 1 

l.B 

0.83 

0. 30 

0.44 

Yes? 

Yes? 

34 

200 X 

1400 Screen Backup 

I 

2 

F 

0. 19 

0. 27 

12.5 

2.0 

0. 33 

0. 25 

0. 36 

Yes? 

Yes 

35 

720 X 

140 

Plain 

I 

2 

F 

0.22 

0. 32 (5) 

20. 7 

3. 3 

0. 30 

0. 62 

0. 90(3) 

Yes? 

Y es 

36 







s 

0.22 

0. 32 (5) 

19.8 

3. 1 

0. 30 

0. 62 

0. 90(3) 

Yes? 

Yes 

37 







s 

0. 17 

0.25 

• 22.9 

3. 6 

0. 35 

0. 62 

0. 90(3) 

Yes? 

Yes 

3 B 





1 

2 

F 

0. 19 

0.27 

22.9 

3.6 

0. 35 

0. 62 

0. 90 (3) 

Yes? 

Yes 

39 





9 

2 

F 

0. 17 

0.25 

18.7 

3.0 

0. 90 

0. 19 

0. 26 

Yes? 

Yes 

40 







S 

0. 17 

0. 25 

19.0 

3. 1 

0.95 

0. 19 

0. 28 

Yes? 

Yes 

41 







S 

0. 18 

0.26 

21.9 

3.5 

0. 30 

0.25 

0. 36 

Yes? 

Yes 

42 

720 X 

140 

Plain 

9 

2 

F 

0. 15 

0.22 

21.9 

3. 5 

0. 29 

0. Z 6 

0. 38 

Yes? 

Yes 

43 

165 X 800 

Plain 

9 

2 

F 

0. 17 

0. 25 

14.2 

2. 2 

0. 90 

0. 21 

0. 30 

Yes 

No 

44 







S 

0. 17 

0. 24 

14.2 

2. 2 

1. 1 

0. 20 

0. 29 

Yes 

No 

45 







S 

0. 17 

0. 24 

17.2 

2. 7 

1.2 

0. 21 

0. 30 

Yes 

No 

46 





9 

2 

F 

0. 16 

0.23 

18.2 

2.9 

1.4 

0. 21 

0. 31 

Yes 

No 

47 





1 

2 

F 

0. 19 

0. 26 

16.2 

2.6 

0. 53 

0. 26 

0. 38 

Yes 

Yes 

48 







S 

0. 19 

0.28 

17.2 

2.7 

0. 53 

0.23 

0. 34 

Yes 

Yes 

49 







S 

0. 18 

0. 26 

20. 5 

3. 2 

0. 54 

0.23 

0. 34 

Yes 

Yes 

50 

165 X 

800 

Pla 

in 

I 

2 

F 

0. 18 

0.26 

20. 5 

3.2 

0. 60 

0.23 

0. 34 

Yes 

Yes 

51 

165 X 

800 

Pleated 

1 

2 

F 

0. 21 

0.31 (5) 

30. 3 

4.8 

0. 14 

0. 62 

0. 90 (3) 

Yes? 

Yes 

52 







5 

0.21 

0. 30 (5) 

29.5 

4.7 

0. 11 

0. 62 

0.90 


Yes? 

Yes 

53 







S 

0. 17 

0. 24 

32.2 

5. 1 

0. 11 

0. 62 

0.90 


Yes? 

Yes 

54 





1 

2 

F 

0. 19 

0. 27 

32.2 

S. 1 

0. 14 

0. 62 

0.90 


Yes? 

Yes 

55 





9 

2 

F 

0. 14 

0. 20 

30. 0 

4.8 

0.26 

0.62 

0.90 


Yes? 

Yes 

56 







S 

0. 12 

0. 17 

30.0 

4.8 

0.26 

0. 62 

0.90 


Yes? 

Yes 

57 







S 

0. 13 

0. 19 

33.4 

5. 3 

0. 30 

0. 62 

0.90 


Yes? 

Yes 

58 

165 X 

800 

Pleated 

9. 

2 

F 

0.08 

0. 11 

33.4 

5.3 

0.23 

0. 62 

0.90 ( 

3) 

Yea? 

Yes 

59 

165 X 

800 

Screen Backup 

9. 

2 

F 

0. 10 

0. 14 

14.8 

2. 3 

1.0 

0. 15 

0.22 

No 

No 

60 







S 

0. 11 

0. 16 

14.8 

2. 3 

1.0 

0. 15 

0. 22 

No 

No 

61 







S 

0. 12 

0. 17 

17.2 

2. 7 

1.0 

0. 17 

0. 24 

Yes 

No 

62 





9 

2 

F 

0. 11 

0. 16 

17.2 

2.7 

1.5 

0. 17 

0. 24 

Yes 

No 

63 





1 

2 

F 

0. 10 

0. 14 

17.2 

2. 7 

0.2 

0.14 

0.21 

No 

No 

64 







S 

0. 10 

0. 15 

17.2 

2. 7 

0.2 

0. 21 

0.31 

Yes 

No 

65 







S 

0. 10 

0. 15 

19.2 

3. 0 

0. 38 

0.21 

0.31 

Yes 

Yes? 

66 

165 X 800 

Screen 

Backup 

1 

2 

F 

0. 10 

0. 15 

19. 8 

3. 1 

0.2 

0. 22 

0.32 

Yes 

Yes? 

67 

500 X 

500 

Plain 

1 

2 

F 

0. 13 

0. 19 

34. 3 

5.4 

0. 13 

0. 62 

0.90 (3) 

Yes? 

Yes 

68 







S 

0. 18 

0. 26 

32.9 

5. 2 

0. 13 

0. 62 

0.90 


Yes? 

Yes 

69 







S 

0.20 

0. 29 (5) 

34. 3 

5.4 

0. 13 

0. 62 

0.90 


Yes? 

Yes 

70 





1 

2 

F 

0. 12 

0. 17 

34.3 

5.4 

0. 37 

0. 62 

0.90 


Yes? 

Yes 

71 





9 

2 

F 

0.21 

0. 30 (5) 

31.0 

4.9 

0. 36 

0. 62 

0.90 


Yes? 

Yes 

72 







S 

0.28 

0.41 (5) 

31.7 

5. 0 

0.45 

0. 62 

0.90 


Yes? 

Yes 

73 







s 

0. 14 

0. 21 

35. 0 

5.6 

0.24 

0. 62 

0.90 ( 


Yes? 

Yes 

74 

500 X 

500 

Pla 

in 

9 

2 

F 

0. 05 

0. 07 

35.0 

5.6 

0.26 

0. 19 

0.2? 

Yes? 

No? 

75 

200 X 

1400 

Perf. Sht. 

9 

2 

F 

0. 08 

0. 11 

10.2 

1.6 

1. 1 

0. 19 

0.27 

No 

No 

76 



No. 

1 



S 

0. 09 

0. 13 

9.3 

1. 5 

1. 1 

0. 19 

0.28 

No 

No 

77 







s 

0. 14 

0. 20 

7.2 

i . 1 

1.0 

0.21 

0.31 

No 

No 

78 







F 

-0. 14 

0. 20 

8.4 

1. 3 

1.2 

0.23 

0.34 

No 

No 

79 





9 

2 

F 

0. 14 

0.21 

10.2 

1.6 

1.2 

0. 30 

0.43 

No 

No 

80 





1 

2 

F 

0. 15 

0. 22 

7.8 

1.2 

1.0 

0. 24 

0.35 

No 

No 

81 







S 

0. 15 

0. 22 

7.8 

1.2 

1.2 

0. 24 

0. 35 

No 

No 

82 







s 

0. 16 

0. 23 

7.2 

1. 1 

1.0 

0.23 

0. 34 

No 

No 

83 

200 X 

1400 

Perf. Sht. 

1. 

2 

F 

0. 15 

0. 22 

8.9 

1.4 

1.0 

0.27 

0. 39 

No 

No 




No. 

‘ 








1 






77 























Table 13 


INITIAL FLOW PULSE TEST DATA (Continued) 


Teat 

No. 

Screen 

Configuration 

Line 

Length 

(m) 

Valve 

Opening 

Impoiied 
Screen>Ullage 
Preaaare Difference 

N/cm^ psl 

Steady Flowrate 
(Approx) 

/Sec X 10^ gpm 

Time 

After 

Start 

(Sec) 

Flew Pulae 
Screen* Ullage 
Preaeure DLtference 

N/cm^ pai 

Gaa Ingeatlon 
Observed 

Babble Oscillograph 
Trap (Flowmeter) 

84 

200 X 

1400 

Perf, Sht. 

1 

2 

F 

0.15 

0. 22 

5.1 

0.8 

O.S 

O.IB 

0. 26 


Yes? 

No 


85 



No. 

7 



S 

0. 18 

0. 26 

S.9 

0.9 

0. S 

0.24 

0.35 


Yes? 

No 


86 







s 

0. 18 

0. 26 

9.3 

1.5 

0.4 

0. 35 

0. 51 


Yes? 

Yes 


87 

200 X 

1400 

Perf. Sht. 

1 

2 

F 

0. 18 

0.26 

6.6 

1.0 

0.4 

0. 28 

0.41 


Yet? 

Yea? 





No. 

2 















88 

165 X 

800 

Perf Sht. 

1 

2 

F 

0.62 

0. 90 (5) 

8.9 

1.4 

0. 20 

0. 62 

0.90 (3) 

Yes? 

Yes? 


89 



No. 

1 



S 

0.62 

0.90 1 

9.3 

1.5 

0. 22 

0. 62 

0.90 


Yes? 

Yes 


90 







S 

0.62 

0.90 \ 

8.4 

1.3 

0. 18 

0.62 

0.90 


Yes? 

Yes 


91 

165 X 

BOO 

Perf. Sht. 

1 

2 

F 

0.62 

0. 90 < 5) 

7.2 

1. 1 

0. 18 

0. 62 

0.90 (3) 

Yes? 

Yes 





No. 

1 















laoprnpyl Alcohol *>' 

Vertical Specimen 














92 

200 X 

1400 

■ Plain 

9 

2 

F 

0. 18 

0. 26 

8.9 

1.4 

1. 75 

0. 27 

0. 39 


No (4) 

No 


93 







S 

0. 26 

0.37 

8.4 

1.3 

1. 0 

0. 30 

0.43 


No 1 

No 


94 







S 

0.21 

0. 30 

11. 1 

l.B 

1.6 

0, 37 

0.53 


No 1 

No 


95 

200 X 

1400 

Plain 

9 

2 

F 

0. 19 

0. 28 

U.2 

1.9 

1.9 

0. 37 

0.53 


No (.) 

No 


96 

325 X 

2300 

Plain 

9 

2 

F 

0. 32 

0. 47 

9. 3 

1. 5 

1.9 

0. 47 

0. 68 


No (4) 

No 


97 







S 

0. 17 

0.25 

9. 3 

1.5 

2.0 

0.42 

0.61 


No 1 

No 


98 







S 

0. 24 

0.35 

14.2 

2.2 

1.6 

0. 50 

0.73 


No 1 

No 


99 

325 X 

2300 

Plain 

9 

2 

F 

0. 36 

0. 52 

10. 7 

1.7 

1.6 

0.51 

0.74 


No (4) 

No 


100 

325 X 

2300 

Plain (5) 

9 

2 

F 

0.08 

0. 12 

(21.9)15) (3.5) 

1. 34 

0. 34 

0. 49 


No (4) 

No (9) 


101 







S 

0.07 

0. 10 

(2). 9) 

(3.5) 

1.70 

0. 30 

0.43 


No } 

No 1 


102 







S 

o.oe 

0. 11 

(29.1) 

(i.3) 

1.45 

0.47 

0.68 


Yes 

No 


103 







F 

0.08 

0. 12 

(39.1) 

(6.2) 

1. 50 

0.45 

0.65 


Yes i 

No 


104 

325 X 

2300 

Plain <5) 

9 

2 

F 

0. 10 

0. 15 

(39.1)(5) (6.2) 

1.50 

0. 46 

0. 66 


Y«« (4) No (9) 


F reon 

114 (Saturated at 0.0 N/ctn^] 

— Vertical Specir 

nen 











105 (6) 

325 X 

2300 

Plain 

9 

2 

F 

0.28 

0.41 

17.7 

2.8 

1.0 

0. 35 

0.51 


No (4) No 


106 (6) 







S 

0.28 

0. 40 

17. 7 

2.8 

1. 1 

0.34 

0. 50 


No 

No 


107 (6) 







S 

0.24 

0.35 

19.2 

3.0 

1. 0 

0.27 

0. 39 


No 

No 


108 (6) 







F 

0.27 

0. 39 

11.8 

1.9 

1. 0 

0. 29 

0. 42 


No? 

No? 


109 







F 

0.43 

0.63 (8) 

15. 3 

2.4 

1.0 

0. 48 

0.70 


Yes? 

Yes? 


no 







S 

0. 31 

0. 45 

11. 1 

1.8 

0.4 

0. 62 

0.90 (3) 

No 

No? 


in 







s 

0.41 

0. 59 

18.2 

2.9 

0. 8 

0. 46 

0. 66 


No 

No? 


112 

325 X 

2300 

Plain 

9 

2 

F 

0.40 

0. 58 

18.7 

3.0 

1. 3 

0.45 

0.65 


No (4) No? 


113 (7) 

200 X 

1400 

Plain 

9 

2 

F 

0.46 

0.67 (8) 

20. $ 

3.2 

0. 26 

0.55 

0. 80 


Yes?(4) Yes? 


114 (7) 







S 

0.44 

0.64 (8) 

13.2 

2. 1 

0. 29 

0. 59 

0.85 


Yes? 

Yes? 


115 (7) 







s 

0. 33 

0. 48 (8) 

18.7 

3.0 

0.44 

0. 62 

0.90 


No 

1 No? 


116 (7) 

200 X 

1400 

Plain 

9 

2 

F 

0. 27 

0. 39 

21.9 

3.5 

2.0 

0. 33 

0.48 


No ( 

4) No? 


T reon 

1 14.Horisontal Specimen 















117 

325 X 

2300 

Plain 

9 

2 

F 

0.62 

0.90 (8) 

18. 7 

3.0 


0.62 

0.90 (3) 

Yes 

Yes (10) 

118 







S 

0.44 

0. 64 (8) 

20. 5 

3.2 

0.4 

0. 62 

0.90 


Yes 

Yes 


119 







S 

0. 52 

0.76 (8) 

8.9 

1.4 

0. 4 

0. 62 

0.90 


Yes 

No? 


120 





9 

2 

F 

0.43 

0.63 

17.2 

2.7 

0. 4 

0.62 

0.90 


Yes 

Yes 


121 





1 

2 

F 

0.50 

0. 72 

21.9 

3.5 

0. 15 

0.62 

0.90 


No 

Yes 


122 







S 

0. 50 

0. 72 

24.9 

4.0 

0. 15 

0. 62 

0.90 


No 

Yes 


123 







s 

0. 51 

0.74 

20. 5 

3.2 

0. 15 

0. 62 

0.90 


No 

Yes 


124 

325 X 

2300 

Plain 

1 

2 

F 

0. 51 

0. 74 

20.7 

3.3 

0. 15 

0. 62 

0.90(3) 

No 

Yes (10) 

125 

200 X 

1400 

Plain 

1 

2 

F 

0. 42 

0.61 

20. 5 

3.2 

0. 30 

0. 62 

0.90(3) 

No 

Yes (10) 

126 







S 

0. 43 

0. 62 

22.9 

3.6 

0. 30 

0.62 

0.90 


No 

Yes 


127 







S 

0. 34 

0. 50 

14.8 

2.3 

0. 30 

0.62 

0.90 


No 

Yes 


128 

200 y 

1400 

Plain 

1 

2 

F 

0. 43 

0. 62 

17. 7 

2.8 

0. 36 

0. 62 

0.90( 

) 

No 

Yes (10) 

129 

200 X 

1400 

Pleated 

1 

2 

F 

0. 34 

0.50 

25.4 

4.0 

0. 28 

0,44 

0. 64 


Yes 

Yet (10) 

130 







S 

0. 39 

0.57 

25.4 

4.0 

0. 34 

0. 48 

0.70 


Yes 

Yes 


131 







S 

0. 46 

0. 70 

18.2 

2.9 

0. 36 

0. 55 

o.eo 


Yes 

Yes 


132 

200 X 

1400 

Pleated 

1 

2 

F 

0. 54 

0.78 

17. 7 

2.8 

0. 40 

0. 59 

0.86 


Yes 

Yes (10) 

133 

200 X 

1400 

Screer 

Backup 

1 

2 

F 

0.41 

0. 59(8) 

20. 5 

3.2 

0. 42 

0. 60 

0.87 


No 

Yes (10) 

134 







S 

0. 26 

0. 38 

25.4 

4.0 

0. 35 

0. 59 

0.86 


No 

Yes 


135 







S 

0. 16 

0.23 

19.8 

3. 1 

0.41 

0.31 

0. 45 


No 

Yes 


136 

200 X 

1400 

Screen Backup 

1 

2 

F 

0. 19 

0.27 

21.7 

3.4 

0.41 

0. 34 

0. 50 


No 

Yes (10) 

137 

200 X 

1400 

Perf. Sht. 


2 

F 

0. 26 

0.40 (8) 

13.2 

2. 1 

0.44 

0. 60 

0.87 


Yes? 

Yo. (10) 

138 



No. 

1 



S 

0.46 . 

0.66 

11.8 

1.9 

0. 31 

0. 62 

0.90 (3) 

Yet? 

Yes (10) 

139 







S 

0. 34 

0. 50 

17.2 

2.7 

0. 60 

0. 47 

0. 68 


No 

No? 


140 

200 X 

1400 

Perf. Sht. 

1 

2 

F 

0.46 

0. 66 

13.2 

2. 1 

0.45 

0.51 

0.74 


No 

No? 





No. 

1 















141 

720 X 

140 

Plain 

1 

2 

F 

0.53 

0. 77 

22.9 

3.6 

0. 26 

0. 62 

0.90 (3) 

Yes? 

Yes (10) 

142 







S 

0. 55 

0. 79 (8) 

21.9 

3.5 

0. 30 

0. 62 

0.90 


No 

Yes 


143 







s 

0. 58 (5) 

0.84 

26.4 

4.2 

0.25 

0.62 

0.90 


No 

Yes 


144 

720 X 

140 

Plain 

1 

2 

F 

0. 62 (5) 

0. 90 

25.4 

4.0 

0.23 

0.62 

0.90 (3) 

No 

Yes (10) 

145 

165 X 

800 

Plain 

1 

2 

F 

0. 30 (5) 

0. 44 

17,2 

2.7 

0. 26 

0.41 

0.60 


Yes 

Yob (10) 

146 







s 

0.52 1 

0. 76 

13.2 

2. 1 

0. 38 

0. 62 

0.90(3) 

No 

Yes 


147 







S 

0. 57 1 

0.83 

18.7 

3.0 

0. 27 

0.62 

0.90 (3) 

Yes 

Yes 


148 

165 X 

800 

Plain 

1 

2 

F 

0.62 (5) 

0,90 

14.8 

2.3 

0.25 

0.62 

0.90 (3) 

Yes 

Yes (10) 

149 

500 X 

500 

Plain 

1 

2 

F 

0. 39 (5) 

0.57 

23.5 

3.7 

0. 29 

0. 58 

0. 84 


Yea 

Yes (10) 

150 







S 

0.59 1 

0.65 

20.7 

3.3 

0. 25 

0. 62 

0.90 (3) 

Yes 

Yes 


151 







s 

0.45 1 

0.65 

29.1 

4.6 

0. 23 

0.62 

0.90 (3) 

Yes 

Yes 


152 

500 X 

500 

Plain 

1 

2 

F 

0. 46 (5) 

0. 66 

27.2 

4.3 

0. 24 

0.62 

0.90(3) 

Yes 

Yes (10) 
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Table 13 


INITIAL FLOW PULSE TEST DATA (Continued) 


Test 

No. 

Screen 

Configuration 

Line 

Length 

(m) 

Valve 

Opening 

U) 

Imposed 
Screen-Ullage 
Pressure Difference 

N/cm^ psi 

Steady Flowrate 
(Approx) 

m^/Sec X 10^ gpm 

— 

Time 

After 

Start 

(Sec) 

Flow Pulse 
Screen-Ullage 
Pressure Difference 

N/cm^ psi 

Gas Ingestion 
Observed 

Bubble Oscillograph 
Trap (Flowmeter) 

l-»2 

- Horizontal 

Speciman 














(9) 


153 

325 X 

2300 

Plain 

1 

2 . 

F 

0.03 

0.04 

1.4 

2.0 

0. 16 

0.32 

0.46 


Yes? 

Yes? 

(10) 

154 









0. 03 

0.045 

1.4 

2.0 

0. 11 

0. 33 

0.48 


No? 

Yes? 

1 

155 









0. 26 (5) 

0.41 

0.7 

1.0 

0. 12 

0. 50 

0.72 


No? 

Yes? 

1 

156 









0. 017 

0. 025 

0. 7 

i.o 

0. 11 

0. 31 

0.45 


No? 

No? 

1 

157 









0. 02 

0. 03 

0. 0 

0.0 

0. 16 

0.46 

0. 66 


Yes? 

Yes 

(10) 

158 

325 X 

2300 

Plain 

1 

2 

F 

0.23 (5) 

0. 33 

0. 0 

0.0 

0. 15 

0. 52 

0. 75 


Yes? 

No 


159 

200 X 

1400 

Pleated 

1 

2 

r 

0. 02 

0.03 

1.4 

2.0 

0. 14 

0. 19 

0. 28 


_ 

_ 

(10) 

160 









0.017 

0. 025 

1.4 

2.0 

0. 14 

0. 19 

0. 28 


- 

- 

1 

161 









0.05 (5) 

0. 07 

0. 0 

0.0 

. 



(11) 


- ■ 

1 

162 

200 X 

1400 

. Pleated 

1 

2 • 

F 

0.23 (5) 

0. 34 

0. 0 

0.0 

- 

- 

- 

(11) 

- 

- 

(10) 

163 

200 X 

1400 

Screen Backup 

1 

2 

F 

0. 03 

0. 045 

1.4 

2.0 

0. 10 

0. 15 

0. 22 


No 

No 


164 









0. 03 

0.05 

1. 4 

2.0 

0. 09 

0. 17 

0.24 


No 

No 


165 









0. 01 

0 . 02 

0.0 

0.0 

0. 12 

0. 06 

0. 09 


No 

No 


166 

200 X 

1400 

Screen Backup 

1 

2 

I 


0. 03 

0.045 

0. 0 

0.0 


- 

- 

(10) 

No 

No 


167 

200 X 

1400 

Plain 

1 

2 

F 

0.01 

0.02 

1. 4 

2. 0 




(10) 

No 

No 


168 









0.05 (5) 

0. 07 

1. 4 

2.0 

0. 11 

0.01 

0. 02 


No 

No 


169 









0. 60 (5) 

0. 87 

0. 0 

0.0 

_ 

0. 60 

0. 87 


? 

9 


170 

0 

0 

X 

1400 

Plain 

1 

2 

F 

0.60 (5) 

0.87 

0.0 

0. 0 


0.60 

0. 87 


? 

? 


LH 2 

- Vertical Specimen 
















171 

200 X 

1400 

Screen Backup 

1 

2 

t 


0.007 

0.01 

1.4 

2.0 

0. 11 

0. 06 

0. 09 


No 

No 


172 

200 X 

1400 

Screen Backup 

1 

2 

I 


0.007 

0. 01 

1.4 

2. 0 

. 

. 

- 

(10) 

No 

■- 


173 

200 X 

1400 

Screen Backup 

1. 

2 

i 


- 

- (JO) 

0. 0 

0.0 

- 

• 

- 

(10) 

No 

- 


174 

325 X 

2300 

Plain 



i 


0.0 

0.0 

1.4 

2.0 

0. 14 

0. 06 

0. 08 


No 

No 


175 









0.0 

0.0 

1.4 

2.0 

- 

. 

. 

(10) 

No 

No? 


176 










- (10) 

0.0 

0.0 



. 

(10) 

No 

No? 


177 

325 X 

2300 

Plain 


2 

F 

“ 

- (10) 

0. 0 

0.0 

• 

" 

* 


No 

No? 



(1) Fa Op«n 0. 01 S Sec; Close 0. 015 See 
S a Open 0. 019 Sec; Close 0. 019 Sec 

(2) ? a Possible Leakage Ingestion 

(3) Severe Oseillation«Transducer Saturated 

(4) Liquid Spill ^ Bubble Trap not used 

(5) Flowmeter Removed (Est, Flow) 

(6) Saturated at 2. 1 N/cm^ (3. 0 Psig) 

(7) Saturated at 0. 7 N/cm^ (1, 0 Psig) 

(8) Excessive AP Inadvertently Imposed 

(9) Oscillograph — AP Transducer 

(10) Line Boiling 

(11) No ^P Response 
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Figure 34. Oscillograph Record Showing Startup Gas Ingestion 
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Figure 36. Run 22 Startup — 0.01 5 Second Valve — Short Line 
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Figure 37. Run 25 Startup - 0.015 Second Valve — Long Line 
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times, but differed from the short-line startup trace (compare Figure 36) 
in taking 0. 1 sec longer to reach full flow. Because the flow acceleration is 
slower, the screen-ullage pressure difference responds to a lower value for 
the long-line tests, compared to the short-line tests. 

It was also anticipated that the combination of 0. 015-sec valve closing time 

combined with the long line would result in waterhammer during shutdown. 

Waterhammer only occurred in seven tests and always at high flow rates. 

The most severe occurred with the pleated 200 x 1400 mesh (Runs 24 to 26) 

and happened with both valve closure times. The shutdown trace for Run 25 

(0. 015-sec valve closure) is shown in Figure 38, and for Run 26 (0. 019-sec 

valve closure) in Figure 39. Although the initial pulse in both tests saturated 

2 

the transducer, the magnitude of the pulse was estimated at 3. 5 N/cm 
(5 psi). The response of the screen-ullage pressure differential and flow 
was similar for both tests. It was not clear why the waterhammer frequency 
was lower for Run 25 with the 0,015-sec valve (Figure 38); the presence of an 
ingested gas bubble in the bubble trap for Run 25 may have contributed to the 
difference. Much less severe waterhammer occurred with the pleated 165 x 
800 mesh (Runs 55 and 58). Although the flow rates for Runs 55 and 58 were 
just as high as for Runs 24 to 26, it is thought that the gas ingestion which 
occurred in the former runs reduced the waterhammer intensity by reducing 
the effective sonic velocity in the pipe. This thesis was further confirmed 
by the 500 x 500 mesh tests (Runs 71 to 74) where, despite high flow rates, 
the severe gas ingestion apparently eliminated waterhammer. The water- 
hammer phenomenon is further discussed in the following section on data 
correlation. 

Except for the minimal screen response to waterhammer with the long line', 
neither the valve open or close time nor line length had any significant effect 
on screen device response beyond that expected and capable of being verified 
by the simulations using the H672 code with the screen model. 

Examination of the initial startup flow pulse data shown in Table 13 revealed 
that many of the tests were subjected to excessive gas pressure, which led 
to gas ingestion regardless of conditions of flow, valve open/close speed, 
or line length. 
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Although it was anticipated that the minimum hydrostatic head imposed on 
the specimen would be about 5 cm (equivalent to 0. 04 N/cm , or 0. 06 psi), 
as discussed in the previous section, it was found. that the minimum head 
which could be reliably imposed (so that bubbles were slowly and continuously 
emitted) was about 0. 17 N/cm^ (0.25 psi). , (See Table 13. ) It is thought that 
this may have been due to pressure drop through the very small remote- 
actuated bleed valve. This imposed head was quite close to the bubble point 
of the coarser mesh screens (see Table 12), which generally led to gas 
ingestion at flow startup. In addition, for some specimens there was minor 
leakage from the plexiglass gas enclosure which led to difficulties in clearing 
liquid from around the screen within the gas enclosure. When additional gas 
(pressure) was introduced to overcome this leakage and surround the screen 
with gas, the imposed pressure was inadvertently excessive (see, for exam- 
ple, Runs 88-91 of Table 13). 

For many runs, gas ingestion occurred at startup as evidenced by gas 
bubbles caught in the bubble trap. Usually the quantity of gas ingestion was 
easily visible and readily measured. As noted from Table 13, apparently 
with some runs (e. g. , compare Runs 24 and 25) not all gas ingestion was 
caught in the bubble trap and was detected by the flowmeter, as discussed 
previously. 

Although startup gas ingestion was observed for many alcohol tests, liquid 
spill was never observed for horizontal specimens even at shutdown where 
waterhammer occurred. On the other hand, for the vertical specimens, 
liquid spill occurred on occasion during startup. Since the bubble trap could 
not be used with the vertical orientation, liquid spill was construed to indi- 
cate gas ingestion for these tests. 

Despite the fact that the 104 alcohol tests were performed with all 13 test 
specimens listed in Table 8, the problem of inadvertent gas overpressure 
which occurred with many of the lower bubblepoint screens (e. g. , Runs 27-58, 
67-73, 88-91) effectively masked the effects of other parameters and reduced 
the usefulness of these tests in providing data for analysis and correlation. 
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F reon 1 14 Tests 

The experiment matrix for the Freon 114 tests included 48 tests using 10 of 
the screen/structure specimens (the 165 x 800 mesh structural backup spe- 
cimens were not tested) as shown in Table 13. The number of F reon 114 
tests performed was somewhat reduced from that planned in Table 9i because 
it was found after about the first 16 tests (see Table 13) that it was imprac- 
tical to continue testing the "long" line configuration. The reason for this 
was that the "long" line could not be kept wet up to the valve because heat 
transfer and line boiling resulted in reverse flow and emptying of the line 
before the run could be started. This occurred even with preflow to fill and 
chill the line before actual test initiation. For the same reason, it was found 
impossible to obtain reliable bubble point data for the specimens using 
Freon 114; therefore, the alcohol bubble point data were extrapolated to 
Freon 114 properties, as shown in Table 12. Another consequence of the 
line boiling was gas backflow into the bubble trap, which made if difficult to 
empty the bubble trap of gas prior to the run. Thus for many runs the bubble 
trap was nearly empty and gas ingestion could not be measured in the bubble 
trap but was inferred from liquid spill from the specimen (e. g. , Runs 119 to 
152). The partly empty bubble trap may also have contributed to gas inges- 
tion observed by the oscillograph for Runs 119 to 152. Therefore, the most 
reliable runs for analysis of hydrodynamic effects were the runs with vertical 
orientation when the bubble trap was not used. 

The first tests (Runs 105 to 108) were with the Freon 114 saturated at 2. 1 
N/cm^ (3. 0 psig), which is equivalent to a fluid temperature of 282. 2 K 
(48 °F). This was quite close to the line temperature, hence significant line 
boiling did not occur. Note that as the line got more chilled down, from 
Run 105 to 107, the screen-ullage pressure difference got smaller (the low 
flowrate of Run 108 was due to low fluid level relative to the screen at 
startup). For tests run with the Freon 114 saturated at 0. 0 N/cm (0. 0 psig), 
the reverse flow caused by line boiling tended to fill the gas enclosure with 
liquid, which resulted in inadvertent gas overpressure (e. g. , Run 109) that 
caused gas ingestion. Even when the line was chilled enough so that backflow 
did not cause gas overpressure, the flow into the partially empty line caused 
a large spike in screen-ullage pressure differential (Run 110) but inconclusive 
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encountered during the Freon 114 tests of line boiling, inadvertent over- 
pressure in the gas enclosure, backflow, surging, and liquid spill also 
occurred to some degree during the tests. Additional problems pecu- 

liar to testing with also occurred. These included a reduced amount of 
instrumentation which included only platinum resistance temperature trans- 
ducers, screen-ullage pressure difference, pressure upstream of the cryo- 
genic ball valve, and a dewar pressure gage. In addition, on occasion the 
screen-ullage pressure differential transducer was either noisy or failed to 
function (see Note 11 in Table 13). The dewar vacuum was lost initially and 
necessitated repair of a dewar window, which resulted in good dewar per- 
formance for the balance of the test program. Although great care was taken 
to slowly chill and fill the dewar, apparent stress concentrations in the plexi- 
glass gas enclosure caused them to crack and leak gas in some specimens. 

The gas enclosure around the 200 x 1400 pleated specimen cracked so severely 
that gas was not in contact with the specimen (see Runs 159 to 162 in Table 13), 
and for the 200 x 1400 plain specimen only partial gas coverage of the specimen 
was obtained (Runs 167 to 170, Table 13). 

It was anticipated that the effect of saturation pressure on line boiling cryo- 
surge could be determined; these effects were apparently masked by the 
degree of line chilldown which contributed more significantly to the cryo surge 
effect. The line downstream of the valves was insulated and became progres- 
sively more chilled with each test, and it was not practical to allow the line 
to warm up between tests. Since the line temperature distribution after each 
test was not well defined, the first test in each series was the optimum test 
for analysis in terms of a known temperature distribution in the downstream 
line . 

2 2 

The dewar driving pressures of 3. 5 N/cm (5 psig) and 6. 9 N/cm 
(10 psig) were selected to give comparable volumetric outflow rates as for 
the alcohol and Freon 114 tests. Because the flowmeter resistance was not 
present, and the downstream resistance was reduced (the line downstream of 
the valve was empty), and the density and viscosity are low, the initial 

LH^ flowrate was higher than anticipated. This led to dynamic pressures 
which were excessive compared to the low bubble point of screens in 
(see Table 12) and generally resulted in startup gas ingestion. 
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gas ingestion. That is, the flowmeter trace on the oscillograph showed some 
"noise, " but not the severe spiking shown in Figure 34. The flowrate for 
Run 110 was probably low due to vapor choking (boiling) in the line; the next 
two runs showed reasonable flow rates, but because the imposed screen- 
ullage pressure differential was close to the screen bubble point, inconclu- 
sive gas ingestion occurred. Obvious specimen movement also occurred as 
a result of the line boiling and surging. 

The differences in screen device performance between tests with Freon 114 
saturated at 0. 0 N/cm or saturated at higher pressures was basically that 
the higher temperature (saturation pressure) fluid produced more reliable 
hydrodynamic effects because the effects of line boiling were minimized. 

On the other hand, lower temperature fluid led to line boiling, partially empty 
lines, and possible large flow surges and pressure spikes as described above. 

Because of the erratic flowrate and pressure behavior caused by line boiling 
and emptying, it was decided to restrict the remaining Freon 114 tests (Runs 
121-152) to the "short" line configuration in the hope that the line could be 
kept wet up to the valve, and the results of downstream line boiling on screen 
response observed under repeatable flow conditions, Unfortimately, the same 
problems occurring during the alcohol tests of imposing excessive screen- 
ullage pressure differential (above the 5. 0-cm head planned) were magnified 
with the Freon 114 tests. Because of the increased density of Freon 114, 

2 

the gas head which could be reliably imposed was of the order of 0. 28 N/cm 

2 

(0.4 psi, compared with 0. 17 N/cm , or 0.25 psi, for alcohol; see above) 
which is near the bubble point of even the 2 00 x 1400 mesh because of the 
reduced surface tension of Freon 114 compared to alcohol (see Table 12). 
Efforts to reduce this imposed head were not successful because of the 
arrangement of the gas -bleed valve, and hence general gas ingestion was 
observed for most of the tests with the coarser screens. This problem had 
been anticipated with the tests which was why only the 325 x 2300 mesh 
and 200 x 1400 mesh screens were tested with LH^. 

Tests 

The experiment matrix for the tests included 25 tests using six of the 

screen/ structure specimens, as shown in Table 13. The problems 
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As was the case in the Freon 114 tests, the boiling at the valve resulted in 
backflow which tended to fill the bubble trap with gas, and which resulted in 
liquid spill as evidence of gas ingestion. Therefore, as was the case with 
the Freon 114 tests, the vertical specimen configuration would have produced 
more reliable data; unfortunately, the pressure difference transducer behaved 
erratically during most of these tests, as shown in Table 13. It was concluded 
that the first series of tests with the 325 x 2300 mesh specimen yielded the 
best test data and demonstrated the cleanest example of line boiling cryo- 
surge and screen response to the cryosurge. 

It was anticipated that modelling of these tests with both the P4557 cryosurge 
code and then with the H672 code would provide adequate verification of the 
screen model as well as the capability of the analyses to model test specimen 
response with LH^. The modelling of the LH^ tests is described in detail in 
the next section. 






DATA CORRELATION 


In order to obtain correlation of the experimental data with the transient 
effects analysis, it was necessary to model the transient test apparatus and 
screen specimen for simulation with the H672 computer code. The objective 
of the data correlation was to verify important aspects of the analytical repre- 
sentation of the screen device. These aspects include (1) modelling of the 
screen device as a uniform composite structure with orthotropic structural 
properties based on a combination of screen and backup material properties, 
and (2) modelling of the gas ingestion, liquid inflow /outflow, and liquid spill 
characteristics of the screen device. The first aspect would be effected by 
comparing the performance of plain screen, pleated screen, screen backed 
up with coarse screen, and screen backed up with perforated plate, all of 
which have different screen fractions, f, and for some, different structural 
compliance. The second aspect leads from the first and would be verified by 
prediction of gas ingestion (or not) and ingested quantity under conditions 
which would occur experimentally. These predictions should be effected for 
conditions of startup, shutdown, and for all three experimental fluids. 

TEST APPARATUS MODEL 

In order to obtain realistic simulation, the test apparatus must be modelled 
accurately in order to achieve the necessary flowrates at the prescribed con- 
ditions of pressure and configuration. The test apparatus model is shown 
schematically in Figure 40, which identifies the heads, Hj, etc., and flows, 
Qj, etc. , throughout the test apparatus. The flow system starts with a tank 
{ 1), modelled by subroutine SR05, connected to a resistance (2), modelled 
by SR02. In actuality, the screen specimen is dead-ended and all outflow is 
through the screen, but in order to be compatible with the H672 program 
requirements, a very high resistance was used so that the flow through it is 
less than 0. 1% of the screen flow. Following the resistance comes (3) the 
"wet" screen, ( surrounded by liquid), and (4) the "dry" screen (surrounded 
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by gas), both modelled by SR26. Next comes a branch (5), SR03, which leads 

to the bubble trap (6), modelled as a dead-ended pipe, SR 14. The branch also 

leads to the flowmeter, (7), modelled as a resistance SR02, then a simple 

pipe junction (8), modelled by SROl. Next is another resistance (9), SR02, 

which accounts for the pipe frictional resistance between the flowmeter and 

control valve (10) modelled by SR 11. Following the control valve is the 

throttle valve (11), also modelled as a resistance, SR02, which includes also 

the frictional resistance of the pipe downstream of the throttle valve. Finally, 

another simple junction (12), SROl, connected to the outlet in the receiver 

tank, (13), modelled by an injector, SR12. The configurational details of 

each component are shown in Table 14. The actual pipe lengths and sizes 

were used, except for the bubble trap and screen devices, in which artificial 

wall thicknesses were used to make these lower modulus pipes modelled 

7 2 

accurately by steel pipes with a common modulus of 2. 07 x 10 N/cm (3 x 

7 

10 psi), as required by the H672 code. The actual tank size and average 
ullage height were used, together with the measured initial ullage pressure. 
The resistance values for the flowmeter and pipes were essentially based on 
frictional fluid flow analysis and modified as necessary to obtain the proper 
baseline flowrate at a baseline tank ullage pressure. Common values of 
resistance were used for all simulations either with the "long" or "short" 
line configuration. The valve open/close characteristics were modelled as 
accurately as possible from the experimental data; a total valve opening time 
of 0. 2 second was used (the effective valve opening time depended on the rela- 
tive resistance values throughout the flow circuit). Details of the screen 
devices depend on the specimen configuration and are discussed below. 

TESTS SELECTED FOR CORRELATION 

To verify the screen analytical model, only certain tests were selected for 
detailed analysis using the H672 code. This was necessitated by the fact that 
H672 is quite complex to set up, and fairly costly to run; further, many of the 
tests were inappropriate for simulation due to experimental problems of 
leakage, overpressure, etc., as discussed previously. 

In examining data from Table 13 to determine appropriate runs for simulation, 
it was noticed that, for many initial runs following specimen or line length 
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Table 14 

TEST APPARATUS CONFIGURATION MODEL FOR H672 CODE 


Components 

ID 

Pipe Parameters 
Length (ft) Type 

Datum 
Height (ft) 

1. 

SR05 

Tank 

T1 

0. 01 

1 

0. 0 

2, 

SR 02 

Resistance 

T2 

0. 13 

1 

0. 0 

3. 

SR26 

Sc reen (wet) 

T3 

0. 30- 

1 

0. 0 

4. 

SR26 

Screen (dry) 

T4 

0. 22 

1 

0. 0 

5, 

SR 03 

Branch 

T5 

0. 33 , 

3 

0. 0 

6. 

SR 14 

Bubble Trap 

T6 

0. 83 

2 

-0. 33 

7. 

SR 02 

Flowmeter 

T7 

1. 08 

4 

-0. 67 

8. 

SROl 

Junction 

T8 

15, 7 0 L- 0. 70S'!' 

4 

-1. 83 

9. 

SR02 

Resistance 

T9 

15. 70 L 0. 70S 

4 

-3, 50 

10. 

SR 11 

Valve 

TIO 

0. 50 

4 

-3. 50 

11. 

SR02 

Throttle Valve 

Til 

2. 50 L 17. SOS 

4 

-3. 50 

12. 

SROl 

Junction 

T12 

2, 50 L 17. SOS 

4 

-1. 75 

13. 

SR 12 

Injector 

T13 

0. 0 

4 

0. 0 

(14. 

SR18 

C ryosurge) 

= 

Long S = Short 




Pipe Types 


Diameter (in) 


Wall Thickness (in) 


1 

2 

3 

4 


Screen - depends on screen mesh and construction 


1.250 
1. 007 
1 . 000 


0. 033 
0. 0033 
0. 035 


Tank: 


Dia (in) Ullage Height (ft) 

16. 9 3. 0 


Volume (ft ) 
5. 58 


Valve: 1. 0 in Dia. Ball Valve, L/D = 3, C = 11.5, Open/Close Time 

= 0. 20 sec. ^ 


Injector: Open pipe; C^ = 4. 15; Back Pressure: 15.48 PSIA 
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changes, flowmeter "spiking" occurred during startup (e. g. , Runs 1, 5, 13, 

19, 24, etc), which may have been due to a partially empty line up to the 
valve. Such runs were basically not suitable for simulation since the H672 
code requires liquid-filled lines up to the valve. 

The first unambiguous and clearly observed and measured gas ingestion 
occurred with Runs 15 to 18. Run 15 was selected for analysis because it 
was at the lowest flowrate of the four runs (higher flowrate would increase 
the tendency for gas ingestion) and was also plain 200 x 1400 mesh. On the 
other hand. Runs 21 to 23 with pleated 200 x 1400 mesh, with 2. 5 times the 
flowrate as Run 15 (6. 25 times the velocity head) showed no tendency for gas 
ingestion; therefore, these runs (all at essentially identical head/flowrate 
conditions) were also selected for simulation. 

The next several series of runs were not suitable because of either leakage 
or excessive initial imposed screen-ullage pressure difference. The 165 x 
800 mesh with screen backup tests (Runs 60 and 6l) were selected because 
Run 61 has a slightly higher flowrate than Run 60; thus, there was clearly 
observed and measured gas ingestion for Run 6l, and Run 60 shows no evi- 
dence of gas ingestion. 

Unfortunately, the 500 x 500 mesh screen tests generally showed gas ingestion 
due to a combination of unintentional gas overpressure plus possible external 
screen leakage. One test (Run 74) had the least ingestion (one bubble of 
2. 5 cm, or 1. 0 in diameter, into the bubble trap, plus a minimal spiking 
region in the oscillograph trace for the flowmeter), so it was also selected 
for analysis. The previous test (Run 73) had more serious ingestion (3. 8 cm, 
or 1. 5 in, in bubble trap) plus general oscillograph spiking, but it too was 
selected for analysis. 

All of the tests of 200 x 1400 mesh (with perforated sheet No. 1 backup) 
showed zero gas ingestion, and therefore the test with the maximum applied 
pressure differential and maximum flowrate (Run 79) was analyzed to verify 
that zero gas ingestion would be predicted. Similarly, all tests with 200 x 
1400 mesh (with perforated sheet No. 2 backup) showed some startup gas 
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ingestion attributed to leakage. Run 85 ingested only a single bubble on 
startup which did not appear on the oscillograph, and Run 86 showed general 
ingestion, so these two runs were also selected. 

All of the vertical specimens showed no tendency toward gas ingestion except 
the 325 X 2300 mesh specimen with no flowmeter. Runs 100 and 101, run at 
low tank (dewar) pressure, showed no ingestion, while Runs 102, 103, and 
104, run at higher dewar pressure, showed clearly observed and measured 
liquid spill on startup, which was construed as gas ingestion. Therefore, 

Runs 101 and 102 were selected for simulation and comparison. 

For the Freon 114 tests, there were few tests which were acceptable for 
simulation because of the problems of gas overpressure (due to higher density 
and lower surface tension of Freon 114, as discussed previously), combined 
with surging and the need for unambiguous determination of gas ingestion. 

The first three tests (Runs 105 to 107) showed no evidence of spill or oscillo- 
graph flowmeter trace spiking. Since all three runs were similar, the high 
flowrate test (Run 107) was selected for analysis. Nearly every screen/ 
structure specimen was of lower bubble point and did not clearly indicate 
transition from zero gas ingestion to observed ingestion (as did the alcohol 
tests above) since general gas ingestion due to overpressure was the usual 
case. 

For the tests, most tests displayed problems of backflow into the bubble 

trap, and definite surging during outflow. The first test (Run 153) indicated 
a clear boiloff pressure spike at the valve, which also caused the differential 
pressure transducer at the screen to respond and follow the surge. Subsequent 
tests often did not display clear line boiloff pressure surges because the line 
was chilled down to some degree. For the first test in each series, the line 
temperature could be estimated with reasonable accuracy, unless preflow 
for operational reasons had chilled the line. Run 153 was unique in that it 
was the best example of cryogenic line boiling pressure surge, the line tem- 
perature was known, the screen specimen had maximum bubble point and 
reasonable initial gas pressure imposed, and all instrumentation was oper- 
ating properly. Therefore, it was the only LH^ test modeled. 
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The 15 tests selected for the H672 simulation, including the observed initial 
gas pressure imposed on the screen, flowrate, and gas ingestion, are sum- 
marized in Table 15. It is believed that the tests selected represent a com- 
prehensive cross-section of specimens and test conditions that will allow 
definitive correlation of the important aspects of the model for transient 
behavior of screen devices as described earlier. 

Table 15 

TESTS SELECTED FOR H672 SIMULATION 

Experimental Observations 




Initial 

AP 

Flowrate 

Gas Ingestion 

Run No. 

Me sh/ C onfiguration 

N/cm^ 

psi 

Kg/ sec 

lb/ sec 

Quality 

15 

200 X 1400 Plain 

0. 17 

0.25 

0. 10 

0. 22 

0. 012 - 0. 016 

21) 







22 [ 

200 X 1400 Pleated 

0. 14 

0.20 

0. 25 

0. 55 

0 

23 ) 







60 

165 X 800 Screen 
Backup 

0. 11 

0. 16 

0. 12 

0. 26 

0 

61 

165 X 800 Screen 
Backup 

0. 12 

0. 17 

0. 14 

0.3 0 

0. 0015 

73 

500 X 500 Plain 

0. 14 

0. 21 

0. 28 

0.61 

0. 22 

74 

500 X 500 Plain . 

0. 05 

0. 07 

0. 28 

0.61 

0. 099 

79 

200 X 1400 Perf. 
Sheet 

0. 14 

0.21 

0. 08 

0. 18 

0 

85 

200 X 1400 Perf. 
Sheet 

0. 18 

0.26 

0. 05 

0. 10 

0. 148-' 

86 

200 X 1400 Perf. 
Sheet 

0. 18 

0.26 

0. 07 

0. 16 

0. 444* 

101 

325 X 2300 Plain 

0. 07 

0. 10 

No flowmeter 

0 

102 

325 X 2300 Plain 

0. 08 

0. 11 

No flowmeter 

0. 003 - 0. 012 

107 

325 X 2300 Plain 

0. 24 

0.35 

0.28 

0. 62 

0 

153 

325 X 2300 Plain 

0. 03 

0. 04 

No flowmeter 

0. 091 


^Includes leakage 
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SCREEN SPECIMEN MODEL, 

Development of the pertinent screen specimen properties and characteristics 
for input to H672 requires knowledge of certain properties which are not 
well-defined, such as percent open area (PO) for the dutch weave screens, 
and effective screen fraction (f) for composite screen/ structure specimens. 
The structural properties of the composite also depend on the screen fraction 
(f), as well as the structural properties of the components of the screen/ 
structure specimen. For modelling purposes, assumptions were made for 
the f of each type of screen/ structure specimen, depending on their construc- 
tion characteristics. As shown previously in Figure 23, the screens are 
slipped over a support tube, which has three longitudinal struts (0.64 cm or - 
0. 25 in wide by 7. 6 cm or 3. 0 in long), and are then bonded to the support 
tube ends. For the plain screens, it was assumed that the struts would 
obstruct flow through the screen, so the effective screen fraction would be 
the screen area less the strut area, which results in an f of 0.84. However, 
the struts would not materially affect the deflection of the screen (since the 
screen would deflect into the space between the struts) so that steel strut 
properties (modulus of elasticity, etc. ) were ignored and only the screen 
properties were used in evaluating the specimen structural properties. 

The pleated screen specimens were made with 0. 5 cm (0. 2 in) deep pleats 
which increased the effective specimen diameter to 4. 3 cm ( 1. 7 in), rather 
than the 3. 8 cm (1,5 in) diameter of the plain screens. The f for the pleated 
screens was found by taking the actual screen circumference (number of 
pleats times 2, 0 times 0. 5 cm) and dividing by the circumference of a 4. 3-cm 
cylinder. 

The perforated sheet used as screen backup material was stainless steel 
0. 064 cm (0, 025 in) thick with 0.478 cm (0. 188 in) holes giving 50% open 
area. It was assumed that the 50% solid sheet kept half the screen from 
deflecting and obstructed half the flow area of the screen, hence the f was 
equal to 0. 5. The perforated sheet did not contact the three struts mentioned 
above, so that their additional blockage was ignored. In addition, it was 
assumed thkt the structural properties of the screen/sheet composite were 
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made up of a 50% contribution by the screen and a 50% contribution by the 
solid sheet in accordance with equations (3) through (7). (See the section on 
Analytical Program. ) 

For the screen backed up with coarse-mesh screen, two different models 
were assumed: 

1. The coarse mesh would make only point contact with the fine mesh. 

The fine mesh could deflect freely into the spaces of the coarse 
screen, giving f = 1. 0 and assuming that the coarse screen properties 
do not contribute to the composite properties, but only the fine mesh 
screen properties would be used. It was also assumed that the PO 
for the fine mesh screen would be used since the coarse mesh screen 
would not materially affect the flow. 

2. The second model assumed that the coarse screen rigidly supported 
the fine mesh screen, but that the coarse screen did not block the 
flow area of the fine screen, hence the screen fraction f would still 
equal 1. 0. However the coarse screen properties would be combined 
with the fine mesh properties in the ratio of (1-PO) and PO, according 
to equations (3) through (7). For this case of the "rigidly" supported 
fine screen, again only the fine screen PO would be used. 

The two models described above represent the two extremes of screen 
structural/flow behavior: the first quite compliant with minimum tendency 
for gas ingestion, and the second quite stiff with maximum gas ingestion 
tendency. 

The percent open area (PO) was shown in the Analytical Program section to 
be important in affecting screen gas ingestion characteristics, yet it can not 
be analytically defined for dutch weave screens. However, for all square 
weave screens, PO can be analytically defined and can be related to the 
square weave screen void fraction, e, by 

PO = 

where the exponent n has values from 2. 4 to 2. 6. It was assumed that it 
would be conservative to use the same equation to define the value of PO for 
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dutch weave screens based on their e (which is known), since the weaving 
process for the dutch weave screens could tend to reduce the PO relative to 
square-weave screens. The larger PO thus assumed for the dutch weave 
screens would have a greater tendency for gas ingestion, hence would be 
conservative. 

The experimental proportional limit stress, elastic moduli, bubble point, 
and flow loss coefficients were used to define the required inputs for all of 
the screen/ structure specimens selected for analysis that are shown in 
Table 16. Different models for the coarse screen backup specimens would 
be resolved through the results of the data correlation. The bubble points 
and flow loss coefficients (and gas properties) also depend on the test fluid 
properties; those shown in Table 16 are for isopropyl alcohol. All the units 
in Table l6 are British Engineering Units because these were used in the 
H672 code. 

CORRELATION RESULTS 

There are a number of problems in running the H672 code which make corre- 
lation of experimental data difficult. First, the computing time step must be 
kept small (~ 0. 0005 sec) in order to be a reasonable approximation of the 
acoustic wave travel time through the shortest important component. On the 
other hand, the acoustic length of the longest component cannot exceed 40 time 
steps, so the system piping must be broken down into relatively short lengths. 
Because of the short computing time step, obtaining a reasonable time span 
for simulation requires many steps and results in large quantities of data. 
Start up of the flow process (valve opening) is straightforward, but it is not 
practical to start the simulation run with an existing head/flow distribution 
and shut down (valve closing) because the heads and flows everywhere in the 
system must be accurately and consistently specified. Therefore, in order 
to simulate shutdown, the simulation must also include startup, followed by 
adequate run time to achieve reasonable stability. To obtain the best com- 
promise between total run time (and volume of data) startup flow stabilization, 
shutdown, and transients following shutdown, the following arbitrary event 
timing was used; valve starts to open at 0. 001 sec, valve is wide open at 
0.201 sec, valve starts to close at 0.25 sec, valve is completely closed at 
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Table 16 


SCREEN/STRUCTURE SPECIMEN CHARACTERISTICS 

(Alcohol) 


Mesh 

Construction 

BPS'' 

EE 

STEM 

ASl 

BSl 

AS2 

DPS 

DS 

DC 

PO 

f 

325 X 

2300 

Plain 

2. 57 

4473. 

60. 3 

21. 35 

2. 27 

0. 187 

0. 000020 

1.5 

0. 0001035 

0. 027 

0. 84 

200 X 

1400 

Plain 

1. 71 

5753. 

46. 6 

16. 57 

2. 01 

0. 145 

0. 000033 

1. 5 

0. 0000791 

0. 027 

0. 84 

200 X 

1400 

Pleated 

1. 71 

3952. 

46. 6 

16. 57 

2. 01 

0. 145 

0. 000033 

1.7 

0. 0000791 

0. 027 

2.40 

200 X 

1400 

Perf. Sht. 

1. 71 

59430. 

46. 6 

16. 57 

2. 01 

0. 145 

0. 000033 

1.6 

0. 00000582 

0. 027 

0. 50 

1 65 X 

800 

Scrn. Backup 
Model (1 ) 

0. 82 

4114. 

67. 25 

2. 023 

0. 264 

0. 0177 

0. 000069 

1.6 

0. 0000920 

0. 11 

1. 00 

165 X 

800 

Scrn. Backup 
Model (2) 

0. 82 

9006. 

67. 25 

2. 023 

0. 264 

0. 0177 

0. 000069 

1.6 

0. 00003954 

0. 11 

1. 00 

o 

o 

If) 

500 

Plain 

0. 80 

1031. 

12.44 

1.037 

0. 140 

0. 0091 

0. 0000833 

1. 5 

0. 000409 

0. 25 

0. 84 


See Table 3 for nomenclature and units. 

Common to all screen specimens were 

LS = 3.0 in. Gamma = 1.67 RTGAS = (Depends on fluid temperature) 



0.45 sec, and run terminates at 0. 75 sec (the total valve open/close time 
was not arbitrarily selected, but was based on the experimental data). The 
simulations that follow show that the flow often continues to increase and/or 
stabilize even while the valve is closing so that usually more than 0.4 sec is 
available to analyze the startup transient, and about 0, 3 sec to analyze the 
shutdown transient. However, these are rather limited times relative to the 
experiment operation times (~3 sec) or fluid boiling transient times (~1 to 5 
sec depending on fluid), so that a certain amount of phenomenological overlap 
is inevitable. Generally, this problem does not arise; exceptions will be 
discussed below. Three basic parameters were chosen for the data correla- 
tion: (1) the screen-ullage pressure differential, (2) the outflow rate, and 
(3) the gas ingestion quantity (converted to effective quality in the screen 
device. ) Correlations for this first two parameters are presented as plots 
generated by the H672 code and the SC4460 plotter, superimposed on which 
are the experimental data taken from the oscillograph records such as shown 
in Figures 34 to 39. The gas ingestion quantity observed in bubble trap or 
via liquid spill quantity is compared to the H672 prediction as a gross number 
(not as a plotted parameter). 

The first test simulated (Run 15) had experimentally shown modest gas inges- 
tion (equivalent quality, X, of 0. 012 to 0, 016) manifested as a single 1.3-cm 
(0. 5-in) diameter bubble into the bubble trap at startup. The H672 simulation 
for this test, run at the proper flowrate, indicated gross, unstable gas inges- 
tion starting at 0.436 sec, which reached X = 0. 14 at 0.471 sec. The unstable 
pressure surges caused excessive screen stress at 0.464 sec, as shown in 
Figure 41, and a gross flow reversal. The apparent stabilization of the pres- 
sure surges after 0. 53 sec is artificial since the code was no longer computing 
properly due to the gross instabilities experienced earlier. As these insta- 
bilities were not in accord with observation, the screen model was reexamined 

The model had assumed that gas ingestion was uniformly distributed through 
the screen device and affected the acoustic velocity in the screen device 
according to equation (27). It was apparent that uniformly distributed gas 
ingestion was not a physical reality for the screen specimen, because essen- 
tially all gas ingestion would occur at the top downstream edge of the screen 
(where the dynamic /friction/gravity head difference is maximum), and 
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Figure 41. H672 Simulation — Run 15 — la' = f{x)l 
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OUTFLOW RATE - KG PER SEC SCREEN-ULLAGE .iP - N/CM‘ 


quickly would move into the bubble trap (see Figure 42) where it would no 
longer affect the screen or pipe flow or acoustic properties. Therefore, 
the screen model was modified so that the quality increase from gas inges- 
tion no longer affected the acoustic velocity in the screen device. The effects 
of the model change on the screen-ullage pressure difference and flow for 
Run 15 are shown in Figure 43. 

The final quality is X = 0. 014, which is excellent agreement with experiment, 
and the pressure differential is stable and does not result in excessive stress, 
while the flowrate after shutdown behaves in much more reasonable fashion. 
Also shown in Figure 4 3 is the good agreement with starting experimental 
data, which verifies that the assumptions made in deriving PO and f were 
appropriate. The pressure surges occurring after shutdown were a common 
simulation problem which is discussed in detail below. 


CR21 



Figure 42. Bubble Ingestion Phenomenon 
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OUTFLOW RATE - KG PER SEC SCREEN-ULLAGE AP-N/CM' 






For the vertical specimen orientation, where the bubble trap is not used, it 
was conceivable that bubble ingestion could still affect the acoustic velocity 
and pressure/flow relationships in the screen device. Run lOZ with the 
325 X 2300 screen specimen in the vertical orientation, demonstrated a very 
small liquid spill of 0. 25 to 1.0 cc (construed to be equivalent to a gas 
ingestion quality of 0. 003 to 0. 012). When the original model (with screen 
acoustic velocity a function of X) was used to simulate run 102, gas ingestion 
started at 0. 448 sec and continued until the screen stress was encountered 
at 0.461 sec, as shown in Figure 44. Note also in Figure 44 the gross 
discontinuities at about 0. 5 sec. When the modified model was used to sim~ 
ulate run 102, the predicted quality was 0. 003 at 0.448 sec, and the predic- 
ted pressure differential and flowrate are shown in Figure 45. Figure 45 
also shows the agreement with the pressure differential data (flowrate data 
was not taken). 

It was apparent that the screen model in which the quality increase due to 
gas ingestion does not affect the acoustic velocity within the screen device 
was physically more realistic and provided better correlation of the experi- 
mental data, hence this model was used for all further data correlation. 

Simulation of Run 101 (companion run to Run 102) indicated that zero gas 
ingestion occurred, as observed experimentally. The simulations compared 
to experimental data for pressure difference and predicted flowrate are 
shown in Figure 46. 

For the 200 x 1400 mesh pleated specimens (Runs 21, 22, 23), the simulation 

of the startup flow and pressure differential (see Figure 47) was excellent, 

and no gas ingestion was predicted (in agreement with the experimental data) 

until after the shutdown pressure spike arrived at 0.484 sec. Note that 

Figure 47 indicates that significant pressure surging was predicted following 

shutdown. In fact, the data showed no such surging. The observed shutdown 

behavior for virtually all tests was poorly simulated by the H672 code. The 

2 

code generally predicted a very large pressure spike (over 69 N/cm , 
loo psia, in some cases) at the valve immediately after valve closure; this 
was not observed by the pressure transducer at that location. In the half- 
dozen or so cases of waterhammer that were observed, the maximum pres- 
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screen-ullage ap-psi 




Figure 44. H672 Simulation — Run 102 — [a' = f(x)] 
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Figure 45. H672 Simulation — Run 102 
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Figure 46. H672 Simulation — Run 101 
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Figure 47. H672 Simulation — Runs 21, 22, 23 
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sure excursions at the valve were of the order of only 3. 5 N/cm (5 T>si); 

(see Figures 38 and 39) and the' response of the screen pressure transducers 
to these modest excursions were at most about 0, 07 N/cm (0. 1 psi). 

Despite considerable effort investigating the reasons for this poor simulation, 
its cause remains unknown; possible, but unconfirmed, reasons include valve 
behavior at closing not in accordance with the valve model, gas bubbles in 
the outlet line, and leakage through the valve at closure, 

"When the large pulses generated by the valve closure reach the screen device, 
they are considerably attenuated (see Figure 47) but cause gas ingestion 
(negative pulse) and liquid/gas outflow (positive pulses). A potential short- 
coming of the screen model is that liquid/gas outflow from a positive pulse 
in excess of the bubble point is in the ratio of the average quality in the 
screen device, assuming uniform distribution of the gas within the device. 
This is probably not precise; the quality near the screen may be much higher 
than the average quality since the gas bubbles ingested may stay near (or 
attached to) the screen. This is especially true for rapidly oscillating plus 
and minus pressure pulses, since the bubbles would not have time to 
migrate far from the screen. Therefore, the time required for the screen 
device to pump itself full of gas (while spilling liquid) may not be accurately 
predicted. In our tests, the observed pressure excursions were quite low; 
therefore, there were almost no cases of observed liquid spill following 
shutdown. 

The simulations for the 165 x 800 mesh with coarse screen backup (Runs 60 

and 61) were able to resolve the choice of the proper model for the properties 

of this form of construction (see previous discussion in section on Screen 

Specimen Model). For Run 60, where no gas ingestion was observed, the 

more compliant model (assuming that only the 165 x 800 mesh screen 

properties were appropriate) correctly predicted zero gas ingestion. For 

Run 61, where gas ingestion resulting in an effective quality of 0. 0015 to 

0.003 was observed, the same model also predicted zero gas ingestion, but 

2 

the pressure pulses came within a few thousandths of a N/cm (psi) of the 
bubble point. "When the initial gas pressure around the screen was increased 
by only 0. 007 N/cm (0. 01 psi), well within the possible error in recording 
the gas pressure, the result was gas ingestion which reached a quality of 
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0.089 by 0.479 sec when the shutdown pulse arrived. The less compliant 

model (assuming that the 165 x 800 mesh screen properties were combined 

with the backup screen properties) properly predicted gas ingestion (to a 

quality of 0. 075 at 0.479 sec) for Run 61, but also predicted gas ingestion 

for Run 60, which had none. Reduction of the initial surrounding gas pressure 
2 

by 0. 007 N/cm (0. 01 psi) reduced the quantity of gas ingestion predicted 
for Run 60, but ingestion was still predicted. Therefore, it is believed that 
the compliant model for the screen backup structure is the model which best 
represents the response of this type of composite screen structure. Fig- 
ures 48 and 49 indicate the good agreement with pressure differential and 
flow data for the simulations of Runs 60 and 61 respectively. 

The simulation for Run 74 for the 500 x 500 mesh specimen also indicated 

2 

that at 0. 05 N/cm (0. 07 psi) initially imposed gas pressure, there would be 

zero gas ingestion, but that if the initial gas pressure were increased by only 
2 

0.007 N/cm (0.01 psi), gas ingestion would occur. The predicted pressure 
differential and flowrate for Run 74 is compared to the data in Figure 50. 

Note from Figure 50 that the flowrate is somewhat underpredicted (at these 
high flows), which may account for discrepancy in gas ingestion prediction. 
The simulation for Run 73, showing general gas ingestion, as observed, is 
shown in Figure 51. 


The final alcohol runs simulated were for the 200 x 1400 mesh backed up by 
perforated sheet (Runs 79, 85, 86). For the tests in which no gas ingestion 
was predicted (Runs 79 and 85), the pressure differential was somewhat 
^ overpredicted but the flowrate prediction was accurate (see Figures 52 and 53). 
For Run 86, with general gas ingestion observed, the pressure differential 
prediction was accurate until gas ingestion occurred (see Figure 54) when 
both the pressure and flow observed were erratic as was often the case for 
tests with general gas ingestion. 

The simulation for the selected Freon 114 test (Run 107) with the 325 x 2300 

mesh screen in the vertical orientation was started by using the P4557 code 

to predict the magnitude of the boiling pressure surge in the line. The 

observed startup pressure surges recorded by the pressure transducer next 

2 

to the valve, were of the order of only 2 N/cm (^3 psi) or less. However, 
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Figure 48. H672 Simulation — Run 60 
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Figure 49. H672 Simulation — Run 61 
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Figure 50. H672 Simulation — Run 74 
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Figure Si. H672 Simulation — Run 73 
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the P4557 code did not predict this large a pressure pulse in the time 
observed (-'^0. 1 sec), but rather a slow rise in pressure toward saturation 
pressure at the line temperature. It is believed that the pressure pulses 
observed were hydrodynamic in nature, resulting from liquid flow into a 
partially empty line, followed by a sudden stop at the restrictive throttle 
valve just downstream from the ball valve. In general, the recorded pressure 
differential at the screen device did not respond to the pulse recorded at 
the valve. The hydrodynamic simulation of Run 107, using H672 without the 
downstream pressure pulse input through the injector /chamber pressure 
spike routine, is shown in Figure 55. Both the pressure differential and 
flow are somewhat overpredicted, yet zero startup gas ingestion was 
predicted as observed. 

Simulation of the ^^2 was started with modeling of the LH^ cryosurge 

using the cryosurge code P4557. The modeling was complicated by the dis- 
crepancy between the P455 7 code assumption of a long horizontal line 
downstream of the control valve, and the actual case of a high-resistance 
throttle valve downstream of the control valve. The long-line equivalent 
of the throttle valve was developed from P4557 simulations of the Freon 114 
tests with the flowmeter installed, matching the predicted flowrate to the data. 

The equivalent line to the L 1 H 2 configuration was then corrected by removing 
the flowmeter resistance and reconfiguring the line. The resulting peak 
flowrate predicted by the P4557 code was close to the value predicted by the 
H672 code when it was run without the cyrosurge option (in order to determine 
the system hydrodynamic flow conditions). The experimental data for the 
initial valve pressure spike and differential pressure for Run 153, as traced 
from the oscillograph record, is shown in Figure 56. The circled points 
shown in Figure 56 are the interface (valve) pressure trace predicted by 
P4557. The peak value matches well, although the predicted rise time is 
somewhat slower because of the discrepancy in modeling the physical 
restriction of the throttle valve. The pressure pulse data was input to the 
H672 code through the injector /chamber pressure spike option, and the 
H672 simulation is shown in Figure 57 for screen ullage pressure differential 
and predicted flow. The reverse flow due to the pressure spike is evident 
from Figure 57. Gas ingestion occurs rather generally until the reverse flow 
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Figure 55. H672 Simulation — Run 107 
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Figure 56. LH 2 Cryosurge Data 



at 0. 37 sec causes outflow from the screen device and cessation of ingestion, 
until the pressure drops below the imposed gas pressure again at 0. 50 sec. 
The screen ullage pressure differential peaks from Figure 56 are super- 
imposed on the H672 prediction in Figure 57 where there is fair qualitative 
agreement. 

The results of all of the data correlation simulations are summarized in 
Table 17, which compares data and predictions for gas ingestion flow and 
screen-ullage pressure differential, and comments on the aspects of the 
screen model verified by the correlations. 

The results of Table 17 indicate that the H672 code properly simulates flow 
startup and operational characteristics of a feed system, but overpredicts 
the shutdown pressure pulses (compared to experimental data). The SR26 
subrouting modelling the screen device was shown to be an adequate 
representation of the characteristics of a screen device in response to 
pressure pulses regarding screen-ullage pressure differential, flow, and 
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Figure 57 . H 672 Simulation — Run 153 





SUMMARY COMPARISON OF ANALYTICAL CORRELATIONS AND DATA 
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gas ingestion. The model of the screen/structure which includes screen 
properties and structural backup material properties as a uniformly- 
distributed composite structure was verified for plain screen specimens,, 
pleated screen specimens, and screen backed up by coarse screen or 
perforated sheet. The results indicate that pleated screen construction has 
maximum performance, and screen/structure methods which use only the 
more compliant screen properties (plain, pleated, coarse screen backup) 
are superior in performance to perforated sheet backup in which the less 
compliant backup material properties are included. The P455 7 code 
adequately simulated the cryosurge and, when modeled with the 

H672 code, produced fair agreement with experiment. 

It is believed that the H672 code with screen model has been experimentally 
verified as a useful design tool which can predict the startup transient 
performance of screen devices under a wide variety of fluids and operating 
conditions. 


' 'i ; • ; ■' 
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CONCLUSIONS 

The conclusions reached as a result of this comprehensive analytical and 
experimental program to provide an analytical design tool, verified by 
experiment, to predict the effects of transient pressure surge on the reten- 
tion performance of screen acquisition systems are as follows: 

1. A detailed survey of current and potential space vehicle screen 
acquisition systems defined important operational and configura- 
tiona;! aspects of twelve systems and defined appropriate categories 
of screen systems as (1) localized— for engine restart, and (2) dis- 
tributed— for long term transfer. 

2. An existing MDAC computer code, H672, for feed system transient 
analysis was adapted for study of screen device performance by the 
addition of screen model subroutines. The MDAC cryosurge anal- 
ysis, P4557, was used to analyze pressure surges caused by cryogen 
boiling in warm lines, and provide input conditions for the H672 
code. 

3. The screen device model developed and programmed for inclusion 
into the H672 code had the following features: 

a. The screen device was treated as a uniformly-distributed 
composite structure incorporating orthotropic screen properties 
combined with backup structure properties according to the 
relative screen area fraction, f. 

b. The screen device flow characteristics depend on the percent 
open area, PO, assumed as a power function of the void frac- 
tion, and experimental flow loss characteristics as defined in 
Reference 1. 

c. The screen device model provides for liquid inflow /outflow; 
liquid leaving the screen; gas ingestion; and when resisted by 
the screen bubble point, no gas ingestion; screen stress analysis 
and quality in the screen device due to gas ingestion. 

d. The acoustic velocity and pressure/flow balance in the screen 
device was assumed as a function of gas ingestion and quality. 
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A comprehensive experimental program of 177 tests using 13 screen 
specimen configurations; 5 meshes, 4 construction methods, and 
2 orientations, together with 3 fluids; alcohol, Freon 114 and LH^, 
verified the following aspects of the screen model by good agree- 
ment between analysis and experiment for gas ingestion, screen- 
ullage pressure difference, and flowrate: 

a. The assumption of screen device properties as a composite 
of screen and backup material properties was confirmed for 
the following structures: 

• plain screen; f = 0. 84 (due to strut blockage); screen 
properties only 

• pleated screen; f = 2. 4 (actual screen area); screen 
properties only 

• screen with coarse screen backup; f = 1.0; fine screen 
properties only 

• screen with perforated sheet backup (50% open area); 
f = 0. 5; combined screen properties and perforated 
sheet properties 

b. The assumption of screen device flow characteristics as a 
function of the assumed percent open area functional relation- 
ship was confirmed. 

c. The screen model operational characteristics of paragraph 
(3c) above were confirmed. 

d. The acoustic velocity and pressure/flow balance in the screen 
device was shown not to be a function of quality, and the screen 
model was modified accordingly. 

e. The cryosurge code P4557, was shown to predict the line 
boiling pressure surge with reasonable accuracy. 

The prediction of the feed system flow shutdown transient by the 
H672 code did not agree well with experiment in that it tended to 
overpredict system pressure surges. 

The H672 transient analysis, including the screen model, has been 
experimentally verified to be a useful design tool, which can 
accurately predict the startup transient performance of screen 
acquisition systems under a wide variety of fluids and operating 
conditions. 
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Appendix A 

SCREEN ACQUISITION SYSTEM 
SURVEY DATA TABULATION 
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No. 1 Application 

: Cryogenic Tug 

Name: . Cryogenic Space Shuttle/Auxiliary Propulsion System 

Fluid: 

LH2 

(Saturated at 17.2 N/cnr (25 psia)) 

Quantity: 

4750 KG (10,444 Lb) 

Flowrate: 

0-4.54 KG/Sec 
(0-10.0 Lb/Sec) 

Environment: 

Acceleration: 

"g's" 

On-orbit ± 0.045 (X, Y, Z) 
Reentry 1.25 (Z) 

0.08 (X) 

Boost 3.00 (X) 


\ 


LO2 

(Saturated at 17.2 N/cm^ {25 psia)) 

17.950 KG (39,536 Lb) 

0-22.7 KG/Sec 
(0-50.0 Lb/Sec) 


On-orbit ± 0.045 (X, Y, Z) 
Reentry 1.25 (Z) 

0.08 (X) 

Boost 3.00 (X) 


Vibration: 


U* 


U 


Thermal : 


ML! Tank/Line Insulation 
Cold He Pressurization 


MLI Tank/Line Insulation 
Cold He Pressurization 


Valve/Pump See Reference 

Characteristics: 


u 


Tankage and 
Acquisition System 
Configuration 

Transfer Line 
Configuration 


See Figures A-1, A-2 


10.2 cm dia x .05-cm 
wall X 3.05 m long 
(4-in. dia x .02-in. 
wall X 10 ft long) 
Aluminum 


Similar to Figures A-1, A-2 except 
17.8-cm (7-in.) diameter channel 


7.6 cm dia x .05-cm 
wall X 1.53 m long 
(3-in. dia x .02-in. 
wall X 5 ft long) 

A1 uminum 


Reference : G. W. Burge and J. B. Blackmon. Study and Design of Cryogenic Propellant Acquisition 

Systems - Volume I, Design Studies, MDAC Report MDC G5038, Deceirfcer 1973 


* U - Unknown 




FINE MESH STEEL 



Figure A-2. LHj, CHANNEL CONSTRUCTION DETAIL 
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No, 2 Application : Cryogenic Tug 

Name : Advanced Space Propulsion Module 


Fluid: 


Quantity : 

Flowrate: 


Envi ronment : 

Acceleration: 

IlgIgH 

Vibration: 
Thermal ; 


LH2 

(Saturated at 14.7 H/cftr (21.3 psia)) 

3626 KG (7,972 Lb) 

1.38 KG/Sec 
(3.04 Lb/Sec) 

0.0001-1.64 

U 

MLI Tank/Line Insulation 
Cold He Pressurization 


L02 , 

(Saturated at 15.9 N/cnr (23 psia)) 

21,118 KG (46,428 Lb) 

8.28 KG/Sec 
(18.24 Lb/Sec) 

0.0001-1.64 

U 

MLI Tank/Line Insulation 
Cold He Pressurization 


Valve/Pump U 

Characteristics: 


U 


Tankage and See Figure A-3 

Acquisition System 

Configuration 

Transfer Line U 

Confi guration 


See Figure A-4 


U 


Reference : G. W. Burge and J. B. Blackmon. Study and Design of Cryogenic Propellant Acquisition 

Systems - Volume I, Design Studies, MDAC Report MDC G5038, Decetrfjer 1973 
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OOMON DOME 


. FLANGE 
■ PRESS 6 VDfT 

- LEVEL SENSOR - CAPACITANCE 


HIGH PRESSURIZATION HELIUM 
(ZOOO PSIA) ) 


y// <1-3 R 

' A 

TOP REGION (16.8 FEET^) 


"SCREEN TUBE (PLEATED) (8-TYPICAL) 

^ plate 

SCREEN (PLEATED) 


SCREEN TUBE (PLEATED) (8-TYPICAL) 
^ BAFFLE 


TO LO^ START TANK 


SOLENOID VALVE 


PRESS 6 VENT 
COMPONENTS 


PRIMARY TRAP REGION W|| 
(10.8 FEET*) |:,j| 


VISOOJET ^ 
- VENT VALVE 


, SECONDARY TRAP 
REGION 


-12 IN. R - 


I II FT3 j 
r 6.S FT3 ' 


r 

JOcN\ / 


' ooc; \ 

I! 

90 C' 


9o:: 




SCREEN TUBE (PLEATED) 
PORT FOR SCREEN TUBE 

— MAHHAN FLANGE 
FEED LINE 


REVILL VALVE 


Figure A-3. ASPM LH 2 ACQUISITION SYSTEM PRELIMINARY DESIGN 
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No. 3 Application ; Cryogenic Tug 


Name ; Centaur D-IS 

Fluid ; LH2 

(Saturated at 14.1 N/cnr (20.5 psia)) 


LO2 

(Saturated at 21.4 N/cm^ (31.0 


Quantity ; 

Flowrate ; 

Envi ronment ; 

Acceleration; 

"g's" 

Vibration; 

Thermal; 


Valve/Pump 

Characteristics; 


2397 KG (5279 Lb) 

5.1 Kg/Sec 
(11.2 Lb/Sec) 


See Reference 


See Reference 

Common Bulkhead Heating Rate 
+19.5 W/fi2 (+6.2 B/Hr-Ft2) 
Warm Helium Pressurization 

Initial Sump/Pump Chilldoi/n 
See Figure A-5 

Boost Pump Startup 
See Figure A-7 

Engine Chilldovm and Shutdown 
See Figure A-9 


Tankage and 
Acquisition System 
Gonfi guration; 


See Detailed Drawings in 
the Reference 


Transfer Line 
Configuration ; 


11554 KG (25450 Lb) 

25.6 Kg/Sec 
(56.4 Lb/Sec) 


I See Reference 

I See Reference 

I 

I -19.5 W/m 2 (-6.2 B/Hr-Ft^) 

I Bubbled Helium Pressurization 

I 

. Initial Sump/Pump Chi 11 down 
I See Figure A- 6 

Boost Pump Startup 
See Figure A-8 

Engine Chilldov/n and Shutdown 
See Figure A-10 

! 

i 

I See Detailed Drawings in 

• the Reference 


Reference; M. H. Blatt and M. D. Walter. Centaur Propellant Acquisition System Study. 
GD/C Report CASD-ilAS-75-023 (NASA CR-134811), June 1975. 


psia)) 
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I 

oi 

3 

o 

rJ 

s 

J 



T\vAe- sec. 


Figure A- 7* LHg BOOST PUMP STARTUP FLOW HISTORY 



Figure A-S. LO 2 BOOST PUMP STARTUP FLOW HISTORY 


V42 


fcic./sec tc»/sec 




Figure A-9. EMGIflE CHILLDOWN/SHUTUOWN LH 2 FLOW HISTORY 



Figure A-10. ENGINE QIILLDOWN/SHUTDOWN LO 2 FLOW HISTORY 
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t-Ca/se-c. icc»/sec. 



No. 4 Application ; Cryogenic Tug 

Name ; Interface Demonstration Unit (IDU) 

Fluid ; LH 2 

(Saturated at 20.7 N/crr (30 psia)) 

Quantity ; 19 KG (42 Lb) 

Flowrate; 3.4 KG/Sec (7.5 Lb/Sec) 


Envi ronment ; 

Acceleration; -1.0 (X) 

II . I ^ II 

g s 

Vibration; U* 

Thermal; Cold and Warm, Overhead and Submerged, 

GHe and GH£ Pressurization 


Valve/Pump 

Characteristics ; U - Valve Malfunctions 

Tankage and See Reference 

Acquisition System 

Configuration; 

Transfer Line See Reference 

Configuration ; 


Reference ; J. B. Blackmon. Design, Fabrication, Assembly, and Test 
of a Liquid Hydrogen Acquisition Subsystem. MDAC Report 
MDC G5360. May 1974. 


* U - Unknown 
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No. 5 Application ; Cryogenic Tanker 


Name : S-IIB LO 2 Tanker 

Fluid ; 

Quanti ty ; 

Flowrate: 


LO2 

(Saturated at 10.3 N/cm^ (15 psia)) 
90,720 KG (200.000 Lb) 

16.8 Kg/Sec (37 Lb/Sec) 


Environment ; 

Acceleration: 

Vibration; 
Thermal ; 


Drag: -10"^ (X) 

Disturbance; ±10“^ (X, YZ) 

U 

Cold Pressurization with Thermodynamic Vent 


Valve/Pump 

Characteristics; 


Tankage and 
Acquisition System 
Configuration; 

Transfer Line 
Confi qu rati on; 


U 


See Figure A-11 


U 


Reference ; M. H. Blatt, et. al. Low-Gravity Propellant Control Using 

Capillary Devices in Large-Scale Cryogenic Vehicles - Phase I, 
Final Report, NAS8-21465. GDC Report GDC- DDB 70-007, August 1970 



ACCESS OIENING 
TRANSITION SECTION 


15 GORE 
PANEIS 

PANELS 

INTERCONNECTED 
AT CORE LINES 
AND AT EQUATOR ’ 


WELD CLEVIS 
. STRIP TO 
TUBE 


WELD 
WEB TO 

tube 




.35 SPACING 
BETWEEN IJNER 
AND TANK WALL 


2?.0 INCH DIA 
SPHERICAL . 
TANK 


SUPPORT GORE 
FRAMES FROM 
TANK WALL 


OUTLET TRANSITION 
SECTION 


CRES TUBE 
FRAME 

SEE ALTERNATE ARRAN®MENT 


SEAM WELD OR 
BRAZE SCREENS 
TO WEB 

■CAPILLARY AND 
BACK UP SCREEN 


VIEW AA 

CLEVIS STRIP 
CONTINUOUS lEG ON F.S. 
AND LOCAL TABS ON N.S. 


WEB STRIP 


BUTT WELD LEG OF 
CIEVIS STRIP TO 
WEB @ ASSEMBLY 



try y y y w y y 



PRIMARY 

TUBULAR 

frame 


SECONDARY 

TUBUIAR 

SUPPORTS 


CAPILLARY 
SCREEN AND 
BACK-UP MESH. 
WELD OR BRAZE 
TO frame 


TYPICAL GORE EIEMEKT 


TYPICAL CONNECTION BETWEEN PANEIB 


Figure A-11* Oxidizer Collector System/Full Liner Concept 
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No. 6 


Application : Cryogenic Tanker 

Mare : Tug-Scale Transfer Module 

Fluid : LH 2 - 

(Saturated at 17.3 H/cnr (25 psia)) 

Quantity : 4600 KG (10,140 Lb) 

Flowrate : .08 KG/Sec 

(.175 Lb/Sec) 


Envi ronrent : 

Acceleration: 

”g's" 

Vibration: 
Thermal ; 


Valve /Pump 
Characteristics: 


Coast: -10"^ (X) 


U 

Cold GH 2 Pump Pressurization 
H 2 Thermodynamic Vent 

U 


Tankage and 
Acquisition System 
Confiqurati on: 


L/D = 2 Tank with Hemispherical 
Ends 

D = 3.78M (12.4 Ft) 

See Figure A- 12 


Transfer Line U 

Configuration ; (3.81 cm dia. x 7. IIM long) 

(1.5 in dia. x 282 in long) 


LO2 

(Saturated at 19.3 N/cm^ (28 psia)) 

22,300 KG (49,160 Lb) 

.39 KG/Sec 
(.86 Lb/Sec) 

-10"^ (X) 

U 

Cold GO 2 Pump Pressurization 
H 2 Cooling of LO 2 Tank 

U 


Spherical Tank 

D = 3.44M (11.3 Ft) 

See Figure A- 12 

U 

(3.81 cm dia. x 2.03M long) 
(1.5 in dia. x 80 in long) 


Reference ; E. C. Cady. Design and Evaluation of Thermodynamic Vent/Screen 
Baffle Cryogenic Storage System. HDAC Report MDC G5479 
(NASA CR-1 34810) June 1975. 



PRESSURE vessel 




Figure A-12. SCREEN LINER MOUNTING METHOD 
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Ho. 7 Application ; Cryogenic Stage 


Hame: S-IVC Stage Restart 


Fluid : 
Quantity : 
FI ov;rate ; 


LH2 

(Saturated at 17.2 H/cm^ (25 psia)) 

19,500 KG (43,000 Lb) (Full) 
[Restart at 60" Full] 

1.27 Kg/Sec (2.8 Lb/Sec) to 
33.7 Kg/Sec (74.4 Lb/Sec) 


Envi ronment : 

Acceleration: 

“g's" 

Vibration; 
Thermal : 


Valve/Pump 

Characteristics: 


Drag: -10“^ (X) 

Idle Mode: .00823 (X) 

Full Thrust: .337 (X) 

U 

100*K (180*R) GH 2 Pressurization 

Orbital Heating 
627 W (2140 B/Hr) 

See Reference 


Tankage and See Figure A-13 

Acquisition System 

Configuration; 

Transfer Line See Reference 

Configuration ; 


LO2 

(Saturated at 17.2 H/cm^ (25 psia)) 

88,900 KG (196,000 Lb) (Full) 
[Restart at 60% Full] 

3.7 KG/Sec (8.1 Lb/Sec) to 

168.7 KG/Sec (372 Lb/Sec) 


Drag: -10”^ (X) 

Idle Mode: .00823 (X) 

Full Thrust: .337 (X) 

U 

150“K-200‘*K (270'’R-360'»R) 
GHe Pressurization 

Cooling: -381W (1300 B/Hr) 

See Reference 


See Figure A- 14 


See Reference 


Reference ; M. H. Blatt, et. al. Low Gravity Propellant Control Using Capillary 
Devices in Large-Scale Cryogenic Vehicles - Phase II, Final Report, 
NAS8-21465. GDC Report GDC- DDB 70-008, August 1970. 
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Figure A- 13. Collector Tube Details 
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MAT'L & ALLOWABLES 
2219 - T62 
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MOUNTED ON OUTSIDE FACE 



Figure A-14'. Oxidizer Tank Reservoir General Arrangement 
and Design Ground Rules 
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No. 8 Application ; Space Shuttle 
Name; Orbital Maneuvering System (OMS) 


Fluid ; 

Quantity ; 

Flov/rate; 


Envi ronment ; 

Acceleration; 

"glgll 


Vibration; 
Thermal ; 


Valve/Pump 

Characteristics; 


Tankage and 
Acquisition System 
Configuration; 

Transfer Line 
Confi guration; 


3688 KG (8130 Lb) 

Startup; 5.41 KG/Sec (11.93 Lb/Sec) 
Settled (Max); 10.82 KG/Sec 
(23.86 Lb/Sec) 

RCS Supply; 6.35 KG/Sec (14.00 Lb/Sec) 


Boost; +3.3 (X), =0.3 (Z) 

Coast; +.0767 (X), +.04 (X.Y) 
-.0230 (X). (Roll, Yaw) 
Settling; +.024 (X) 

Entry; -1.6 (X), +2.5 (Z) 

See Figure A-15 

Propellant Temperature: 

4.4°C - 37.8®C (40®F - 100°F) 

See Figure A- 16 
(Plus Additional Data) 

See Figure A- 17 

(Plus Complete Detail Drawings*) 


Complete Detailed Drawings Available* 


MMH 

2221 KG (4896 Lb) 

3.28 KG/Sec (7.23 Lb/Sec) 
6.56 KG/Sec (14.46 Lb /Sec) 

3.97 KG/Sec (8.75 Lb/Sec) 


I +3.3 (X). -0.3 (Z) 

I +.0767 (X), +.04 (X.Y) 

I -.0230 (X), (Roll, Yaw) 

+.024 (X) 

i -1.6 (X), +2.5 (Z) 

See Figure A-15 


4.4°C - 37.8“C (40°F -lOO^F) 

See Figure A- 16 
\ (Plus Additional Data) 

See Figure A- 17 

(Plus Complete Detail Drawings*) 

I 

i Complete Detailed Drawings Available* 


Reference ; D. A. Hess and G. F. Orton; Space Shuttle OMS Propellant Acquisition. 


* MDAC-E Drawings 73T740000, 73T740010, 730501001, 730551001, 730551002, 730561001, 73A000014 
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Figure A-16. GALLERY TRANSIENT PRESSURE AND FLOWRATE 
DURING OMS STARTUP 
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Figure A-17. OMS PROPELLANT ACQUISITION SYSTEM 
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No. 9 Application ; Space Sliuttle 
Name ; Orbiter Reaction Control System 


Fluid: 




MMH 


Quantity ; 609 KG (1343 Lb) 

Flowrate : Per Thruster 

.907 KG/Sec (2 Lb/Sec) 


Envi ronment : 

Acceleration; 

"g's" 


Vibration; 
Thermal : 


Boost; +3.3 (X), -0.3 (Z) 

Coast: +.0767 U) . ±0.4 (X.Y) 
-.0230 (X), (Roll, Yaw) 
Settling; +.024 (X) 

Entry; -1.6 (X), +2.5 (Z) 

See Figure A-15 

Propellant Temperature 

AA% - 37.8‘’C (40«F - 100“F) 


381 KG (840 Lb) 

Per Thruster 

.567 KG/Sec (1.25 Lb /Sec) 


+3.3 (X). -0.3 (Z) 
+.0767 (X), ±.04 (X.Y) 
-.0230 (X), (Roll, Yaw) 
+.024 (X) 

-1.6 (X). +2.5 (Z) 

See Figure A-15 


4.4“C - 37.8“C (40“F - 100“F) 


Valve/Pump U* U* 

Characteristics; \ 


Tankage and See Figure A-18 

Acquisition System 
Confi quration; 

Transfer Line U* 

Configuration : 


I See Figure A-18 



Reference: SS/RCS Surface Tension Propellant Acquisition/Expulsion Tankage Technology. 

MMC Report No. MCR-75-171, August 1975. 

RCS Propellant Tank Design and Performance Definition. RI Supplementary 
Engineering Memo No, MDAC-075-104, 5 September 1975. 


* Not Fully Defined to Date 
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Ho. 10 Application ; Space Shuttle 


Name: Advanced Fuel 

Cell Reactant Supply System 

Fluid; 

LH2 - 

(Saturated at 20.7 N/cnr (30 psia)) 

Quantity: 

(Per Tank) 

41.7 KG (92.0 Lb) 

Flowrate: 

.032 to 1.48 Kg/Hr (.07 to 3.27 Lb/Hr) 


(Saturated at 20.7 N^cm^ (30 psia)) 

3543 KR (781.0 Lb) 

.28 to 10.34 Kg/Hr (.618 to 22.8 Lb/HR) 


Envi ronment : 

Acceleration: 

"g's" 


Vibration: 
Thermal : 


Valve/Pump 
Character! sti cs: 


Tankage and 
Acquisition System 
Configuration; 

Transfer Line 
Configuration ; 


Boost: +3.3 (X), -0.2 (Z) 

Coast: +.0767 U). ± .04 (X.Y) 
-.0230 (X). (Roll. Yaw) 
Settling; +.024 (X) 

Entry: -1.6 (X), +2.5 (Z) 

See Figure A-15 

Cold GH2 Pressurization 
Line Tenperature 20®K - 222*K 

(36*’R - 400° R) 

Not Applicable 


See Reference 

With Full Pleated Liner 


.64 cm Dia. x .056 cm wall x 15.24M long | 
(.25 in dia. x .022 in wall x 50 ft long)' 
Aluminum ' 

( 


+3.3 (X), -0.3 (Z) 

+.0767 (X). +.04 (X.Y) 

-.0230 (X), (Roll, Yaw) 

+.024 (X) 

-1.6 (X), +2.5 (Z) 

See Figure A-15 

Cold GO? Pressurization 
Line Temperature 98° K - 222° K 

(176°R - 400° R) 

Not Applicable 


See Reference 

With Full Pleated Liner 


.64 cm Dia. x .056 cm wall x 15.24M long 
(.25 in Dia. x .022 in wall x 50 ft long) 
Aluminum 


Reference ; E. C. Cady. Design and Evaluation of Thermodynamic Vent/Screen Baffle 

Cryogenic Storage System. MDAC Report MDC G5979 (NASA CR- 134810) June 1975. 
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jjo. 11 Application : Improved Agena 

Name : Primary 5 Secondary Propulsion System 


Fluid: A-bO 

Quantity : 1G33 KG (3600 Lb) 


2858 KG (6300 Lb) 


Flowrate; 


Envi ronment : 

Acceleration: 

"g's" 

Vibration: 

Thermal 

Valve/Pump 
Characteristi cs: 


Tankage and 
Acquisition System 
Configuration: 

Transfer Line 
Configuration : 


Primary 9.1 KG/Sec (20 Lb/Sec) 
Secondary .091 KG/Sec (.20 Lb /Sec) 

15.9 KG/Sec (35 Lb/Sec) 
.159 KG/Sec (.35 Lb/Sec) 

. 

Boost +9.0 (X) 
Coast -0.012 (X) 

Boost +9.0 (X) 
Coast -0.012 (X) 

u 

U 

U 

U 

Mod. 8096 Pump - U 
Pneumatic Belloivs Pump - U 

Mod. 8096 Pump - U 
Pneumatic Bellows Pump - U 

See Reference 

See Reference 

See Reference 

u 


Reference; W. Heller and E. A, Cadv/allader. Positive Expulsion. CPIA Publication 210. 
Hay 1971. 

R. 0. Sloma. Capillary Propellant Management System For Large Tank Orbital 
Propulsion Systems. CPIA Publication 190, Septentier 1969. 
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No. 12 Application ; Satellite 
Name; Orbit Adjust System 


Fluid: 

N 2 H 4 

Quantity: 

^1090 KG (2400 Lb) 

Flowrate: 

.25 to .47 KG/Sec (.55 to 1.032 Lb/Sec) 

Envi ronment: 


Acceleration; 

"g's" 

+.02 (X) + .00013 (Y.Z) 

-.000044 (X) 

Vibration: 

(shock) 

± .005 g-sec (X,Y.Z) 

Thermal ; 

Propellant Temperature 


4.4®C - 48.9°C (40°F - IZO^F) 

Valve/Pump 

Characteristics: 

See Figure A«19 

Tankage and 
Acquisition System 
Configuration; 

See Figure A-19 

Transfer Line 
Configuration: 

See Figure A-19 

Reference; S. C. 

DeBrock. Development and Flight Experience With A 


Capillary Propellant Management System For A Three- Axis 
Stabilized Vehicle. 


S. C, DeBrock, et. al. A Survey of Current Developments In 
Surface Tension Devices For Propellant Acquisition. AIAA 
Paper No. 70-685. June 1970. 
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Figure A-19 


SATELLITE CONTROL SECTION ORBIT ADJUST SYSTEM 



Appendix B 

TECHNICAL DESCRIPTION 
LIQUID PROPULSION FEED SYSTEM 
DYNAMIC ANALYSIS PROGRAM, H67 2 

INTRODUCTION 

The description of propellant feed system and rocket engine operation during 
severe transients inherent in startup, shutdown, pulse operation, or 
throttling modes has long resisted successful analytical attack. Essentially 
the difficulty lay in the nonlinearity of component operation during these 
periods, coupled with feedback effects of the individual components on the 
system as a whole. 

Previous analytical efforts to describe engine transients are of two types. 

The first consists of linearizing the pertinent fluid dynamic equations 
and the relations that describe individual components such as valves and 
pumps. The linearized differential expressions are then transformed by 
La Place relations to a system of simultaneous equations. Solutions of these 
relations yield stability polars, indicating areas of operation for which the 
engine system may operate in a stable mode. Such solutions include com- 
ponent feedback effects, and generally yield a satisfactory picture of the ‘ 
systems operation about its design point. Unfortunately, the linear restraint 
placed on the pertinent equations restricts the method to small perturbations 
about the engines steady- state point. 

The second technique does not resort to linearization of the equations. 
Nonlinear properties of valves, turbines, and combusters are input in 
curve form, and the entire system is iterated until convergence is achieved. 
The engine system is treated for each time interval in a quasi steady-state 
manner. Valve resistances are evaluated at their current degree of opening 
and turbine and pump torques at their instantaneous speeds; utilizing these 
values the engine is steady- state balanced for this particular juncture in time. 
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Feedback effects are of necessity ignored, and are compensated for to some 
extent by continually modifying line component descriptions with observed 
test data. Granted the availability of sufficient experimental feedback, 
component modifications can be formulated which will effect an accurate 
simulation of the engine system. While this procedure can achieve a trust- 
worthy simulation of a: particular engine system, the results are not trans- 
ferrable to new engine systems, nor can the effects of altering component 
parameters sensitive to liquid feedback effects such as surge tanks, valves, 
and injector, be assessed accurately without further experimental tests. 

This appendix describes a new approach. Nonlinear description of engine 
components is retained at discrete junctions; however, the flow equations 
between components are solved by an exact technique, the method of 
characteristics. The result is the formulation of a pseudo wave theory. 

In this analysis, liquid rocket feed networks are portrayed as variable 
sets of components located at discrete junctions within the system. Pressure 
and flow information is transmitted from component to component by 
characteristic waves, generated and appended by perturbations and boundary 
conditions at the components. 

The momentum and continuity equations describing the one- dimensional flow 
of liquids in an elastic conduit are derived and the characteristics solutions 
to the differential equations are developed with modifications for structural 
dynamic effects. These equations are then set up in finite difference form 
for computer solution. 

The simulations of various feed system components have been developed. 
These include feed tanks, valves, surge tanks, spring loaded accumulator sj 
other compliant devices such as bellows and screen channels, pumps, and 
liquid injectors. The component simulations were designed to use the 
characteristic wave equations together with the specified wave modification 
associated with each component to determine the pressures and flows in the 
system. 

A general computer program was written which uses the technique described 
above. Emphasis was placed on overall computational flexibility; the major 
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result is a program in which system changes can be made by data deck 
modification alone. 


TECHNICAL SUMMARY 


The characteristic equations governing transient flow of a slightly com- 
pressible liquid in a linearly elastic conduit are derived in Reference B-1*, 
and listed here for convenience. (In order to simplify the following dis- 
cussion, structural feedback effects and acceleration effects are neglected 
here. ) Along a left running plus characteristic, defined as in the nominal 
direction of flow or as moving downstream in a piping system, the equations 
are: 


dH a dV 

dT g dT 


0 


dX 

dT ' ^ 


( 1 ) 


For a right running minus characteristic, defined as moving against the 
nominal flow or upstream, the equations are: 

dH a dV • 

dT - g dT ^ (2) 

dX 
dT ~ 

It is more convenient to handle flow in terms of volumetric flowrate Q, rather 
than velocity V. This has no effect on the characteristic directions, but it 
does modify the derivatives along the characteristic. In terms of Q, the 
equations along the plus and minus characteristic respectively are; 


dH a dQ 

dT gA dT " 

dH a dQ 

dT " gA dT “ 


(3) 


B-1. W. H. Robinson. Liquid Propulsion Feed System Dynamic Analysis 
Program. Volume I, Technical Description Program H67Z. MDAC Report 
MDC G0451A, February 1972. 
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In order to get a physical picture of what Equation 3 represents, a simple 
pipe junction is illustrated in Figure B-1. In finite difference form, 
Equation 3 can be written for a plus characteristic as: 

AH + AQ = 0 (4) 

The speed of sound, a, is a function of liquid modulus K, considered a 
constant in lieu of temperature fluctuations in the line, and the wall 
thickness and elastic modulus of a particular pipe section. This is a con- 
stant for a single section of pipe, but may vary from section to section 
generally as the thickness and diameter of the pipe vary. The pipe cross- 
sectional area A may vary from section to section, but remains constant 
for any particular section. Consequently the group a/gA is a constant for 
any section of pipe, but is variable from segment to segment. If the 
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properties peculiar to waves moving with the flow are denoted by the 
subscript X, and those against the flow by Y, then: 


SX 

SY 



( 5 ) 


Equation 4 expanded in finite difference form, and utilizing Equation 5 is: 
+ SX (Q^ - Qj) = 0 

( 6 ) 

H 2 + SX • + SX • Qj = CX 


Equation 6 is the basic characteristic equation for a wave traveling in the 
direction of flow. In like manner it can be shown that the characteristic 
equation moving against the flow is: 

H 2 + SY ' = H 3 + SY • Q 3 = CY (7) 

The constants CX and CY described by Equations 6 and 7 show that the waves 
move unmodified through a section of pipe. This result is the direct con- 
sequence of the assumptions peculiar to rocket engine systems made in the 
derivation. This is not to say that the waves move unchanged through the 
entire system; they do not. The waves are modified at discrete junction? 
in the system, but pass from junction to junction unmodified. 

Equations 6 and 7 show no explicit time dependence. This again is the 
consequence of the simplifications made in the derivations. Time dependency 
is implicit, however, in the selection of the correct heads and flows to 
compute the constants CX and CY. Figure B -1 also shows the time distance 
diagram of the pipe system. It is assumed that the flow velocity is much less 
than the speed of sound. Therefore, the slope of the characteristic lines 
are +a for the plus characteristic, and -a for the upstream moving character- 
istic. The direction of flow is taken from left to right. In the diagram the 
pipe dimensions in Sections 1 and 2 are identical, although this is not a 
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prerequisite for solution. Thus both slopes are identical, but opposite in 
sign. These characteristics signify that wave perturbations in the system 
propagate with the liquid speed of sound. 

In the illustration, the lengths of the pipe system are such that it requires 
1. 2 computing intervals for the wave to travel from Station 1 to Station 2, 
and 2. 5 time intervals from Station 3 to Station 2. Consequently, the head 
and flow conditions at Station 2 at time t are a function of conditions at 
Station 1 at time t-1. 2At, and Station 3 at time t-2. 5At. The characteristic 
Equations 6 and 7 are computed as: 


CX = (t - 1 . 2At) + 

SX - Qj 

(t - 

1. 2At) 

CY = (t - 2. 5At) + 

SY-Q^ 

(t - 

2, 5At) 


The constants CX and CY are now specified and from Equations 6 and 7, the 
following relations also hold: 


CX = H 2 (t) + SX-Q^ (t) 
CY = (t) •+ SY- (t) 


These now are two equations for the two unknowns, and Q^, at time t. 
Conditions then at Station 2 are uniquely specified at time t as: 


©2 (t) 


( CX - CY l 
(SX - SY) 


(t) = CX - SX . (t) 


The calculation proceeds in the same manner for each junction in the 
pipe system, merely repeated once each time interval for each component 
contained in the system. The result is a lattice network of discrete points 
in time and space for which all the flow properties are known. 
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The concept of the characteristic carrying information from junction to 
junction, independent of all segments of the system other than those 
immediately adjacent to the junction under computation suggest the 
possibility of a modular or building block approach for the computational 
solution. In such an approach pieces can be added or removed from the 
system at will, without disturbing components not adjacent to the point 
of addition. 

To illustrate the technique, the schematic of the pipe shown in Figure B-1 
is reproduced in Figure B-2 along with a schematic showing the addition of 
a branch line at Station 2, For this case a third characteristic is needed; 
the one which runs from Station 4 to Station 2. If, for illustrative purposes, 
it requires 3. 1 computing intervals for a wave to reach Station 2 from 
Station 4, the characteristic equation for this section is: 

CZ = (t - 3. lAt) + SY^2* % - 3. lAt) 


I 

I 

f 

1 


I 

4 

I 


I 

I 

t 

I 



Simple Pipe Syitem 



Figure B-2. Simple Pipe Schematic 
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Where SY ^2 is the group -a/gA for the pipe section from Station 4 to 2. 

It is completely irrelevant, for this calculation, what the pipe system 
configuration is beyond Station 1, 3, and 4, as these regions do not con- 
tribute to the present calculation. They will of course have made significant 
contributions to the conditions existing at these stations, but these con- 
tingencies are assumed to be known prior to the initiation of the computation. 

The conditions at Station 2 for the tee has a system of four unknowns: H 2 
and the flows entering or leaving the junction at each of the pipe segments. 

A system of four equations is available to afford a unique, closed solution. 
For this discussion, the flows entering the junction are denoted by a double 
subscript, the first indicates the station from which the flow originated, and 
the second the station at which the flow is entering. Further, the flows are 
assumed to move from left to right, and from Station 2 towards Station 4. 
Using this notation, the system of equations describing the conditions at 
Station 2 is the three characteristic equations: 


CX = H2(t) + SX • 

CY = H2(t) + SY • Q23(t) 
CZ = H2(t) + SY^2 ■ Qz 4^^^ 
Plus from continuity 


12 


" ^23 


'24 
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Simultaneous solution of these four equations yields: 


■ CX 

CY 

CZ 1 

SX " 

SY " 

SY42 

■ 1 - 

1 - 

1 1 

SX 

SY 

■ 

fSI 

V* 

>< 

to 


^12 = (CX - H^) / SX 
Q23 = (CY - H^) / SY 

Q24 = ' SY42 

It is important to note the ease, from a computational point of view, with 
which the system has been radically changed. As far as Stations 1 and 3 
are concerned, the system change from the simple pipe at Station 2 to a tee 
connection on which a great number of additional components may be 
connected involves only the addition of a single flow variable, 023* 
other parameter designation for these stations remain the same. Only the 
sequence at Station 2 was changed to handle the tee connection. 

The convenience with which modules may be added to or subtracted from the 
system allows a variety of feed configurations to be analyzed by a single 
computational scheme. For example, engines may be added to or deleted 
from a system by converting simple junctions to tees. Lines could be 
shortened, surge tanks or spring accumulators added, all with a minimum 
of inconvenience. 

A computer program, H672, has been written utilizing the modular building 
block approach described above. The program can be conveniently sub- 
divided into four segments for description purposes; 

A. Main program — performs data initialization tasks and supervises 
the computational sequences. 
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B. Line routines — computes pressure and flow parameters in pipe 
segments by the characteristic method. The routines for a 
simple junction and tee were described previously and are 
representative of the routines grouped in this segment. 

C. Component routines — these routines simulate various components, 
such as valves, tanks, pumps, injectors, and combustors, 
contained within the system or at its boundaries, 

D. Auxiliary routines — these routines assist the line and component 
routines by handling such routine tasks as information retrieval, 
averaging, curve reading, and iteration selection. 

A few additional words concerning computational features of the character- 
istic method are appropriate. The selection of length increments is critical 
to the stability of the solution. In Figure B-1, the lengths were selected 
such that a wave required a minimum of one computing interval in traversing 
the distance between stations. This procedure must be rigidly adhered to 
if the stability of the solution is to be insured. 

The reason lies in the method of computation. For each increment in time, 
the flow properties are computed at each predetermined discrete junction in 
the system. Thus, in the example, when computing the properties at Station 2, 
we are assured of values at Stations 1 and 3 at times no later than t-At. 
Generally, within the system, this criteria can be assured; however, at the 
boundaries it may be impossible. Valves may be located less than one 
computing interval from an injector, and the only recourse is to shorten the 
computing interval, which of necessity lengthens the computational time to 
cover a fixed time interval. The selection of time interval of computation 
is then a general function of the system pipe segments. 

In calculating the properties at discrete points in the system, iterations 
are reduced to small closed loops within a single component routine. This 
reduces computational time and increases the stability of the solution. A 
large number of check cases have been run, and stability was never a 
problem even when large surges due to engine ignition spikes were 
simulated. 
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APPENDIX C 

SR 26 SCREEN DEVICE MODEL DEVELOPMENT 

The first step in the development of the screen-channel model is the deter- 
mination of the equation for acoustic velocity in the screen device (including 
properties of the liquid, gas bubbles, channel, or supporting structure and 
screen). It is assumed that initially the screen "pipe" is full of liquid. The 
equation for acoustic velocity in a slightly compressible liquid in an elastic 
pipe is well-known (Reference C-1). However, complications arise when the 
pipe is a nonuniform or composite structure, when the properties of the 
screen (hence of the pipe structure) are orthotropic, and when gas bubbles 
exist in the liquid. There are two kinds of sc reen- channel configurations 
which are amenable to analysis: the first is a uniformly distributed screen- 
channel structure which can be characterized as a pipe having properties 
based on the proper combination of screen and channel properties. (The 
second is the special case of screen only. ) 

Fine-mesh twilled-weave screens are basically orthotropic in nature, i. e. , 
they have different properties in the right angle directions of the weaving 
axes. The properties of importance are the screen thickness, e^, elastic 
modulus in the axial (Eg^) and circumferential (Eg 2 ) directions, and Poisson's 
ratio jJLgj^ (axially). Poisson's ratio in the circumferential direction can be 
derived from by: 






( 1 ) 


for orthotropic structures (Reference C-2). The channel support structure 

properties, E » M- » ® fraction of the channel area consisting of 

P P P 

screen, f, are also necessary. It should be noted that the screen fraction, f, 
may not be the actual fraction but, depending on the method of screen- 
channel fabrication, may be an "effective" screen fraction. 


C-1. J. Parmakian. Waterhammer Analysis. Dover Publications, New York, 
1963. 

C-2. M. H. Schneider and J. T. Hofeditz. Buckling of Fiberglass Cylinders 
Under External Pressure. ASME Paper No. 64-WA/UNT-12, 

December 1964. 
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Definition of the acoustic velocity in a pipe results from derivation of the 
continuity equation for one-dimension, unsteady flow and accounts for the 
elasticity of the fluid and the pipe. The pipe is shown in Figure C-1 and 
is assumed to be a uniformly distributed orthotropic screen-pipe structure 
with composite properties based on the fraction of screen area to total area: 


f = 


( 2 ) 


and screen and pipe properties (see Figure C-1) as follows (Reference C-2): 


The equivalent thickness, e, is; 


e = e ( 1 - f ) + f e 
P s 


(3) 


The longitudinal elastic modulus is; 


' = E (1 - f) + fE 

1 P £ 


(4) 


and circumferential elastic modulus is; 


E, = E ( 1 - f) + f E 
2 p s. 


(5) 


Similarly, the longitudinal Poisson's Ratio is; 


u (1 - f) E^ + f E 


P s 


1 “1 


E ( 1 - f ) + f E 
P s 


( 6 ) 


and circumferential Poisson's ratio is: 


^^2 ^ E ( 1 - f ) + f E 

P s- 


(7) 
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It will be noted that if the screen is not constrained by or attached to the pipe 
f = 1, and if the screen is completely constrained by the pipe, f = 0. 

Referring to Figure C-1, the deformation of the element of the composite 
pipe shell produced by a change in the longitudinal and circumferential 
stresses is 


AR = 


^A(r-> 


/A(T2 Ao-A 
\^2 " ^1 ) 


^ R 


1^2 


Act, ^iJA(r^^ 


and the change in the axial length of the element is 


6x = dx 


P2_^\ 
\^1 ■ ^2 / 


( 8 ) 


( 9 ) 


where AcTj^ and A(T 2 represent the change in the longitudinal and circumferen- 
tial stresses, respectively, produced by a change in pressure. The volume 
enclosed by the newly stressed element is tt (R + AR) (6x + dx) and the change 
in length compatible with the change in volume is 

2 2 
TT (R -I- AR) (6x + dx) - IT R dx 

ttR^ 


After neglecting terms which are very small when compared with those 
retained, the total change in length of the element is found to be 

6x+^^^dx (11) 

XV 

The change in the longitudinal stress is dependent upon the ability of the pipe 
to move in a longitudinal direction. Seyeral typical cases are considered in 
which the movement of the pipe is restricted to varying extents 
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PIPE PROPERTIES ELASTIC MODULUS, Ep 

POISSON'S RATIO, /ip 
WALL THICKNESS, 6p 

SCREEN PROPERTIES ELASTIC MODULUS, AXIAL DIRECTION, Eg^ 
ELASTIC MODULUS, CIRCUMFERENTIAL, Es, 

POISSON S RATIO, AXIAL, jUSi 

SCREEN THICKNESS, eg 

Figure C-1. Orthotropic Structure Nomenclature 


■A. For a pipe anchored at the upper end, free to move in a longitudinal 
direction throughout its lengths, and without expansion joints, the 
longitudinal and circumferential stresses produced by a pressure 
change of p dH are 


Aa 


1 


p dH D 
4e 


A^2 = 


p dH D 
2e 


( 12 ) 


from which the change in length from (11) is 


dx 


/p dH D 

1^4e Ej 


- I"2 


p dH D 
2eE2 



p dH D 
ITE^ 


|i^ p dH D\^ 


p dH D / 

^ -/>l) , p dH D A 


e E^ V 

4 / ^ e V 

" 2 / 


dx 


(13) 
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B. For a pipe which is anchored against longitudinal movement through- 
out its length. 


p dH D 

- l"2 2^ 


A.2 = ^ 


( 14 ) 


from which the change in length is 


j / p dH D p dH D\ , ^ , /p dH D p dH D\ 

'^’‘(>‘2TTe7 


p dH D 1^2 

\ . p dH D / 


e E^ \2 P'2 

•V e E, ( 

'2 j 


dx 


( 15 ) 


C. For a pipe with expansion joints between anchors (or free at both 
ends): 


^ CTj = 0 


Acr- 


p dH D 


2 2e 

from which the change in length is 


( 16 ) 


dx 


/ p dH D\ 

(-^2-2FE^) 


p dH D \ P dH D _ p_dH 

^^2 2e E, ) ° 2e E, " e E 




dx 


( 17 ) 


Therefore, the total change of length produced by the pressure change p dH 
is : 


dx C, + dx C, 

e Ej^ 1 ® ^ 


( 18 ) 
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where for case 


1 - 2 p. ' 

A. ^ 

p.~ 

S' = — -|ii|X2 


C. Cj = 0 



The change in volume of the original element of fluid, dx in length, because 
of the elasticity of liquid under the action of a pressure change p dH, is: 


p ttR^ 
K 


dH dx 


(19) 


and the corresponding change in length of the element of liquid is 


P TT 

K TT R^ 


dH dx 


^dH dx 


( 20 ) 


The total change in length of the element of liquid caused by a pressure 
change p dH, when both the compressibility of the liquid and the deformation 
of the pipe are considered, reduces to 


P 


dH dx I ^ 


D ° ^2 \ 

+ TE7 + T^j 


(21) 


Since H is a function of x and t, and dx/dt = V, it follows that 

The total change in length of the element produced by a pressure change is 
then 


1 ,°Cl ,PC2 

K e E, e 


\0 


8H , ,, 3H 

;V 


) 


dx dt 


(23) 
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The total change in length is also 


3x 


dx dt 


( 24 ) 


so that 


av 

9x 



+ 


e I \dt 8x 


) 


(25) 


or 


+ V — = 8V 

8t 8x g 8x 


where the acoustic velocity a is 


a 



D C. 


e E, 



(26) 


The extension of the basic pipe analysis to liquids which include gas bubbles 
is given in Reference C-3. The experimental data of Reference C-3 is best 
correlated by assuming adiabatic bubble expansion and contraction. The' 
final equation for the acoustic velocity in a liquid-gas mixture in a composite 
orthotropic pipe is: 


X 


(1 - x)' 




+ X { 1 - x) 






_!s_\ 

K g) 


Py P 
(1 + 

g g 


\ ® ^1 ® ^2 


T 1/2 


(27) 


C-3. B. R. Hanks and D. G. Stephens. Helium Injection to Reduce Resonant 
Frequencies in Propellant Lines. Journal of Spacecraft, Volume 6, No. 11, 
October 1969, pp 1202-1204. 
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where x is the volume fraction of gas; C and C- are the acoustic velocities 

g ■L' 

in the gas and liquid, respectively; p and p. are the corresponding densities; 

S i-* 

is the absolute pressure in the pipe; y is the ratio of specific heats for the 
gas; g is the gravitational constant; K is the compressibility of the liquid; 


D and e are the pipe inside diameter and equivalent thickness; and E 2 are 
the orthotropic elastic moduli; and and are constants as given in 
Equation 18 depending on the form of pipe constraint. 


The complete screen device will be divided into a series of sections, one of 
which is shown in Figure C-2. The screen area is defined by the screen frac- 
tion, and the length of the section 


A = TT • D • L • f 
s 


(28) 


The upstream and downstream conditions of head and flow are known from the 
previous timestep so that 


H4 



CR21 


Figure C-2. Screen Device Model 
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CX = HI + SX • Q1 

CY = H3 + SY • Q5 (29) 


where 


SX 

SY 



P 


P 


Some of the Q4 flow, flows through the screen (Q6) because of the head 
difference: 


H2 - H4 


A'Ob 



s 


(30) 


where A' and B' are the experimentally determined liquid flow-loss 
coefficients for the screen, and H4 is the tank pressure. 


The rest of the flow, Q4 - Q6, expands the pipe slightly such that the change 
of volume is : 


AT (Q4 - Q6) = L 2ttR AR 


(31) 


where the appropriate time increment is that required for the wave to traverse 
the section length: 


AT 



(32) 


and where AR from Equation 9 is: 


AR 



^1 / 


(33) 
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For the pipe free to move at the edges of the section. 


Ao-j^ = 0 


so that combining Equation 31, 32, and 33: 


A (Tp 

Q4 - Q6 = 2 A a'-^ 
p E^ 

(34) 

Taking a pressure (force) balance on the pipe, considering the 
in the pipe: 

stress change 

Pl- H2 • - A^) = PlH4(A^ ‘ ^ 

(35) 

where A is the open area of the screen and P = A /A , thus 
o o o s’ 


AcTp e 

®-H4 = ,D(1.P ), 

E 0 

(36) 

Solving for Atr^ and substituting (34) into (36) gives 


e Ep (Q4 - Q6) 

H2 - H4 - -r-7 r~^ ^ 

2A a’pT'n'D(l-P)f 
p o 

(37) 

From Figure 3 and the characteristic identities: 


Q2 = (CX - H2)/SX 

- 

Q3 = (CY - H2)/SY 


Q4 = Q2 - Q3 

(38) 


Equation 38 is substituted into Equation 3 7 to eliminate Q4 and the result is 
combined with Equation 30 to eliminate H2 and give a quadratic expression in 
only Q6 and H4. Solving this quadratic explicitly for Q6 (flow through the 
screen) gives 
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where 

ZZ 

a 

b 

c 


1 + 

SX SY 
A'/A 

s 

B'/A ^ 
s 

2e E-/a' (1 - P ) f 

2 L. o 


From Q6, H2 is computed from Equation 30 and the flows from Equation 38. 
Q3 and H2 then become the upstream values used to compute the new condi- 
tions in the next section. 


With gas surrounding the screen, a negative pressure pulse in excess of the 
screen bubble point will result in gas ingestion with the inflow Q6' repre- 
sented by a new form of equation (30): 


H4 - H2 


A" Q6' 


, B" (Q6')^ 

A 2 
s 


(40) 


where A" and B" are the flow-loss coefficients corrected to gas properties. 
Q6' is found from Equation 39 with the appropriate values of a and b. 

In order not to permanently deform the screen, the equivalent stress in the 
screen, Ao-2, defined by Equation 36 for the H2 solution, must be below some 
critical stress, e. g. , the proportional limit of the screen. Equation 36 is 
used to check this condition for the screen device. Another criterion which 
must be checked when gas surrounds the screen device is whether liquid 
expelled through the screen by a pressure pulse leaves the screen. 


181 



The velocity of the liquid through the screen pore is: 


V 


Q6 


A P 
s o 


(41) 


and the velocity of the liquid globule (if r^^ > r^) is approximately: 


= V 

L s 


/ 

fe) 


(42) 


where r^ is the globule radius and r is half the screen pore diameter, other- 
i-f s 

wise V, equals V . 

L 8 

The globule volume is: 


QL = yir r^ 


(43) 


and for n pores where: 


n = 


A P 
s o 

2 

TT r_ 


(44) 


then the globule radius is: 


QLSUM 3 2 

A P ■ 4 
s o 


-I 1/3 


(45) 


The quantity QLSUM equals QL * n and also equals ZQ6 * AT and is updated 
with each time step. 

The liquid globule kinetic energy is 


K. E. 


4 3 

^L 



(46) 
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and the surface tension energy is approximately 


S. E. 


3 

O’ IT r 

s 



( 47 ) 


Assuming negligible potential energy, the criterion for liquid breakaway is 
for the globule kinetic energy to exceed the surface tension energy. Com- 
bining Equations 4 6 and 47 and solving for the critical velocity, 
which results in liquid globule breakaway gives: 


V 




3_ 

4 



1/2 


(48) 


If V from Equation 42 is larger than from Equation 48 the liquid 

globule will leave the screen. In addition, if the quantity of liquid outflow . 
is sufficient to wet the entire screen surface, it was assumed that screen 
pore surface tension forces no longer exist, so that the liquid also leaves 
the screen, The quantity of liquid (QLSUM) held by the screen is kept cur- 
rent by correction (if any) in each time step. If QLSUM is positive, then 
inflow during a negative pressure pulse is liquid, until QLSUM goes to zero. 
This may retard gas ingestion for several time steps. 

The model for gas ingestion and quality assumes that the quality in the device 
is uniform (i. e. , the gas bubbles are not concentrated near the screen, but 
are distributed uniformly throughout the liquid in the screen). Throughout the 
H672 program, the liquid density is used to compute head, so that for gas 
inflow, the volumetric flow rate is corrected by the gas-liquid density ratio. 
The quality in the device is equal to the gas volume divided by the liquid 
volume (screen device volume minus gas volume). If outflow from the screen 
device through the downstream valve is occurring, the gas volume is dimin- 
ished by the amount of entrained gas that outflows (conversely, if inflow 
occurs, and the downstream quality is not zero, the gas volume is similarly 
increased). If the quality in the device is not zero, and a positive pressure 
pulse in excess of the bubble point occurs, it is assumed that the gas volume 
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is diminished by the outflow times the current quality. The gas volume 
(QSUM) is kept current by correction (if any) in each time step. The acoustic 
velocity a', from Equation 27 is updated in each time step by the quality 
from the preceding time step. 
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APPENDIX D 


000035 

OOOODr 

90000^ 


OOOOOB 


900025 

900036 

900063 

900067 

90007? 

900075 

OOOlCl 

900104 

OOOlU 

gOOlZi 

®00i2(t 


900i3i 

900134 

900134 

000134 


900134 


subroutine S 2 RPRP (SljRXiI) 

c 

C THIS SUBROUTINE" IS USED TO PREPARE THE SUKX DATA BLCCkS 
C rOR COMPUTATIONS IN THE SR26 SCREEN DEVICE SUBROUTINE 

C 

DlMEflSlON SURXl20i4 ) ,h1(5) 

C0MH0N/An/2.CNS7<5 > iRH0(2) 

COMMON/TURai'l/XLiTHPP.WTTfXKiVO 

C 

C 

READ (t>ilOO) Hl» (SURX(k|.I)iKsiil 6 >,N 

write <6,1015 <Hl(K>",Kii>l,5).l,<SURX(K‘, I).Ksi,i 6 ),H 

SURX(3,l) a sURX(3,I)*eNST(l) 

SURX(5,I) 8 3URX(5,l5*CNST(.l).*CNST{l>- 
SURX<6, 1 5 ..',a SURX < 6 , 1 5 *12 i . ; ■ — 

SURX(H, ifV.;a-SURX(i'l,4 )’/l2. 

SURX<iii»l) a- EURX(12, 1 J/CNST(i) 

SURXdPii) a SURXU5, 1 )/l2. . 

suRX(i’/,n 8 sSurx( 3; i)^suRX(2, n 
SURX(l«, I ) a.:SQRT(5URX<l, I ) *CkST ( 4 ) »Si;RX t 2 , I 5 > 

SURX<i 9 ,I) a RH'0 <N)/SURX(i, I )/SURX{3* I )/CKET(4) 

SURX(20,n a SURX(17, n/XK/CNsT(4> 

SURX(3il) a sURX<3, I )/RH0(N) 

RETURN 

lOfl fCRMAT (544,6FlC.3,/,20X,6ri0.7,/,2oX‘,4F-if ,7,113) 

inl format <lHc,//,25X»4aHDESlGN CATA FCR SCREEN DEVICE LOCATED aTiIXi 
15A4i//iS4X,9HSEQUENC£ , I4«///,iX,8HGAKMa a , Fg . 3 . 5X , BnRTGAS a ,FH 
2.3 i 5X,7HPGAS = • F6 , 3 • 5X# 6KBPS a , F6. 3'. 5X , EE a , F n , 3 , SX , 7HSTEM * 
3 |Fll,3 »//, iVi 6 HAS 1 a ,F9,6,5X, 6 HES 1 a , T9 , 6 • 5 X , 6 HAS 2 a »F9.ft,5Xi6 
4IIDPS a ,F9,7,5Xi5HCS a , F 6 .2 , 5X, 5HDC a , F 9 , 7 , // , 1 X, 5 HPO a ,F6,4,5.X 
5.4HF , .,Fs.6,5X,5HLS a , F 6 . 2 , gX , 5HPV- c ,F6,4,5y,4HN a ,11) 

END 
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000006 


-0000.06 


00CD06 

000006 

000006 

000006 


00 000 6' 
000007 
OOOOlO 
000020 
000022 
000024 
000032 
000034 
000034 
000036 


000041 

000047 

000065 

OOOlQl 

OOOIC.7 

000110 

000i26 

000143 

000151 

000153 

.000157 

000163 

000165 

000167 

000170 

000171 

000172 

000173 


APPENDIX E 


SUP R PUT I N E S P 2 6., <,s y Jj Jl 


This. SUBROUTtNE COMPUTES THE PARAPETEPS 
FOR A SCREEN EPVICE THE P6ED ClflCUlT 


DIMENSICK SURX<20.4) I J(20«50)i AP(20) .AS(2d>i>!( 20) rCLSUH(20)i 

1 OA(20).Q6(20>. AA Si(20).AAS9(90).6BSi (20)iCCSl (20>. 

2 yU<5(i)itHpP(3C>AWTT(60),)iK(2>iyo(2>iZZ(EO)»RCCT4(2C)i 

3 WDOT6(20).STE{20),VSUM(;2S),QSUK(20>,RU:2C),VL(20>,VCt:2C) 

COMMON/APV? I CNST ( 5 PHP( 21 > NPL (.4 . N P M M > . D ( 4 ) . T I PE » KFT t hPA6'T . 

lSX(64).SY(64)iAREA(4:ltBL0CK(60)iFV(2l.F(50>«Kl(3Cl»NCAT(60)* 
2^<4 ()i6q)iGUoi6C )*FS(6c)iHOOT(6C)i>HC(6o),GHC{6o’ >^(60liTt(5C)i 
.- 3 ° at { 5 c > , aCL ( 4(i'.5.->.XOOT( 40.5c).Sl6(?> 

common/seeker/mjblk 

COMmON/LINEB/HDN(4) 

COMMON/TURRlN/.yyiTHPPiWTTA_XKilO. 

EQUIYALEKCE <Ki(l>,Ji),{Kl(2)# J2 >i (KK3). J3)i {Xi<4>i J4)i <K i(51i is) 
1* (Kl<6).J6»i <Ki<7)iJ7)i (Kl(8>, JB)i <Ki(9>.j9)» (Hi(lo)iWlO)i CKHll)i 

3U16|i (Ki,<i7)! JI7 ) I ( Ki <18 ) « Jl8) i (><1 ( 1« > • J19 > i < Ki ( 2 0 > * U20 1 ^ 

4 1 <Ki< 21 )*vl2l’« <Xl(22l*J22>»<'<l<23>i J23)» CK1(24 >im24),(Ri(2?)iU2&), 
5<Ki<26) . J261.I («1<271* ‘ ‘ .»-> 

••• control block conversion loop 


12 

1,1 

13 


21 


22 


NjBLK '« 'N 

DO 5c l«iil8 

Kl(I)8j(I»N) 

N2bnDAT( j2» 

N4,NPAT< J4) 

GO TO <lliI?)|Jie 

S = -1 , 

Go TO 13 

S = 1 .. . 

IF (TIME»TTCJ9) ) 2C.30.30 

•••SCREEN Piyicfc Parameter s prior to wave arrival ••» 

GO TO (21,?2),J17 

H < KPT, J3 ) »HHD ( J2 ) *ACL (KPT|jli_6 ) •DA T ( J16 > .ACL <K PT,N2)«P A T(N2 > 
G(KPT, J6)»!5HD( J5 )*ArEA( ji4}*(XU0T(KPt. Jl6)..XCOT{Kf>T,Ngn 
Q<ypT, J7)«a<KPT, J6) 

GO TO 23 


23 

24 


23 


wu 1 1; _ 

H ( KPTVjs ) “HHD < J4 > * ACL { KPTiTiT>'*ijaT( Ji6 ) •ACl'CK'PTrN 4T»DAtTN4T 
0O<PT, J7l«S»0HD< J8 >* area (J 14)*<XCCT (KPT, Jie>»XCOTiHPT|N4)) 

Q(KPT,J61«0(KPT»J7) ; 

IF (TlHEJ 24,24,25 
HHr(J3}»H(KPT, J3) 

OHnc J6 )sO.<KPT, J6> __ 

3HnCJ7)sGHn( J6’> 

X(Ji) B 5, 

VSllM(Jl) 4 0,C 

OSUM(ji) 4 0,0 

OLSUM(Ji) « 0,5 

KPPac 

GO TO 4n - . - 
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000173 

000i77 

000203 

000206 

000207 

000211 

000214 

000216 

OOO227 

000233 

000236 

000241 

000256 

0002-61 

000263 

000270 

000314 


000337 


000345 


000354 


000366 


000431 

000432 

000444 

00Q455 

000467 

000470 

000503 

000513 

000523 

000532 


C •••SCREEN DEVlCfe PARAMETERS ApTbR INITIAL LAVE ARRIVAL ••• 

C - 

31 CALL PAST(TlMc'.T(Jil».2) 

CALL PASTITIME.TI Jj0)ii } 

Ap(Jl) a LTTIJiZ) 

SXAV s SX/J12) 

SVAV a SV|J133 

HVAP*PV(Ji 5)^CNST(, )/RH0<Wiii) 

APOLD 8 AP(Ji) 

AS(Ji) 8 SURX(i5»Jj )*CNST{3)*SUKy ( 11 , Ji >*SLRXJ14»M1 ) 

AASi(Ji) ? SURX(7iJi>/A,S(Wl) 

AAr2(ji) a S'-'RX(9.Jl>/As<Jl> 

BBPl(Jl) , SURX/6.Ji,/AS(J 1)/AS( jl, 

CCSi(ji) a SURXis, Ji5 82./’8H0(Jl5>/CNST I3)/CM§T<3)/ai"§LRXU3i JDJ 
1 /AP< Ji)/S'-'RV(i4iJi) 

SXAV 8 SXAVaAPI Jl>/A^0LD 

SYAV a SVAV*AP« Ji)/APOLD 

Z2<Jl> s i7 ♦ CCSi(Ji)/SXAV- CCSi {.D/SYAV 

CY ■ RCCNTHi J^ ,2USYAV^S^RC0N(Q» je ,2UACL<KpT| »16}*DAT<wl6) 

1 -RCCMACL.N4i2)*DA7(N4) 

CX 8 RcCN(HiJ 2.1>*SXAV*RC0N(«, J5,j >*ACL{KPT,Li4>^pAT(*ifc> 

1 •RcCMACL'.N2.i>*04T<N2) 

C 

c 

c CHECK IE Screen segment surrounoei by liquid or las 
c ■ . 

IF (SURX(16.J1» .GT: 00 to 14C 

c 

C SCrEEn SURRC'JnCEO 3y gas 

c 

c check if Screen pressure asove or beldl ullage pressure 
c 

131 IF {H(KPAST,J3J-SURX(3.JI)) I50*l58,l40 

c 

C CHECK IF CElTAP BElO'« POUT 

c 

151 U (H(KPASTiJ 3)-SURX<3, Jl>*SURX(4,„in l60'»l70il7C 

c 

C LlPUlD inflow or outflow 

c 

14'j 06(J1) a 2,*(SURX<3,Jl>-(SgRX{3.U})tCCSliJl)*5X('5?CAV9CCSlj4/l>»CY/ 

1 SYAV)/ZZUi))/(«AASl( JU.CCSl(Jl) /ZZ(>l)«SQHt UAASl(ji) 

2 •CCSi( Jl)/ZZ<J?n^*2*ABS{4,*BBSl{Wl>*«SLRM3i J1>*(SLRX 

3 13* Ji>*CCSi( Ji>*CX/SXAV-CCSi( J i)*CY/SYaV jAZZ(Jj))) )) 

IF (ri6(Jl);LY, 0.0> GO TO 142 

R<KPT,J3> • SURX<3. Jl)*AA5lC Jl)*06(Ul>* 

1 RBSl(Jl)^Q<S{dl>*0§(Ul> 

IF <H(KPT,J3)».SURX<3iJi)»SgRX(4;jjj) j 4 3 i 4 5 V 1 4 6 
146 QSUMIji) 4 QSLiMI Ji)*Z/SLRX(i5i Ji>/AR£A<i„14>»X{Ul} ) 

145 GO TO 141 

142 HIKPT.JSI 8 SUPX(3i Jl>*AASi<ji>*C6(Jl)« 

1 0BSK Ji)*Qft(Ul>^Q6(Jl ) 

141 IF (H{KrT, J3) (LT.HVAR) H<KPTiU3) • HVAP 
I5(KPT,J7» 8 (CT-H(KPT, J3))>SYaV 
QCKpT.Ufil 8 <CX-H(KPT, J3>)/SXAV 
04(J1) a 0 {KPT‘iU6>.0«KPT»»I7) 

c 

c check liquid leaving screen 
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000340 

000344 

000334 

000362 

000362 

000366 

000372 

0006^0 

000623 

000633 

000636 

000645 

000646 

000660 

0007CO 

000702 


000703 

000724 

000731 

000733 

000742 

000742 

000746 

000754 


000763 

000765 


001030 

001031 

001043 

001044 

001057 

001062 

001063 

001065 


00H34 

001135 


C 

CSURXfl6i Jl) ,GT;Oi05 70.148.,.. .. 

IF (W(KFT, J3l-SURX(3iJi5«SURX{4i Jl>> I47il47il5l 
15t OCSUH{Ji) a QLSUM(Ji)*(56(»(1)*Z'^<1.»X(Ji>) 

GO TO 

147 OLKUM(Ji) B f)LSUM( Jl)*06f J1)*Z 

152 IF <nL5U^'(Jl>,LE.0.o’ OLSUM(Jl> b g,o 

IF (OLSUR(JiI,r,T,SURX(io,Ji)/SOHT{SUPX<l3»Jin*Ait'(,l>)l)lSVR(*)l»»Oj 
’^U(Ji>b(3,/i6,*SURX<i0i J l)*SURX(iO, J i)/AS(Ji)/SUHX(i3i J l>* 

1 CLSUMUi ) >**.333 

IF (RL(J1);GT,SURX(15,J1)/2.) QO TC 143 
VU(Jl) e 06( Jl»/AS( J1)/SURX(13» Jl) 

IF (RLt Ji I ,UE, 0 i 3 l RL(Jl) “ C. 1 *SlRV (lO* JH/ 2 * 

GO TO 144 

143 Vt(Jl) B c6(Jl)/As{ Jl)/siJRX(l3,Ji )>4.*SURX<10|»;1J*SUKXU0| Jl>/ 

1 RL< J1)/RL( J1 ) 

144 Vc<jl) » S3hT(Si6<jl5).(§URX{10#Jl)**3)*3,/32,/RFj:(Jl5)«CNST(4)/ 

1 {RL(Ji)*»4n 

IF (VL( Ji ) ,LP,VC(Ji ) ) C.O TO 14« 

Ol,PuM(Jii s OjO 
C 

C check ScRfc'FN STRESS 
C 

145 STE(Ji) a (H<KPT, J3)«SURX(3i Jl) )*P(.0(J15 >i»CnST(3>*5UhX{Hi 

1 <l.*SllRX(i:^i Ji))*SURX<j4iJi) 

IF (ABS(STE( Ji) ) .LT,SURX(6*'Jl)> CC TO i49 
STE( Jj ) • STEUl>/SuRX(6iiJl) 

WRITE (6i 2§1> Jl.STE(Jl) 

2ni FORMAT(5^>3HWg;, I3H7H SCREEM PRpFCRT I CNAL tIHIT tXCfeEDtp » STRESS 
1 ' iFii.3iTX,17HTH£S HaX, STRESS) 

149 IF (OSUKNO.LF.O.i'I GO TO 43 
VSUM(Ji) f yS0M( J1)*3<KPT,J8)*Z 
IF <VSUK(Jf )-SDHX<f5, Ji)*6REA< jl4) ) 3o5»3o5i3ll 
C 

c Gas inflck 

16'; IF {C!LSlJf'( Jl) .LE.oVd) go TO 163 

' 56 (J.l> B 2 ;»<SiJHX( 3 , Ji).(SyRX( 3 ijj l + CCSlC JJ>*pX/SXAV»CCSlIx'lIi'CY/ 

1 Sva'M/ZZC Ji) )/(-AASi<vii).CC*i( Jl)/ZZ( Wl>BSQRTr(AASi( Ji> 

2 -ACCSlCJlI/ZZC j; >)**2*AB5(4,»BBSl(v’l>*(SURM3»,,)l)"(SUhX 

3 (3i JD*ccSi< Ji )*gx/sxAV-ccsi ( ji)ACY/syAV)/zzjji> )) >) 

IF (U6(w1 ) ;L7, 0.0> go to 164 

H<KpT,J 3 ) a SURX(3«J1>*AAS1<J1>*C6(J1>* 

1 GBS1< Ji )*3 a< J l>*W6< ) 

GO TC 161 

.161 H<KPT,J3) b SgPXCl, J1)*AAS1< J1)*C6 { J1 >- 
• 1 RBSK J.i >*Q6( J1)*Q6( Ji ) 

QLSuf-'(Ji) » OLSUM( Ji)*c,6<*JlI*Z 

GO TO 161 

163 OLSuf'(Jl) e 0.0 

Jl> 8 27 * f SURX< 3 , Ji)-(SURX{ 3 » Ji >*ccsi ( Ji )*CX/SX aVbCCSi(»< 1 >aCY/ 

X SVAV)/ZZ< Ji>)/(-AAS 2( Ji).CCSi( Jl)^ZZUl>"SORT<( AASg( Ji) . 

2 *CCSi( Jl)/ 2 Z< J?))** 2 *AriS( 4 ,*BBSlNl>*<SUFiM 3 idl)»isUf'X 

3 <3* Jl)*CCSi<0i>*CX/SXAV-CCSi{Ji)*CY/SYAV)/22Ul) )*SlihX<i7 

4 » J1)/RHO(Ji5)))> 

IF (Q6<Jj ).LT, 9,0 > GO TO 162 
H(KpT,Ja) B S0RX<3» J1 )*AAS2( jl)*C6( Jl)* 

1 RBSl( Ji)*Q 6 ( Jl)* 06 (Ji >*Sl:RX( 17 i Jl>/RFC( Ji 8 ) 
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001154 

001155 

001l7« 


001206 

001227 

001234 

001236 

001245 

001245 

001245 

O0i255 

001264 

001272 


001275 

001277 


001342 

001343 


001355 

001356 

001371 

001374 

001376 

001406 

001416 

001425 

001^33 

001434 

001436 

001437 

001456 

001^65 

001475 

00l5o5 


001514 

001522 

001530 

001534 

001536 

001541 

001542 


GO TO i6i 

162 H(KPt.J 3> » SLIRX<3*Ji)*AA52{Jl).*o,6.(JlO« 

1 9BSKJi >«06(Ji)*Q6(Jl)*SURX(17i Jl)/RHt{Jl5) 

ior If (H(KPT, J3 >.L'T,HVaP> H(KPT,J3> ■ HVAP 

c 

c check Screen stress 
c 

169 STE(Ji) ■ {H(KPT,j3)»SiJRX(3*JlJ>*Pj.,0(J15)*CNSTI3V*Y^'.RMlli»ll>T 
1 <1.*SURX<i3*Ji>)*SURX(i 4. Jj) 

If {ABS(STE(Ji>).LT.SURX(6iJin CC TC i69 
STE(Ji) ■ STE{Ji)/SUPX(6iJ1> 

Write (6,2o2v ji.steud 

2o? fOPMAT(5X,3HN0:, 13.47H SCREEN PRCFCRTIONAL LIMIT EXCEEDED " BTRbSS 

1 “ •Fll.3iiXil7HTlME? HaX, STRESS) 

166 CONTINUE 

0IKPT»J7) ? <CV*H<KPT, J3))ySYAV 
Q<KRTiJ6> b (CX-H(KPT,J3))/SXAV 
04(J1) 8 Q(KpTivI6),q(KPT*J7> 

60 TO 180 
C 

C NO INFLOW (RESISTED SY BUBBLE POINT) 

C 

17rj If {°LSUM< JD ,LEi0.'0I 50 TO 173 

06 (Ji) B 2T*<5URX(3. Jl)-«SgRX(3*JH*CCSl( Jl)*CX/SX4V>CCSK»ll>*CT/ 

1 SyAW)/Z7(Ji)>/<-AASi(Ji).CCSl(Jl)/Z2U'l)-60RT( UASil^^U 

2 ♦CC5l(Jl)/27{Ji))**2!»A5S(4,TBBSl(gi)*ISLRM3iJa>»{Sl!hx 

3 (-^iJl)*CCSi(Ji)*CX/SXAV.CCSi(Ji>*CY/SVAV)/2Z{Uj) )))> 

If <'S6(J1);LT, 0.0> 50 TO J74 

W(KPT,J3> » SURXtSiJD^AASltJDBfiiUD* 

1 BBSi( Ji )*(Ja(01)*06(Ui > 

GO TO i71 

17-1 H(KPT,J3) 5 SLRX(3,Jl)*AASl(Jl)*ce(Jl)» 

1 BB'6l( Ji )*QA(dl>*Q6<Ul > 

OLSUM(Ji) 8 aLSUH( Ji)4!56< JlJ*Z 

If (OLSuP(Ui),LE,o;o> 50 TO 173 

17i If <H(i<PT,J3),LT.HVAP) H(KPT,J3> » HVAP 
Q(KPTiJ 7) 3 (CT-H(XPT, J3) J/SYAV 
0(WPT,J6) * (CX-M(KPT, J3))/SXAV 
Q4(Ji) 8 Q(KPTiU6)*>Q(KPT»*I7> 

50 TO 100 

173 0LSUM( Jl> □ o,C 
06(Jl) 8 oT 

H(KPT, J3) a (SURX(3, Ji)*CCSl(Jl)a CCX/SXAV-CY/SVAV ) ) / ( 1 , ♦ CCSi (gi ) / 

1 SX*V-CCSi( JlJi^SYAV) 

IF (H(KPT, J3) .LT.HVaP) H(I3P7#J3) ? HVAP 
0(KPT,J7) a (CY«H(KpT, J3) IXSYAV 
0(KPT,J6) = (CX-H(KPT,J3>)/SXAV 
04(Jj) a Q(KPT'. J6)*Q(KPT. J7> 

C 

C quality computation 

V 

lel VSUMljl.) » VSUM{Ji)40{KpTf J8)*z 

If <VSUP(Jl )-SURX( is,Ji)*AREA<Jia)) 300i300|3iC 
3e.! IF (Q(KpT,ja>) 381,3/3,301 
3n IF (VSUP(Oi) ,LE.C.:> GO TO 303 
IF (OLSuYUl» .GT.oIqI '■’0 TO 305 
IF (OLSbMUl) ,Lr.oilo> 50 TO 305 
IF (06(J1);GT, 0.0) 50 TO 305 
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001545 OSllM( jl)B0SUMUn*<l.*-5(KPT» J8>*2/SUPX(15, jl)//iRfc/K > )«Q6( Jl>*2 

001561 GO tG 3fl2 

001562 3(l3 OSl!M(jn e QSUH ( Jl ) ♦ ( l , .0 { KPT , JH > *2/SURX U5 , J1 )/ARfc A U14 ) ) 

001375 GO TO 3o? 

001576 3n5 IP < QlSu»' ( > 11 ) , CT , j ; 3 ) GO TO ?02 

001601 If (06(J1>IGT| C.0> to 302 

001603 OSUM(Ji) » CGUMUi). 06( Jl)*2 

OOI606 GO TO 3a2 

0016 07 3fj3 X( Ji)BQSUMIJl)/(CNST(3)/4t*SiJRX(ll, Ji J*SURMlli J1 J*SUfiX{i5» Ji)(* 

1 OSUM(Jl)» 

001622 IF <X{ Ji) .LE.oToI X(Ji) • 0,0 - 

001625 IF (X(Ji),L'T, 1 Tc) GO TO 40 

001630 X(Ji) E 1,0 

001631 . write <6.2331 Jl 

0Oj 637 2Cj F0RMAT(5X,3HN'0.. l3.4>^'H -'SCR6E!!l DEVICE FULL CF GAS ? QUAUTT EXtEEES 

1 1 . 0 ) ^ , 

001637 GO TO 40 ' 

.004642 3i;i VSlJMt Jt B'.Cr KPj i ja)*Z 

001647 GO TO loY :• tl-'y-;- ■ . 

001650 3ii VSUM(jir« C<KPTrJ6)«Z:*-':. . 

001655 GO TO 305 

C **B PRESSUflEVF*liCwRA:TE,HQLO CONVERSION SEGKENT 

C ■ • 

001656 4 'i BuoCK< Ji*32> * X(Ji) 

001660 8l0CK<Ji* 36) * AP<Ji) 

001662 BLOCK(Ji*43) s 86 <Ji) 

001663 IF <H<KPT, J3) .LT.HVAP) H<KPT,J3> = HVAP 

001673 DELTAP b ( H< KPT. J3 ) -Sg^X <3i Jl ) >*R^'C ( Jl5 ) -'CNST ( 1) 

001703 B^oCk(J 1*44) f OEltAP 

001705 Bi. 0CK(J1*66> a OLSUM(Ji) 

001707 PS(J3) B H<KPT. J3)*RH0( J15)/CNST(J) 

001715 WD0T<J6)»Q<KPT'. J6)*RW0< Ji5) 

001722 WD0T(J7)b0(KPT', J7)*RW0(J15) 

001726 return 

001727 END 

•following Variables ecuivalenced but not refewencec 
J i9 

J2b 

J21 

J22 

JS3 

J24 

J25 

JS6 

JS7 

J26 

J29 


J30 
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Appendix F 

TECHNICAL DESCRIPTION 
CRYO-LINE PRESSURE PROGRAM, P4557 

When a cryogenic liquid flows into a warm transfer line, boiling of the liquid 
will occur until the line has been chilled to liquid temperature and liquid 
transfer can occur. If the liquid flow is initially rapid, an unstable situation 
will develop with surging flow rate and line pressure. Liquid will flow into 
the line due to the lag in boiling caused by a finite heat transfer rate; when 
boiling does occur, a relatively high velocity liquid flow will have been 
established. The vapor production rate will accelerate rapidly as boiling 
causes liquid breakup and entrainment, with an expanding heat transfer area, 
and will exceed the rate at which vapor can escape from the line, causing a 
line pressure surge which will slow and reverse the liquid flow. When the 
liquid has been forced back into the supply vessel and/or has boiled away, the 
pressure will recede and liquid flow will again enter the transfer line, repeat- 
ing the process. The intensity of the flow and pressure surge will diminish 
as the line is cooled; the maximxim pressure surge occurs on the first cycle. 

The actual flow process during cooldown is very complex; however, a simple 
model is adequate to analyze the first cycle, which will predict the peak pres- 
sure and the time scale of the cycle. The fluid flow model assumes a uniform 
slug of liquid up to the liquid-vapor interface; all vapor production is concen- 
trated at the interface with a discontinuity of density, velocity, and pressure 
at that point. A separate heat transfer model accounts for the expanded heat 
transfer area of the boiling liquid due to the expanded volume of the two-phase 
fluid. 

Other assumptions of the model are as follows: 

A. The vapor flow is adiabatic. 

B. The line is horizontal, straight, and constant cross-section area. 

C. The inlet valve effective area varies linearly with time from zero 
initially to the line area when fully open. 
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D. A step function initial temperature distribution for the line, with the 
initial section at liquid temperature and the downstream section at 
ambient temperature, is set to give approximately the same total 
available heat to chilldown. 

E. The boiling heat flux has a constant empirical value for each fluid. 

F. The wall temperature is constant. 

G. Quasi- steady- state Fanno flow is assumed for the vapor stream. 

The analysis describes the pressure distribution in the pipe as liquid flows in 
from the supply tank and gas flows out at the delivery end. The pressure on 
the liquid side of the liquid- vapor interface is 

= P - APt - AP, - AP (1) 

2 o L f V ' ' 


which is the supply pressure Pq decreased by the liquid momentum pressure 
drop APl,, the pipe friction pressure drop APf and the valve pressure drop 
APy. The latter two terms are given by 



( 2 ) 


where is the liquid velocity, pj^ the density, D the pipe diameter, X 2 the 
distance from the pipe entrance to the interface location and f the friction 
factor; and 

= 7 Px. [(A/A,)^ - 1 ] (3) 

where A is the pipe cross-section area and Ay the valve area, which varies 
with time according to a specified schedule. 


The liquid velocity is found from the fluid momentum equation for a constant 
density liquid and constant area pipe, which is integrated over a time interval 
from tj to tj^i with an average value of APl to give 






2 AP, 


n 1/2 


( 4 ) 
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which provides a time -stepped solution for the velocity when the other terms 
are known. 


When vaporization of the liquid is occurring at the interface, a momentum 
pressure drop results fronri the velocity change, giving 


P 


3 



W 

<“v 3 - "l 2 > 


( 5 ) 


where is the mass rate of vapor formation and 11^3 is the vapor velocity 
at the interface. The pressure ratio across the vapor stream is 


P3 

M4 

r 

1 + 

(V- l )/2 

“4' 

P4 

- M3 ] 

1 + 

(V- l )/2 

M^ 

4 


( 6 ) 


where the pipe exit pressure P4 is equal to the receiver vessel pressure P5 
for unchoked flow. To utilize the latter equation, the following expression for 
Fanno flow in the vapor relates the two Mach numbers; 



where f^ is the friction factor for the vapor flow. This equation is solved for 
M4 when M3 is known. The latter is found from the vapor velocity at the' 
interface given by 



where W^, the total mass rate of vapor production, is 


W = ZW 

V V . 

J 


( 9 ) 


where ZW'^. is the sum of the vapor production rates from boiling for individ- 
ual segments of the liquid from the heat transfer analysis. 
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The pipe and liquid are divided into small segments of equal length AX for the 
heat transfer analysis. The liquid interface travels from X 2 j to X 2 j+i across 
one of these segments during each time step At. Since the velocity Ul is a 
variable, the time step also varies during the solution. The temperature 
Tj+i of each nonboiling segment is calculated from its previous temperature 
Tj and a forced convection heat transfer rate 

%b = ^ - V ''(^) (^) ■ 

in the temperature equation 



T. 

J 


*^nb 


At 


C Dp 
P 


( 11 ) 


When a liquid segment reaches saturation temperature Tg-^^-, determined at 
P 2 , an empirical boiling heat flux q^^ is applied to the segment with the tem- 
perature remainiag at Tg^^j. The mass of fluid in a boiling segment will 
occupy a larger volume due to the production of vapor and will have a larger 
heat transfer area than a nonboiltng segment. Ignoring the loss of liquid in 
boiling, the new heat transfer area is 


A 

Vl 


W At TT D 

V. 



( 12 ) 


where Psv Is the density of saturated vapor. The vapor production rate for a 
segment is 


W 

V . 

J 


R- 

b J 


(A + A \ 

\ Vl ^j) 

2 Ah 


(13) 


where the average value of Aq is used, Ahy is the heat of vaporization and 
the effect of quality Y on two-phase boiling heat transfer is given by the factor 
R for which 




R = exp - 12. 4 (Y 


.1. 715 


- 0. 254) 


(14) 
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represents a fit of empirical data. A convenient expression for quality is 



( 15 ) 


which uses the approximation for the two-phase heat treinsfer area. 


While the heat transfer analysis treats the expanded volume and heat transfer 
area of the boiling liquid segments, tbe preceding fluid flow analysis assumes 
a viniform slug of liquid from the supply vessel to the liquid vapor interface 
with all vapor production occurring at the interface plane. 


The solution proceeds by determining the pressure distribution in the pipe 
for each time step, assuming a constant feed pressure Pq and iterating on 
APj_, to achieve a discharge pressure P4 equal to the receiver vessel pres- 
sure P5. (When the discharge flow is choked, the solution condition is 
M4 = 1 and P4 2; P5. ) .For. a chosen value of APl, Equation 4 gives the liquid 
velocity and Equations 1 ' through 3 give the interface liquid-side pressure 
With Tgg^^ determined at P2 and the time step determined, the heat transfer 
Equations 9 through 1 5 are evaluated giving Wy. Equations 8, 7 , 5 , and 6, 
in order, then give the values of M3, M4, P3, and P4, Based on the value 
of P4, the trial value of APj_, is corrected and the iteration repeated. 
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Appendix G 

CENTAUR SIMULATION MODEL 


\ 

f 



FLUID: 


1 

SR05 

TANK 




2 

SR26 

(WET) SCREEN 

LH2 


LO 2 

3 

SR26 

(DRY) SCREEN 

SAT P. = 19.17 PSIA 


29.38 PSIA 

4 

SR02 

RESISTA.TCE (SUBCOOLER) 

p =4.35 LB/FT^ 


68.8 LB/FT^ 

5 

SROl 

JUilCTIOd (SUMP) 

K = 11558 PSI 


108168 PSI 

6 

SR03 

BRANCH 

a = 1.23 X lO”'^ 

LB/FT 

8 X 10"^ LB/FT 

7 

SRll 

ODD VALVE 

VAPOR Y = 1 . 4 


1.4 

8 

SRll 

EVEil VALVE 

VAPOR RT = . 40 = 

30900 

• 180 = 8690 

9 

SR12 

ODD lilJECTOR 



10 

SR12 

EVEN INJECTOR 





TAiNK : LH2 LO 2 


VOL = 1265.4 FT^ 375.7 FT^ 

DIA = 120.0 liJ 120.0 Id 

HALL = .014 Id .018 IN 

ULLAGE P = 20.5 PSIA 31.5 PSIA 


Id9 


PRIBCEDINO 


BLAN& m 





SCREEN : 




MESH = 325 X 2300 
BPS = 2.29 FT 
AS! = 1.36 
BSl = 2.27 
AS2 = .13 
DPS = .00001 FT 
DS = 6.34 IN 
P0 = .06 
DC = .2967 
LS = 48.0 IN 
EE = .077 
STEM = 15.75 
F = 1.5 

PU = 1.0 or 0.0 


325 X 2300 
.942 FT 
1.13 
2.27 
.0534 
.00001 FT 

6.0 IN 
.06 
.2808 

4.1 IN 
.091 
15.75 

6.2 

1.0 or 0.0 


PLUMBING: 




SUMP 
COMMON 
ODD BRANCH 
EVEN BRANCH 


12-IN. DIA X .016 X 
33 IN. LONG 
3-IN. DIA X .016 X 
40 IN. LONG 

2.5-IN. DIA X .016 x 
62 IN. LONG 

2.5-IN. DIA x .016 x 
62 IN. LONG 


16.7-IN. DIA x 
12 IN. LONG 


2.5-IN. DIA X 
50 IN. LONG 

2.5-IN. DIA X 
57 IN. LONG 


SUBCOOLER RESISTANCE : Cy = 38.7 

AP = 1.2 PSI 


45.4 

1.4 PSI 


VALVE: 


Cy = 253 242 

OPEN/CLOSE TIME = .05 SEC .05 SEC 


INJECTOR: 


Cy = 61.2 


31.9 


.016 X 

.016 X 
0.16 X 
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